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PREFACE 


‘We have attempted in this book to: present the subject of partial differ- 
ential equations at a level which can be appreciated and used by the 
engineering student, Only an elementary knowledge of ordinary differ- 
ential equations is a prerequisite to ite understanding, The guiding 
light throughout the book ia the practical solution of partial differential 
equations, little attention being paid to questions of existence, uniqueness, 
and convergence. It is felt that students (upper clasemen or beginning 
graduate students) are not interested in mich considerations, and that 
‘these topics might obscure the fundamental techniques rather than add 
to the clarity of the exposition. Furthermore, since all illustrative prob- 
lems are based on physical systems, physical reasoning can be used to 
guarantee intuitionally the existence of a solution and make plausible 
the fact that the function we formally obtain is actually a true solution. 
In Chapter I we derive, in typical engineering fashion, some of the 
more common partial differential equations associated with such phenom- 
ena as vibration, heat flow, electricity, and elasticity. The second 
chapter is devoted to an exposition of fundamental facts concerning 
Fourier series. These ideas are, of course, useful in themselves; but 
their main value from our point of view is their use in solving partial 
differential equations. _In Chapter III we apply the techniques of Fourier 
analysis developed in Chapter II to the equations derived in Chapter I 
and solve a number of equations by the method of sepsration of variables. 
Certain problems, for example those involving nonperiodic functions, 
are not amenable to the method of Fourier series, hence in Chapter IV we 
extend the Fourier method to arrive at the Fourier integral, which enables 
us to solve a wider class of problems. The Fourier integral also leads to 
the transform method” of solving partial differential equations. In 
Chapter V we return to the method of separation of variables in conjunc- 
tion with changes in the coordinate systems, and are led to certain ordinary 
differential equations of practical and theoretical importance, namely, 
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Legendre's, Besscl’s, and Mathicu’s equations. A digression must be 
made to consider various elementary properties of the Legendre, Bessel, 
and Mathieu functions. Armed with this information we can solve, by 
the method of separation of variables, the partial differential equations 
which gave rise to these special ordinary differential equations. It seems 

iate in a book of this kind to consider certain elementary aspects 
of the general theory. Chapter VI gives the student some insight into 
the general structure of differential operators. 

Other methods of solution of partial differential equations, such as the 
Laplace transform and numerical procedures, have not been included in 
this book, both in order to limit its size and scope and because excellent 

Part of the material covered in this book has been the basis of courses 
‘taught by Professors M. G. Salvadori and R. J. Schwarz for the last ten 
years in the School of Engineering, Columbia University, and by the 
author in the College of Engineering, New York Univerity. I am in- 
debted to Professors Salvadori and Schwarz, who read the entire manus- 
cript and offered many valuable suggestions, and to my father, Mr. W. A. 
Miller, for the preparation of the many excellent diagrams. 


KS M. 
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Chapter I 


DERIVATION OF PARTIAL DIFFERENTIAL 
EQUATIONS 


1, Introduction 

Partial differential equations arise in many physical problems 
and, as in the case of ordinary differential equations [25],* a great 
many unrelated physical phenomena give rise to the same equation 
if the variables and constants are appropriately interpreted. In 
setting up the partial differential equations governing certain 
physical phenomena we shall restrict ourselves to elementary 
equations requiring little knowledge of the specialized fields of 
engineering. However, we shall assume that the reader has some 
familiarity with the fundamentals of mechanics, heat, electricity, 
and fluid mechanics. 

The list of equations derived in this book is by no means 
exhaustive : it is merely representative. At the end of this chapter 
a list of other important equations is given, with references to 
books containing their derivations. This category of equations 
requires a more specialized knowledge of the various branches of 
engineering. The equations considered in the text might be 
called the “ classical partial differential equations of mathematical 
physics.” They will be derived in typical engineering fashion, 
that is, with more appeal to physical and qualitative reasoning 
than to mathematical rigor. 

After the chapter on derivations (which may be omitted by any 
reader who is interested in the solution but not in the origin of the 
equations), we shall develop in Chapter II a powerful tool, Fourier 
series, which will enable us to solve some of the equations of 





‘* Numbers in square brackets refer to the references at the end of this book. 
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Chapter I. In Chapter III the methods of Fourier series will be 
applied to obtain explicit solutions for a wide variety of physical 
problems. Later chapters will deal with various extensions of the 
tools of Chapter II which are necessary to solve wider classes of 
partial differential equations than those amenable to solution by 
Fourier series. 

It will be necessary to recall a formula from the calculus (see 
[24] Section 5.9), which will be used many times in the various 
derivations of this chapter, and which can be easily derived from 
the Taylor series expansion of a function f() about the point z = a. 
‘The expansion is given by 


fe) = fle) + Faye —ay + F) ES 4 prey EGP + 4 


where f’(a) means the derivative of f(z) with respect to z evaluated 
at the point z =a, f(a) is the second derivative of f(z) with 
respect to 2 evaluated at the point s = a, etc. If, in the above 
equation we let 

a=ath ad a=, 
we obtain 


Sle +h) = fle) + fee + ree + ree er 


On the aseumption that A is small in absolute value compared 
with unity, the terms of the series involving high powers of k may 
be neglected, and we may write as a first approximation 

Set) = fis) + Sr. a 
Thus the value of f(x) at the point x +h (near to x) is given 
by the value of the function at the point x plus the value of its 
first derivative evaluated at x, multiplied by the distance 
[(e +4) — = = Aj between the two points. 

A typical example of the use of Equation (1) is the following. 
Suppose the value of the partial derivative of a function u(z,y) 
with respect to x at the point x, that is u/as, is known. Then the 
value of au/@x at the point x + dx is given by Equation (1) as. 


(3). = (@),+ &@)erea 
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oe 
‘Ou ou ou 
()... at ae 
‘The subscripts x and x + dx indicate that the value of the 
derivative of u(x, y) with respect to x must be evaluated at the 
points (x,y) and (x +'dx,y), respectively. However, the sub- 
script x is generally omitted, our use of it in the above example 
being merely for emphasis. 


2. The vibrating string 
As a fit example of a partial differential equation we shall 
derive the equation satisfied by the small transverse displacements 
of a vibrating string. ‘The sketch in Figure | represents a tightly 
, 


stretched string of length L, fixed at the end points, and subjected 
to a constant tension of T pounds. Our problem is to determine 
what equation governs the motion of the string under various 
initial conditions and to find a function u(x,t) which gives the 
deflection of the string at any point x and at any time f. 

In order to obtain a simple differential equation we shall consider 
a more or less “idealized” string. Presumably, any physical 
string will reasonably well approximate the “ideal string,” and 
hence the function u(x,t) obtained from the simplified equation 
will reasonably well represent the actual displacement. We make 
the following simplifying assumptions. 

1. The string is perfectly flexible, that is, it cannot resist 
bending moments. 
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2. The deflections u are small compared with the length L of 
the string. 

3. The slope at any point of the deflected string is small 
compared with unity. 

4. The tension T is constant at all times and at all points of the 
deflected string, and is large compared with the weight of the 
string. 

5. The horizontal displacement of the particles of the string 
is negligible compared with the vertical displacement, that is, we 
have pure transverse vibrations. 

6. The motion takes place only in the xy-plane. 

If we check through the above six assumptions we see that they 
are all reasonable and consistent. By this we mean that no contra- 
dictions will arise from the acceptance of the truth of all the above 
assumptions. (An example of a pair of inconsistent assumptions 
would be to assume that the displacements u are large compared 
with the length of the string, but that the horizontal displacements. 
are small.) 

Consider a small differential segment ds of the string (see 
Figure 2), and let w be the weight per unit length of the string. 
From elementary mechanics, 2 F = ma; in words, the sum of the 


y 
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forces in any direction applied to the string segment is equal to the 
product of the mass of the segment by the acceleration in that 
direction. In particular, the sum of the forces in the y direction 
is 
BF, = —Tsina + TsinB — wds, 
Since the deflections are small, we may replace sin « and sin 8 by 
tana and tanf, respectively. Also, by the third assumption 
dy <de and ds = dx to first-order approximation. Hence the 
sum of the forces in the y-direction may be written in the form 
ZF, = —Ttana + Ttanf — wdx. 


‘The slope at any point x of the deflected string is, from the calculus, 
ou 
tana = rg 

and by Equation (1), tan f is given to first-order approximation by 


ou Oo (tu ou eu 
anp= 24 Z(H) a= B+ Fede 


Hence & F, becomes 


ou au am 
Br = TH+ Th + 1H ae — wee 


(2) 
Oe 
= Trae — wdx. 


‘This sum of forces in the y-direction must be equated to the 
product of the mass m and acceleration a, in the y-direction. The 
acceleration a, in the y-direction is the second derivative of 
distance with respect to time, namely, 

ou 

Po 


age 
and hence 

° 

ma, = Sart, @ 


where g is the gravitational constant. By Newton's second law 
we obtain, from Equations (2) and (3), 
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or, dividing by dx and multiplying by g/t, 
sT Ou _ aw 
fo ot E> Oe ® 

‘The constant g7/w has the dimensions of velocity squared, and it 

is customary to write 


where a has the dimensions of velocity. Also, by the fourth 
assumption, g may be neglected in comparison with ou/at*, and 
Equation (4) assumes the standard form 

amt = *. (5) 
This is the partial differential equation satisfed by the vertical 
motion s(x.) of the vibrating string. It is also called the one- 
dimensional wave equation. 


3. Boundary conditions 

In solving a physical problem by means of differential equations 
we are frequently interested in the solution of the differential 
equation which satisfies a given set of initial or boundary con- 
ditions. ‘The boundary conditions for the case of the vibrating 
string are readily determined by inspection as 


1 9 = 0, 
2 (L,t) = 

a ) 
3 2am) = 
4. w(x,0) = fiz). 


a 
* 5; 4s, 0) means the derivative of ws, £) with respect to# evaluated at t= 0. 
‘This notation is consistent with that used for ordinary derivatives (see Section 1) 
‘where we write (a) to indicate the derivative of f(s) with respect to = evaluated 
st = 4, Other more or less standard notations for ux, 0) are 


fuen + ud |g 6 uke 0). 
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‘The first two boundary conditions state that the string is fixed at 
the end points, x = 0 and x = L for all time ¢. The third states 
that initially the transverse velocity of the string is a given function 
of x, namely g(x), while the fourth indicates that the original shape 
of the displaced string (that is at time ¢ = 0) is another specified 
function f(x). [If the first two boundary conditions are to be 
satisfied, f(0) = f(L) = 0.] 

‘A question arising in connection with Equation (6) is whether 
this is the requisite number of boundary conditions, that is, 
whether one or many or perhaps no functions exist which satisfy 
the differential equation and these four boundary conditions. By 
physical reasoning it seems reasonable that in this example, at 
least, there exists a unique solution, but in general no conclusive 
answer can be given to the question of whether the boundary 
conditions determine a unique solution, many solutions, or no 
solution at all. Certain remarks, however, may be made in this 
connection. First, we cannot get new boundary conditions by 
differentiating known boundary conditions. For example, 

Zuo =f) or Ze us0) = ete) 
are not additional boundary conditions. Second, if the initial 
position f(x) and velocity g(x) of the particles of the string have 
been specified, the acceleration, say 





oles 0) = He), 


and higher time derivatives cannot be specified. For, once we 
have specified u(x, 0) and 


a 
2 as.0), 


the higher time derivatives are uniquely determined by the dif 
ferential equation. This situation is readily illustrated in the case 
of the elementary ordinary differential equation #+kx = 0. Here, 
once the initial position and velocity, say (0) = a, (0) = b are 
specified, we obtain 2(0) — — ka, 3(0) — —hb,ete. from the dif- 
ferential equation itself. If we were to specify (0), #(0), and 
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(0) independently, we would have no solution to the problem, 
unless (0) were fortuitously chosen as — kx(0). 

A differential equation together with its boundary conditions is 
said to form a differential sy:tem, and the modern theory of partial 
differential equations is concerned with differential systems and 
not with differential equations alone. In the classical theory of 
ordinary differential equations it is frequently not too difficult 
task to find the general solution of the equation involving various 
arbitrary constants and then determine these constants so a8 to 
satisfy the boundary conditions. The equivalent problem in 
partial differential equations is much more difficult, and we must 
look at each differential system as a new problem even though the 
differential equation part of different systems may be identical. 
Considering differential systems rather than differential equations 
alone is not too great a restriction, for if the solution is to have 
physical meaning it must satisfy some boundary conditions. 

To summarize then, our problems will be concerned with 
differential systems, or as we shall say, with boundary value pro- 
blems, and we shall solve particular partial differential equations 
subject to certain specific boundary conditions. 


4. The vibrating membrane 

As a second example we shall derive the equation satisfied by the 
transverse motion of a vibrating drumhead or membrane. Such 
a two-dimensional membrane is illustrated in Figure 3.* Our 
problem ia to find the equation satisfied by the deflection u(x,y,1) 
of a thin membrane oriented in the xy-plane. As in the case of 
the vibrating string certain “idealizing” assumptions will be 
made. 

1. The force per unit length, S, is constant and everywhere 
normal to the boundary of the membrane. 

2. The total force on the boundary is large compared with the 
weight of the membrane. 

‘* ‘The cross in the circle indicates that the x-axis is into the paper, while « dot 
in the circle (see, for example, Figure 4) indicates that the axis is out of the plane 
of the paper. 
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3. The membrane is so thin that it cannot resist bending 
moments. 

4. The deflections u are small compared with the lateral 
dimensions of the membrane. 

5. The slopes are small compared with unity 

6. The lateral displacements are negligible compared with the 
vertical displacements. 


& 





ys to) 


y 
Fig. 3 


Consider a differential element of the membrane of area dx dy 
as indicated in Figure 3. If at a fixed time ¢ a plane parallel to the 
xz-plane is passed through this element, the cross section of the 
element dx dy will appear as in Figure 4. Similarly a plane passed 
parallel to the yz-plane will yield a section as in Figure 5. If w is 
the weight per unit area of the membrane and S the constant 
normal force per unit length, the sum of the forces (E F,) in the 
a-direction is given by 

LF, = —Sdytana + Sdytanf — Sdrtany 
4+ Sdetan3 4 w dx dy, ” 
since by virtue of the fourth assumption, slopes are small and 
sin a, sin, sin y, sind are approximately equal to tana, tan, 
tan y, and tan 8, respectively. 
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From the calculus, 


tana = 3 and tony = 


a 
3 





Fig. 4 
and from Equation (I), 
tnp = 2+ 2 (2) as, 
Ou a (eu 
b= sla)e- 
moa 5+ (a) % 





Fig. 5 


Q 
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Using these results, Equation (7) becomes 
Ou (ou Pu 
EF, = —Siy + Sy (S4 Hac) 
-Ss ax Sdx 
Bt S(t ee 


Simplifying, and noting that w dx dy is negligible by the second 
assumption, 5 F, assumes the form 


) + wdedy. 


BR = Saray (Fs+ 3): 


at OF 
From the fundamental law of mechanics, EF, = ma,, where a, is 
the acceleration in the z-direction. The mass of the elementary 
rectangle dx dy is 


o 
m= aedy, 
: dy, 
and its acceleration is 
Ou 
a = 


Hence 


ou 
S de dy (% zat +%)-4 © ae dy = 
Dividing by dx dy and multiplying rsh by g/t, this equation 


becomes 
gS (Be OY Oe 
wo lat + ae) ~ ae 


‘As in the case of the vibrating string, gS/wo has the dimensions of 
velocity squared, gS/w —a', where @ has the dimensions of 
velocity. The equation of the vibrating membrane then assumes 
the claseal form 

o(% + =) — @) 
and is often called the wave equation in two dimensions. 

In order to determine a differential system we must associate 
boundary conditions with Equation (8). Typical boundary con- 
ditions might give the initial displacement and velocity of the 
membrane and state that the boundary is held fixed for all time. 
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These conditions can be expressed mathematically as 

1. fs, flz),t) = 0 

2 Zu(ay.0) = Hey) 0 

3. u(,7,0) = $y) 
where y= f(z) is the equation (in the 2y-plane) of the curve 
bounding the membrane (see Figure 3). Condition one states that 
the boundary of the membrane is held fixed and is analogous to the 
first two boundary conditions of the vibrating string problem. 
The second boundary condition asserts that the initial velocity 
of the membrane is a prescribed function of x and y, while the 
third boundary condition gives the initial displacement of the 
membrane. Note that (x,y) is not quite arbitrary, for if the first 
boundary condition is to be satisfied, 4(x,y) must vanish on the 
boundary. 

Our results can be extended to the case of a three-dimensional 
vibrating elastic body. In this case the differential system assumes 
the form as the «te a 

a (toe ge) = (10) 
with appropriate boundary conditions (see [17]). Equation (10) 
is also called the wave equation in three dimensions. 


5. Longitudinal vibrations of a bar 

Another example of a physical phenomenon giving rise to a 
partial differential equation is that of the longitudinal vibrations of 
a bar or of a long helical spring. The simplifying assumptions 
usually made regarding the bar are : 

1. The material is homogeneous, isotropic, and elastic; plane 
sections remain plane. 

2. The material satisfies Hooke’s law : o = Ee, where o is the 
stress (psi), E is Young's modulus (psi), and ¢ is the unit strain 
(inches/inch). 

3. The displacements are small in comparison with the 
dimensions of the bar. 


Google 


Ar. 5 DERIVATION OF PARTIAL DIFFERENTIAL EQUATIONS BEY 


Referring to Figure 6, we let 
dz — the length of a differential element of the bar, 
displacement of 2 given section of the bar at a given time f, 
— unieatrain — (new length) — (old length) 
¢ = unitstrain = (oer Tenge) — iene) 
u 






Fig. 6 


AB 
(@+ de) + (wt du] — (e+u) = det du, 
old length = AB = dx. 

Hence ¢ = Meta) — de _ 

oe 

Let A(x) be the cross-sectional area of the bar at x. This is shown 
to be a rectangle in Figure 6, but could be any shape. Let F, be 
the total force in pounds at the point x, assumed positive when 
tensile (see Figure 7). The force F, is equal to the stress o multi- 
plied by the cross-sectional area, 
= 0A, 
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From our second assumption ¢ = Ee, and with « = 6u/ax, 
ou 
Fy = 0A = Bed = BAX ay 
Hence the force Fs» at the point x + dx is, by Equation (1), 
ou é ou 
Farvae = EAR + 3 (BAP) de (12) 


Therefore by the fundamental law of dynamics, namely the sum 
of the forces in a given direction is equal to the mass times the 
acceleration in that direction, we obtain 
o au 

Five — Fe = ($44) - (13) 
where w is the weight per unit volume of the bar and conse- 
quently (zo/g)A dx is the mass of the element shown in Figure 7; 
and @u/at? is the acceleration. From Equations (11) and (12), 


Fetas — Fe = &(# 3) ax 


and from Equation (13), 


a Ou w, hu 
z (24%) de = 2A ae 
Dividing by dx and multiplying by g/d, the equation assumes 
the form a 
£2 (pg). &. 
wA Bx (za 2) = oe aay 


which is the equation governing the longitudinal vibrations of the 
bar. If we make the additional assumption that EA is constant, 


@ ou vg Ou 
Z («4 2) = 1S, 
and Equation (14) takes the simple form 
eB uo 
wast ~ OF 
The constant at = gE} has the dimensions of velocity squared, 
and by its use Equation (14) now assumes the classical form 


Ou 
am (15) 
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which we recognize as the one-dimensional wave equation. We 
see that Equation (15) is identical with that of the vibrating string 
derived in Section 2. Here, then, is our first example of two 
different physical phenomena giving rise to the same differential 
equation. 


6. Transverse vibrations of a beam 


A beam of length L loaded by a variable load of eo(x,t) lb/ft is 
sketched in Figure 8 and referred to a y-axis positive downward 
and an x-axis positive to the right. We wish to find the transverse 
vibrations (x,t) of this beam due to a set of initial conditions. 





Fig. 8 


The weight of the beam itself will be assumed to be included in 
the load (x,t). The “ idealizing” assumptions made to set up 
this equation are those of simple beam theory (that is, the assump- 
tions of the classical Bernoulli-Euler theory), from which it may 
be shown that (see [25] Section 2.10) 


M = —EIC, (16) 


where M is the bending moment (ft-lb), £ is Young's modulus 
(psf), is the moment of inertia of a transverse section of the 
beam (ft), and C is the curvature of the beam (ft). Since from 
the calculus ({24], Section 2.7), the curvature (reciprocal of the 
radius of curvature) is given in Cartesian coordinates by 


y" 


cm Tey 
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and by the small slope assumption of elementary beam theory, 
x <h 
we have to first-order approximation, 
c= 





Equation (16) then assumes the form 
= 227 
M=-HZ. ay 


The quantity E/ is known as the “ flexural rigidity” of the beam. 
‘The equations of motion of an element dx of the beam are given 

by Newton's second law : 

ZF, = may 

EM = Ja. 
The first states that the sum of the forces in the y-direction is 
equal to the mass times the acceleration in the same direction, 
while the second states that the sum of the external moments about 
any point is equal to the moment of inertia of the element times 
its angular acceleration. Since we are considering pure transverse 
motion, the angular acceleration is assumed to be zero, and the 
above equations reduce to 

iF, = my XM =0. 

Let us apply these equations to the element dé illustrated in 
Figure 9, in which, with the coordinate axes of Figure 8, shears Q 
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and moments M are positive as indicated. The sum of the 
forces in the y-direction, © F,, is 
~. 2 a 82: 
BR = 04 (0+ Lax) + wae = Bae + was, 


‘The mass of the element is (w/g) dx and the acceleration in the 
y-direction is ey/at’, Hence 


Beet wie = Bae 


or, dividing by dz, S 

2 
+0= FR. (18) 
In order to eliminate the shear Q and obtain an equation involving 
only y and its derivatives, we utilize the moment equation Z M = 0. 
‘Taking moments about the point x of Figure 9, 


Oe 


em nam — (M+ Mae) + (9+ Bar) ee + oanF 


= Mat gar + Ras + of. 

be 
But terms involving de® arc infinitesimals of higher order ({24], 
Section 1.6) with respect to terms involving dx, and hence can be 
neglected. Therefore 


am 
EM = —Z dx + Ode = 0, (19) 


or canceling dx, 
a 2 Om (20) 

Differentiating Equation (20) with respect to x, 

2Q _ aM 

os ~ Oe? 
and substituting this expression in Equation (18), we obtain an 
‘equation involving only y and the moment M, namely, 

aM wo aly 


ee tem Fae (21) 
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Finally, differentiating Equation (17) twice with respect to x, 


ae = - as Fa (2133), 


and substituting in Equation (21), we obtain an equation involving 
only y: oa fis 7 
© 
-m (#2) tons 
‘This is the partial differential equation satisfied by the transverse 
motions of a beam. If the load w(x,t) can be neglected with 
Fespect to the inertia forces (w/g)(d*y/at*), the above equation 
assumes the form 
way | Bay 
go Oe (# ) = mae 
If, further, the flexural rigidity EJ is constant, we obtain the 
additional simplification, 
2 ty oe 
Fo t Han = 0 
or, on letting a? = Elg/w, 


wad O. (22) 


The quantity a? is not a constant unless w(x,t) is constant. We 
leave as an exercise for the reader the problem of determining 
the dimensions of a and of verifying the dimensional correctness of 
Equation (22). 

Common sets of boundary conditions associated with Equa- 
tion (22) state that the beam is “ simply supported ” (Figure 10) 





Fig. 10 


Google i 


Art. 6 DERIVATION OF PARTIAL DIFFERENTIAL EQUATIONS oe) 


or “ built-in” at both ends (Figure 11), When the beam is 
“ simply supported” these conditions are 


L 04) = 0, 
2 Lt) = 0, 
3. pa ¥(0,t) — 0, 


4 Zrono @ 


5. x20) = fle), 
6 20) = ee). 





Fig. 11 


‘The first two boundary conditions state that the end points 0 and L 
are fixed for all time t. ‘The third and fourth conditions result 
from the fact that a simply supported beam has zero moment and 
hence zero curvature at its free ends. The fifth gives the initial 
displacement of the beam and the sixth the initial velocity. 


The boundary conditions for the built-in beam are 
1 ¥0) = 0, 
2. Ls) = 0, 


a 
3. £0) = 0, 
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a 
4 EHL) = 0, 
5. 9x0) = fla), 
a 
6 F3K20) = als). 


(24) 


As before, the first two conditions state that the end points 0 and L 
are fixed, while the third and fourth result from the fact that the 
slope of the deflection of a built-in beam is zero at the end points. 
The fifth and sixth give the transverse displacement and 
velocity of each particle, respectively. 





1. Heat flow in one dimension 


The problems that appear in the next few sections will involve 
the flow of heat. These equations are of interest to the chemical 
engineer and physicist as well as to the applied mathematician. 
Our first application will be to the flow of heat in one dimension. 
Suppose we have a long thin bar or wire oriented along the 
x-axis, and we wish to determine the temperature distribution 
u(x,t) in the bar. We shall assume that the bar is perfectly 
insulated laterally, and that heat flows in one dimension only 
(the x-direction). ‘The following physical laws will be used : 


1. The heat contained in an element of mass m is proportional 
to the mass and to the absolute temperature, that is, 


Q = om, 25) 
heat (calories or Btu), 
specific heat, 





where, Q 
c 
m = mass, 
& — absolute temperature. 
2. The rate of heat flow through a surface is proportional to 
the area of the surface and to the outward normal gradient of 
temperature, that is, 


ie = — ka & calfec or Buujsee, (26) 
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where, A = cross-sectional area of surface in a direction 
perpendicular to the flow of heat, 
éu/én = gradient of temperature perpendicular to the 
surface, 


hk = thermal conductivity. 


‘The negative sign in Equation (26) follows from the fact that heat 
flows into the body if the temperature to the right of the boundary 
is lower than the temperature to the left, that is, if 2u/an is negative. 





Fig. 12 


Figure 12 represents the insulated bar. If we consider a small 
segment dx of the bar, the amount of heat Q, in this elementary 
element is, from Equation (25), 

Q, = (Ade) du, 
where 8 is the density (mass per unit volume) of the material of the 
bar, and hence A de 8 is its mass. The rate at which heat is stored 
in the element of volume A dx is therefore 2Q,/2¢, 





= cb ae 2 


ot 


If we assume that the positive direction of heat flow is the x-axis 
we find from Equation (26) that the rate at which heat flows into 
the element is 


2Q, ou 
Ba -Ms a) 


and the rate at which heat is flowing out of the element is, by 
Equation (1), 


Qs, 2 (2, ou 8 (, aw 
Bee (Bt) ae = na — 2 (ea ee) 
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The rate of heat flowing in, minus the rate of heat flowing out, 
is the rate at which heat is stored in the elementary volume A dx. 
Hence 





or, by Equations (27), (28), and (29), 


ou é 
Ab ds = x (4 





Dividing by de, we obtain the equation 
ou é Ou 
cape = 2 (4% 5 
If we assume k and A to be independent of x, we may also write 


Au Bou 


ae he @) 


which is the standard equation for the flow of heat in one dimension. 


8. The three-dimensional heat equation 


The above results are readily generalized to the flow of heat in 
a three-dimensional body. Figure 13 represents an elementary 


norma! 





(xy,2) kz 


Fig. 13 


parallelepiped dx dy dz in the three-dimensional body. If we 
let u(x,y,z,t) be the temperature at any point x,y,z and at any 
time t in the volume, then by the fundamental law of Equation (25), 


22, 
et 





dx dy da) 3 (Gl) 
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is the rate of heat stored in the elementary volume, where 8 is the 
density of the material, From Equation (26),.the rate at which 
heat flows into the solid through the three faces at x,y,z, normal 
to the coordinate axes, is given by 

2Q, _ u ou au 

fe OREO HNO a hee, (32) 
and the rate at which heat flows out through the three faces at 
x +dx, y +dy, 2 -+dzis given by 


20 eu, Ou 
oe — haedy (24 Has) — kay as ( 


a 
au Ou ) 
— kdzdx|—+— 5 dy}, 
te (4 
where we have assumed that the thermal conductivity & is inde- 
pendent of x,y, and z (and hence can be taken outside the 
differentiation signs). As in Section 7, the rate of heat that flows 
to the body, aQ,/ét, less the rate of heat that flows out of the 
elementary volume, 2Q,/2t, is the rate of accumulation of heat 
@Q,/at in the body. Expressing this fact in formulas gives 


221 _ Qs _ 20s 
ot at 












Substituting Equations (31), (32), and (33) in the above equation, 





eu Au Au au 
8 dx dy de = (% ++ Be) as dy ds 
or 
eu au eu 8 Ou 
mw ty tm ae (34) 


which is the standard form of the heat equation in three dimensions. 


The differential operator 
a a é 
ast + a t oe 
which was met before in the case of the vibrating elastic solid, is 
called the three-dimensional Laplacian, and is symbolized by V*. 
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Using this standard notation, the heat flow equation may be 
written as 


ou 
We = i 
If the temperature u is independent of time, éu/at is zero, and 
Ve = 


For the case of heat flow in two dimensions, for example, in a 
thin flat slab, insulated on both faces where we assume the flow 
takes place only in the ay-directions, the equation assumes the 
form 


au, Ou Sw 
oe t oF Rat @5) 


This equation is also written 
& 
vu = oR 


Be 


where here V? is the two-dimensional Laplacian. Whether V? 
eee yf. # a. @ 
a 
Biptm ' Bte 


will be clear from the context. 


9. The transmission line equations 


‘As an example of partial differential equations which arise in 
clectrical engineering we shall derive the transmission line or 
telegraph equations. Suppose we have a transmission line which 
we shall consider a distributed constant-parameter system. By 
this we mean that the resistance, inductance, and capacitance vary 
linearly with the space coordinate x, that is, any small section de 
of the line will contain the same amount of resistance, inductance, 
and capacitance a8 any other element of length dx, (in contra- 
inction to a lumped parameter system in which we assume that 
we have a finite number of resistors, coils, and capacitors con- 
nected by impedanceless wires). In a distributed parameter 
system the voltage and current will vary along the length of the 
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line, and in general, will also vary as a function of time. The 
problem is to determine the equations satisfied by the current, 
i(x,t), and voltage, e(x,t), distributions in the line. Referring to 
Figure 14, let 

R = series resistance in ohms per loop mile, 

series inductance in henrys per loop mile, 

shunt capacitance in farads per mile, 

= shunt conductance in mhos per mile, 





L 
Cc 
G 





x 4dr 
Fig. 14 


According to Kirchhoff’s first law, the drop in potential around 
a closed loop equals the rise in potential. Considering an elemen- 
tary section dz, the rise in potential at the point x is e(s,¢) while 
the rise in potential at the point x + dx is 


(x,t) + eax 
Hence the net rise around the loop in the clockwise direction is 
a ae 
(x,t) — [on +a] =-_e 
The drop in potential around the loop is equal to the sum of the 
resistive drops and inductive drops. The drop in the resistance 


Rar is 
iR ds, 
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and the drop in potential due to the inductance L dx is 


a 
Lat 

Hence 

ee “ a 

Rede = Ree + LS de, 

and, on canceling dx, 

ae ae (36a) 

ox oe 


The current i(x,t) also changes from the point x to the point 
x +x since some of the current is diverted through the shunt 
conductance and capacitance. By Kirchhoff’s second law, the 
sum of the currents flowing toward the node at x equals the 
sum of the currents flowing away from the node. The current 
flowing down through the capacitance is 


ae 
(Cas), 
and the current flowing through the conductance is (G dx)e. 
Hence 
4 ee ss a 
iat) = (C dx) + (Gale + Hut) + Fide, 
or 
a Ge 
-Ae = Gedx + Cy 


and on dividing through by dx, 
ai ae 
RT Gt CH (366) 
Equations (362) and (368) form a pair of partial differential 
equations. Both equations involve both unknown functions 
ext) and (x,t). It is a simple matter in this case to climinate ¢ 
(or ¢) and obtain an equation on ¢ (or) alone. This is effected by 

differentiating Equation (36a) with respect to x, 

@e 
~ ot 


a a 
= Rath a 
ow 


ée 
= — RGe- ROS +L G7) 
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(on substituting Equation (368) for i/ax) and differentiating Equa- 
tion (368) with respect to ¢, 


Substituting this result in Equation (37), 


Be ée ae 
LO Ga + (LG + RC) + RGe — FS. (382) 
In a similar fashion we can eliminate ¢ and obtain 
a i 
CL Ge + (RC +16) 5 + GR = 5 (388) 





Special forms of the transmission line equations can be obtained 
by making various assumptions regarding the parameters G, R, L, 
and C. If we assume that the resistance and conductance are 
negligible, which is reasonable at high frequencies, then 
Equations (38a) and (38b) reduce to 
ae a ai 
we Cb on = ae 
which are called the high-frequency line equations, and are of the 
same form as the one-dimensional wave equation. 

If we assume that Z and C are negligible, which is reasonable 
in the steady state at very low frequencies, we obtain the case of a 
direct current transmission ine, namely, 

Oe 
oe 


(9) 





= RGe, z = GR. (a0) 


Note that these are now ordinary differential equations. 





Finally, for a submarine cable we may assume L = 0 = G since 
the leakage conductance is very small and the frequencies are low 
enough to make the series inductance negligible. Thus we 


obtain the submarine cable equations, 
ae _ a a 
ROR = 5a RCE ae (a1) 


which are of the same form as the heat equation in one dimension. 
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10. Potential 


Another important equation of mathematical physics may be 
derived by considering the potential of the forces of attraction of 
two masses. If M and m are two point masses r units apart (see 
Figure 15), Newton’s law states that these point masses attract 

” 
Tab 






Foams OF eo 


(uy) 
Fig. 15 


each other with a force f proportional to the product of their 
masses and inversely proportional to the square of the distance 
between them, 


mM 
foo", 


where G is the constant of proportionality. In particular, if m is 
a uni¢ mass, 
GM 


fa 

Consider the function u — GMjr, called the force function, or 
ite negative — GM/r, called the potential function. The derivative 
of the potential function with respect to r is precisely the force f. 
We shall show that the potential function satisfies a certain partial 
differential equation. 

Instead of considering a single mass M we shall consider a 
number of masses My, My, ..., M, and hence a force function 





(42) 


where, if the coordinates of the Ath mass are (4y,,¢,) and the 
coordinates of the unit mass are (x,y,2), 


re VERA GOT WEF. 
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Equation (42) can be generalized. In the case of a continuous 
distribution of masses in some volume V, 


ueya) = 6 [ff 2%? ay (8) 
, 


where (a,b,c) is the density of the mass distribution, dV is an 
element of volume, and 


= VE=P FORT EH 


The partial differential equation satisfied by u(x,y,2) is derived 
from Equation (43) as follows. Differentiating u with respect to x, 
we have 








x 2G. i) plade) & =O) av, 


me ~6 fff pla,b.e) [7-325 |e. 


If we similarly calculate 2/292 and 2%u/2e® and add, 
vu = —6 fff paso (3-3) av =0. 
7 


Hence the force function u(x,y,2) (or its negative, the potential 
function) satisfies the three-dimensional Laplacian, 


Viu = 0. (44) 








In the theory of electrostatics, Coulomb’s law states that the 
force attracting or repelling two small charged bodies is 


r=Kk®, 
F 


where g and Q are the charges in coulombs on each body, respect- 
ively, r is the distance between them, and K is a constant of 
proportionality whose value depends on the medium in which the 
force is exerted. Since this law is of the same form as Newton’s 
law for masses, the electrostatic potential $ must also satisfy the 
three-dimensional Laplacian, 


Vig = 0. 


Gooal 


Q 
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The same law applies to magnetostatics. Two magnetic poles 
of strength m and M, respectively, attract (or repel) each other 
with a force directly proportional to the magnitude of their 
magnetic charges and inversely proportional to the square of the 
distance r between them, ti 


rane 


‘The constant of proportionality & depends on the medium in which 
the magnetic charges are embedded. Hence the magnetic 
potential yp also satisfies 

Vy =0. 


11, The equation of continuity 


‘The so-called equation of continuity, which essentially states the 
conservation of mass, plays an important role in the theory of 
hydrodynamics. If a fluid (liquid or gas) is enclosed in some 
container V, a vectorial velocity can be associated with each point 
in the fluid with components u, v, w in the x-, y-, and z-directions, 
respectively. The equation of continuity is satisfied by these 
‘components. 

lay + dyz4d2) 


(xyz) 


tx+dxy,2) 

Fig. 16 
Consider a small volume dx dy dz of the fluid (Figure 16) and 
let p = p(x,y,2,t)(b d fluidyat any point (x,y,z) 


at any time t. ‘The mass crossing the surface dy dz (perpendicular 
to dx at the point x) per unit time is 
u(p dy dz), 
while the mass crossing the surface dy dz at the point x + dx in 
the x-direction is 
pu dy dx + £00 ae dy ds. 
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Hence the net gain in mass per unit time from the flow in the 
x-direction is i 
— Fy los) de dy dz. 


Similarly the net gain in mass due to flow in the y- and z-directions 
is 


— Flo) ds dy de and — Zon) de dy ds, 


respectively. ‘The total accumulation of mass per unit of time is 
the sum of these three contributions, namely, 


@ @ a 
= [ZO + Foo + Zoom] aeayes, 45) 
But the total mass in the elementary volume is 
pdx dy da, 
and the rate change of this mass is 
& 
ds dy de. 


Equating this rate change of mass to that given by Equation (45) 
and canceling dx dy dz, we obtain the equation of continuity, 


2+ Zon + Zo + om = 6) 





If pis a function of t alone and not of x, y, 2, we obtain the simpler 


equation 
do | bw 


ep (ou 

a + Pat yt Ge. 
On the other hand, if p is a function of x, y, z and not a function of t, 
Equation (46) becomes 





HoH) + Slee) + Flom) — 0. 7) 


Finally, if p is a constant, the equation of continuity takes the 
special form 
au ao 


a0 
at y 


=o (48) 
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Equations (46), (47), and (48) involve three unknown variables, 
viz., u, v, and w. These equations are used in conjunction with 
the equations of motion of a fluid (Section 14), which give us three 
additional equations for u, 0, w, and the pressure p. Hence we 
see that the equations of hydrodynamics form a system of differ- 
ential equations. Theoretically v, w, and p could be eliminated 
to obtain an equation on u alone. This situation has been met 
before in the case of the transmission line equations. 


12. Tidal waves in an incompressible fiuid 


‘The motion of an incompressible fluid in a long channel furnishes 
another example of an equation in fluid mechanics. The deriv- 
ation of this equation involves some new techniques not 
encountered in previous examples—yet leads to the familiar wave 
equation in one dimension. 























rot 
}——— x —+| urdu |-4 


° x K+dK x xed 
Fig. 17 





A long channel of width h is oriented along the x-axis, Figure 17, 
and is filled to a height a by an incompressible fluid (for example, 
water). The fluid is disturbed by waves (as indicated by the 
dotted line in Figure 17) whose wavelength is assumed very long 
in comparison with the depth a of the fluid in the channel. Let 
u(sy) and 2(x,y) denote the components of displacement of a 
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point (z,y) in the x- and y-directions, respectively. We wish to 
derive the partial differential equation satisfied by the displace- 
ment x. Due to the wave disturbance an element of volume 


hds dy 

at (x,y) will be transformed into the element 
hx dy’ 

at (2',y’) (see Figure 17). Now 





vaesu 
and 
eu 
Std = etde put du = xt de tuts dx, 
Hence 7 
dx = de = Mas. 
a 
Similarly, 
+ ov 
Yaya ee 
and 


de! dy! = deay(1 





On the assumption of “long” waves the slopes du/éx and év/2y 
will be small, and hence to first-order approximation (neglecting 
second-order terms), 


dx’ dy’ = dedy (1+ 4 






a (49) 


But in an incompressible fluid the volume of an element cannot 


change, so that 
he de’ dy’ 





h ds dy. 


We conclude, therefore, from Equation (49) that 


au we 
ae + & = 


From the assumption that the waves are very long compared with 
the depth of the fluid in the channel we are justified in assuming 
that the horizontal displacement is independent of the depth y. 





(50) 
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Hence, by integrating Equation (50), 
va -yHs4, co) 
where A is a constant of integration. But at y = 0, » =O since 
there can be no vertical displacement at the bottom of the channel. 
Hence A = 0. Let o* be the component of vertical displacement 
at the surface (see Figure 17). ‘Then, from Equation (51), 
ot = a (52) 
‘The pressure on the surface A dy at the point x + dz is greater 
than that on the surface A dy at the point x since the former surface 
is at a greater depth. ‘This pressure excess dp is 
dp = 8g de®, 
where 5 is the density of the fluid and g is the gravitational constant. 
But 
eres oe 
aut = Fas. 
o* 
dp = bg Ze 
The force in the x-direction on the surface h dy at the point x 
exceeds that on the surface h dy at x + dx by the amount 


ae 
— bg Sas (hdy), 
From the fundamental law of mechanics, EF, = ma,, 
. 
(Gh de dy) OH = — 8p Bh de dy 
oF fu aut 
we hae 
Substituting @0*/2x from Equation (52) in the above equation we 
obtain the desired result, viz., 
Fe gg Mt 
af — 8° G8 
which we immediately recognize as the wave equation in one 
dimension. 


Hence 


(53) 


4) 
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13, The equations of elasticity 

Numerous partial differential equations appear in the theory of 
elasticity, but their derivation requires a specialized knowledge of 
the subject. We shall limit our treatment here to the derivation 
of one type of equation which can be deduced almost entirely from 
the principles of elementary mechanics. Other equations of 
elasticity will be found in Section 14. 


% oy 


Fig. 18 


Consider a thin elastic plate upon which no body forces are 
acting. On each side of an elementary rectangle, dx dy (Figure 18) 
normal components of stress ¢ and tangential components of 
shear 7 will act. Since the element is assumed to be in equili- 
brium, the sum of the forces in the x- and y-directions, 2 F, and 
=F,, must be zero. Hence 


BE, = ody + (04+ Be te) dy — ay de 


) a =0 





or 


(55) 
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Similarly, the equation EF, = 0 leads to 
Boy, Otay 
a + ee 
Equations (55) and (56) are known as the equations of equilibrium. 
Equating to zero the sum of the moments © M, about the center 
of the rectangle dx dy, we have 


(roe) + [(roe + See ay) 








=0. (6) 





— (rr) S [ (r+ Fe #) | F =0, (7) 


or, neglecting infinitesimals of higher order and dividing by dx dy, 








gee erec 
‘The equations of equilibrium then become 

ae Oey _ ooy Otay _ 

e+ Fm io, Bey Fao, (58) 


The two Equations (58) involve three unknowns oz, 0, try. SO 
far we have utilized only the mechanics of rigid bodies. To find 
a third equation we must also consider the elastic deformation of 
the body. A statement of the fact that the elastic body remains 
continuous after deformation leads to the equation 


Breeton) + Glerta) = 0, (9) 
known as the equation of compatibility (sec (28]). 


One of the techniques used in solving these equations (two of 
equilibrium and the equation of compatibility) consists in introduc- 
ing a function ¢(x,y) called the Airy stress function, with the 
Properties that 

an HB mi AE 
If os, 0,, and 72, are given in terms of such a function, the equations 
of equilibrium are automatically satisfied, for 


£(B)- (0%) -« $(B)- (03) -° 
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In terms of 4 the equation of compatibility becomes 





a a 
E(F+t) +8) -° 
or 
Om 
ae + aaa + at =o 
‘The operator 
We = WEDD — a + 2 ara +e 


that is, the operator ¥? operating on V? is the square of the Laplacian 
in two dimensions. The equation satisfied by the Airy stress 
function may therefore be written as 


Vis = 0, (61) 


and is known as the biharmonic equation. Once the equation 
V4 = 0 is solved subject to the given boundary conditions, 
Oz, y) Tay Can be determined from Equations (60). 


14, Additional equations 


Since many partial differential equations cannot be derived 
without a specialized knowledge of particular branches of mathe- 
matical physics, their derivation would be out of place in a book 
of this nature. However, for the sake of completeness some of 
the more important equations appearing in the various branches 
of engineering will be listed here. References where more 
detailed information concerning these various equations may be 
found are given in connection with each equation. 

We have considered only linear partial differential equations, 
that is, equations in which the unknown function and its derivat- 
ives appear only to the first power. However, the reader should 
not infer from this that nonlinear equations are unimportant; 
many of the equations of hydrodynamics and relativity fall into 
this category. 


Google 


38 DERIVATION OF PARTIAL DIFFERENTIAL EQUATIONS Art. 14 
1. Equations of Com; 
(49Ve, + Beto +o) =0, 


in Three Dimensions : 





(1 +») Vey + Breeton te) =0, 


(+9 Ve + Beater te) = 0, 


a 


O49) Vin + oa 





Cet oy +o) = 0, 
a 
(1+ 9) Wires + gay Ce + oy +0) = 0, 


a 
(1+) Virey + ae py (e+ Oy +01) = 0 


where » is Poisson's ratio, c,, 04, , are the normal stress compo- 
nents, and ry, Ter» Tey are the tangential stress components (see [27], 


(28). 


2. Bending of an Elastic Plate : 


ou _ 4g 


Hu Ou 
wet aa t yD (62) 


where u is the deffection, q is the loading per unit area, D is the 
flexural rigidity (see [29]). 


3. Buckling of an Elastic Plate : 


Ou Hu Au 
& tage t oe 
1 ou au au 
= 5 (Na + MH + NV Ze), @) 


where u is the deflection, D is the flexural rigidity, N,, N, are 
the normal forces components, and Ney is the tangential shear 
component per unit of length (see [29]). 
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4. Torsion of a Prismatical Bar : 


Bit B21 ™ 


where ¢ is the stress function a4/2y — r,s, 04/2 — —tyy where 
ta, and ry, are the x- and y-components of the tangential stress 
(see [27]). 


5. Equations of Motion of a Fluid : 


au, me ae au 1a 
Rt Rt Ate TX ae 
2, 8, 1a ae 1% 
at nt esS tea’ sy 
oo ow ow 


ow 1 
Rt Rte tee = 2 ae 


where X, Y, Z are the components of the external force, u, 9, @ are 
the components of velocity, p is the pressure, and 8 is the density 
(cee [16]). 


6. Laminar Flow Equations: 


me 1m 
aR 4 ot 1B my 
om Lm 
Ce 
Ou oo 
B+ 2 a0, 


where u, v are the components of velocity, 8 is the density, p is 
the pressure, and j: is the coefficient of viscosity (see [16]. 


7. Dynamical Equations for an Incompressible Fluid : 
eu Cae 
o@ = ox - 2+ avy, 


Ps # 
be = BY — 5 + wv, 


ao * 
5% 52-2 + ue, 
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where X, Y, Z are the components of the external force, u, v, w are 
the components of velocity, p is the pressure, 5 is the density, and 
pis the coefficient of viscosity (see [16]). 
8. Vibration of a Microphone Diaphragm : 
ou Lae ( ou a) = 





a@ftSE— (65) 


where w is the displacement of a diaphragm in a capacitor micro- 
phone, and a and b are constants (see [18], [19]). 


9. Oscillation in a Loud Speaker Horn : 


eu ou 1 Ou 
oe tha oa =o (@) 


where u is the velocity potential, ¢ is the speed of sound, and f isa 
constant (see (18), [19]. 


10. Maxwell's Equations : 

















+ = 4m, 


ox 


where E,, E,, E, are the components of electric field strength, 
H,, H,, H, are the components of magnetic intensity, D,, D,, D. 
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are the components of electric flux density, B,, B,, B, are the 
components of magnetic flux density, c is the speed of light, 
y is the conductivity, and p is the charge density (see [26]). 


PROBLEMS 


1. Verify the dimensions of @ in the equation 
a 
oF 
for the vibrating string and prove that the equation is dimensionally 
correct. 
2. Where, in the derivation of the equation of the vibrating string 
(Section 2), were the various “idealized” assumptions used ? 
3. Write the boundary conditions for the vibrating string under the 
conditions that the string is initially undisplaced, fixed at both end points, 
tightly stretched, and given an initial velocity 





Vz sin e 
(L is the length of the string.) 


4. A string of length Z has ita right end fixed, while the left end 
point is given a displacement f(t), starting from rest. Write the cor- 
responding boundary conditions. 

5. Assume that the string ia fixed at both end points and given an 
initial displacement as in Figure 19 below. Determine the appropriate 
boundary conditions. 





y 
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6. Check the dimensions of the equation of the vibrating membrane 
and of the vibrating elastic solid. 

7. A membrane has the shape of a square of side a units. It is given 
an initial displacement « = <y(x — a\(y—a) and has no initial velocity. 
Write the boundary conditions for the transverse vibrations of the 
membrane. 

8, Determine the boundary conditions for the longitudinal vibrations 
of an elastic beam if the initial displacement of the beam is f(x), the 
initial velocity is g(x), and the beam is built-in, 

9. Derive the equation for the longitudinal vibrations of a bar under 
the assumption that the cross-sectional area varies linearly with the length 
of the bar. 

10. If the initial displacement of a cantilever beam is a given function 
of =, say 4(z), and the initial velocity yx), determine the boundary 
conditions satisfied by the longitudinal vibrations of the beam. 

11. Consider a shaft of length L which vibrates torsionally. Under 
the assumption that all cross sections of the shaft remain plane during 
rotation, derive the equation satisfied by the angle of twist @. Assume 
the cross section of the shaft to be a circle of area A. 

12, Determine the boundary conditions satisfied by the angle of twist 
of a shaft vibrating in torsion under the assumption that at t = O the 
angle of twist is a known function of x, the angular velocity is a known 
function of x, and the shaft is built-in at both ends. 

13. Write the boundary conditions for the torsional vibrations of a 
shaft as in Problem 12, but assume that the shaft is built-in at the end 
x = 0 and is free at the end x = L. 

14. Check the dimensions of the equation for the transverse vibrations 
of a beam (Section 6). 

15. Let a beam of length L, buil at x = 0 and free at the end 
* =L, besubjected toa force F = t atthe end x = L. Determine 
the boundary conditions satisfied by the transverse displacement of the 
beam. 

16. If the beam of Problem 15 is free of external forces, but its end 
= = 0 is subjected to a motion y = Yosin t, write the boundary 
conditions satisfied by the transverse motion of the beam. 
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17, The force per unit area on the end x = 0 of a semi-infinite 
bar is f(t). If the bar is fixed at x = 00 and the initial displacement 
and velocity of each section are zero, set up the differential equation and 
boundary conditions for the longitudinal displacement u(x,t) of the bar. 

18. Consider a uniform circular disk whose entire surface is insulated. 
Assume that its temperature at ¢ == 0 is a function only of the distance 
from the center of the disk. Starting from basic physical laws, show 
that the temperature u(r, t) of the disk satisfies the equation 

ou ‘Ou 1 aw 
nox ( * +3) 

19. A thin bar is conducting heat along its length and is also radiating 
heat from its surface. The temperature of the surroundings is T. Assume 
that the radiation obeys Newton’s law, that is, the time rate of change 
of temperature due to radiation is proportional to the temperature 
difference between the body and its surroundings. The bar has a 
cross-sectional area A and a density 8. Derive the differential equation 
satisfied by the temperature. 

20. Consider the flow of heat in a long thin bar (one-dimensional 
case). If the bar gains heat due to radioactive decay, determine the 
equation satisfied by the temperature s(x,i). Assume that the rate at 
which heat is gained per unit volume by radioactive decay is proportional 
tore, 

21. Consider the one-dimensional heat equation in the case where 
heat is gained by the decomposition of the material in the bar. Derive 
the differential equation satisfied by the temperature if the rate at which 
heat is gained per unit volume by decomposition is a constant, A. 

22. Derive the equation of the vibrating string in the case in which 
an external normal force of F(z) ibs/ft acts on the string. 

23. Write the boundary conditions for the one-dimensional heat 
equation if the ends of the bar (of length L), instead of being held at 
2 constant temperature, have no temperature gradient across their surfaces 
(adiabatic case), Assume that the initial temperature distribution is f(x). 

24, Derive the telegraph equations for the case of a direct current 
transmission line, that is, one in which the inductance and capacitance 
may be neglected. 

25. An electric transmission line of length / tially in the relaxed 
state, that is, the currents and voltages in the line are zero. At f = 0, 
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a voltage E = By sinwt is applied to the sending end (x = 0). The 
other end of the line is connected to a pure resistance R. Write the 
boundary conditions satisfied by the voltage (x,t). 

26. It was shown in Section 10 that the function I/r, whete 


r= VE—F TO 
‘satisfied the three-dimensional Laplacian, 





Des 


ea, @ 
Vegtot a 
‘Show that the two-dimensional Laplacian 
a @ 
V-wat oe 
is satisfied by the function log I/r, where 
r= Ve—aFP + oe 

27. Derive the equation of continuity of hydrodynamics on the 
assumption that the density p varies linearly with x and is independent 
of time and the y- and s-coordinates. 

28. A horizontal cantilever beam of length L and slender rectangular 
cross section carries a uniformly distributed load on its upper surface 
(see Figure 20). Determine the boundary conditions satisfied by the 
stresses oz, Oy, and Tey. 


GM LEA 
ee 
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Chapter II 
FOURIER SERIES 


15. Preliminary remarks 

In the calculus the reader considered the problem of expanding 
a function f(x) in a power series ([24], Section 5.9) and obtained the 
result that 


Fe) = 10) + fox + FO 


We chose an expansion in powers of x because these are simple 
functions to handle and are amenable to numerical calculations and 
the operations of the calculus. For example, one way to inte- 
grate e** is to expand e" into a power series and integrate term 
by term. 





asf LOs +. 6) 


fx) 





Fig. 21 


In many engineering problems periodic functions are used. If 
f(z) is a periodic function of period 7, that is, f(x + 7) = f(x) 
(see Figure 21), it would appear logical to try to expand j(x) in the 
simplest periodic functions, namely, sines and cosines, rather than 


45 
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in powers of x. We shall therefore attempt to obtain a repre- 
sentation for f(x) in the form 

Slt) = Ag + a 008 x + ay 008 2x + 

+ 6, sin x + by sin 2e + 





= Ay +S ay coons + DS bg sin ns, 


Such a representation is known as a Fourier series, or the Fourier 
series expansion of the function f(z). 

‘The reasons for wishing to obtain such an expansion of a periodic 
function are analogous but more far-reaching than those given 
above for an ordinary power series. Fourier series expansions 
will enable us to solve various partial differential equations and 
will lead to important generalizations like the Fourier transform. 

Our immediate problem is then to obtain an expression for the 
coefficients Ay, a,, and 5, of the Fourier series a8 was done in the 
case of power series. This requires the knowledge of certain 
integral formulas derived in the following section. 


16. Some integral formulas 


In order to proceed in the most efficient fashion it will be 
convenient to recall certain definite integrals, 


[008 ms coe ne dt = 18m 
[sin ms sin nx de = Sag) m,n integers, (68) 
[7 sin mz cos nz dx = 0 
where Sq is the Kronecker delta defined as follows : 
bun = 0 if mem, 
Sma 


The formulas above, Equations (68), are valid for all nonnegative 
integral values of m and m not both zero. They are known as 
orthogonality conditions. For n = m = 0, it is clear that 





1 if m=n 
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7,008 Or 00s Ox dx = f" dx = 2n, 
J, sin Os sin Ox dx = ["0 de = 0, 
J, sin Os cos Ox dx = "0 de = 0. 
Since the proof of all the above formulas rests on the trigonometric 
identities 
co ms cos ne = J-c08 (m+ m}x + 4-008 (m —n}s, 


sin mz sin nx = cos (m — nx — 4 cos (m+ n)s, 


sin mz cos ns = {sin (m+ nlx + 4 sin (m —n)s, 
we shall establish only the first of Equations (68). 


[i conms cos ne de = + con (mt abede + A. co (wn —e de 


1 sin(mtme * ) 1 
DT mtn 1,72 om 

=0 if mt+n40 and m—n#0. 
But m + nand m— nare different from zero if m and n are distinct 
nonnegative integers, Hence 





ae 





fico me cosmeds = 0 if mon. 
If, on the other hand, m = 2, then 
J, 603 mx cos ms ds = |" cost me de 
= pt (Aeon 2me) 
= f. (4S) « =m 
‘These two results establish the validity of the first of Equations (68). 


17. Calculation of the coefficients 


We now are in a position to derive formulas for the coefficients of 
the Fourier series expansion of a function j(x). We shall assume, 
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for the present, that the function f(z) is periodic of period 2m 
(the period if sin x and cos x), and later shall extend our results to 
functions with periods of arbitrary length. 

Let us write the Fourier series expansion of f(z) in the form 


fia) = 3 + Saconne + Fb, sin nx. (69) 
(The reason for using $2, rather than ay will become apparent 
shortly.) Multiply both sides of Equation (69) by cos Oz = 1, 
and integrate from —m to +7: 


[fe as = Sizeode + Sas foo ns con ede 


+ Btu ft in mx con ox as. 


But by the orthogonality conditions (Equations 68) all the terms 
of both summations are zero. Hence 


fifa) a fijeoe = ag 


a= + fife) ae. (70) 


To continue this process we shall multiply both sides of Equa- 
tion (69) by cos mx (m > 0) and‘integrate from —r to 7, viz., 


and 


[[Aayeoame de = La [" come ds + 3 a, [" cos ux cosmzde 


+ 6, |", sin nx cou ms de. 
By the orthogonality conditions, 
Jf, :08 ms de = 0, 


Dn [200 me con me de = Taq, 


at 
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6, [*, sin nx cos mx dx = 0, 
zef 
that is, all terms of the Fourier series vanish except the one involv- 
ing cos ms. Hence 
[fe cos mx dx = aq, 

and on changing m to m, 

1 

a= & [i Ae) coe ne de, wim AB ace 


Note that since we used $a in our Fourier series expansion, 
4, (see Equation 70) is also given by the above formulas, namely, 


1 Lp 
a = = Ji A) cos Ox de = = [pea 
In order to calculate the b,’s we multiply both sides of Equation (69) 


by sin mx and integrate from —n tom. 


ff) sin ms ds = $e fain me de + pr 


+35, sin ne sin me de. 
a 
By our orthogonality conditions, 


fi,sin mx de = 0, 


Say [7 con ne sin ms de = 0; 


at 


yo J, sin nx sin ms ds = thy. 
Hence 


J Ae) sin ms de = thy, 
and on changing from m to 1, 
b, = 2 J" fe) sin nx de. 


Google 


50 FOURIER SERIES Art. 18-19 


18, Recapitutation 
Summing up the results of the previous section we may say: 
If f(z) is a periodic function of x with period 2m, then 
fis) = Lay + 3 (aconns + b, sin m2), 
a 
where 
aes 1h fen com ne ae n= 0,1,2,-.5 
am 





Lf fe oi nx as B= 1,23, 006 


‘The interval (—n,n) is known as the full Fourier interval. Since 
the function f(x) ia periodic of period 2m, once it has been defined 
for —m <2 < =, it is determined for all values of x. 


19. Example of a Fourier series 


In this section we shall consider an example of a Fourier expan- 
sion of a given periodic function. The function to be considered 
is known as a square wave and is defined analytically in the full 
Fourier interval by the equations 

fe) =—-4 —a4<2<0, 

fe) = +h O<x<z, 
where h is a constant. This function is graphically represented in 
Figure 22. By Equations (71), 


aim LJP Ae) cos ne a 


os + 1 fice ne de 








_ hea? (‘22) =) Ld, ake 
since sinQ = sinnw = sin(—nn)=0. Hence all the a,, 
= 1,2, .., are zero, The case » = 0 cannot be evaluated 
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by the above formula since we cannot divide by zero, but a direct 
calculation yields 


a= 2 fteyds =P a+ l fraa 0. 


th) 





Fig. 22 
To calculate the },’s we write 


8, = +f" fe) sin ne ae 
HAC msinne de + 1 [isin ne ax 


2h 


i Ll — con mr, 


Zi — con (— ey] —  fo0e we — 1) = 
since cos (—#) = cos. Also, since 


cosO = +1, 
cos 7 







cos 2m = +1, 


and in general 
cos nr = (—I)", 
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we find that 
2h 
4 =20-cn 
Thus if is odd, 
be 





2h 4h 
Bie 
if nm is even, 


5, 





2h 
mn qo. 


and we conclude therefore that 
@=0 2 =01,2.5 
4 =0 1=2,46,... 


== n= 1,3,5, 





we) - & 





‘The function f(x) can also be written in the form where the dummy 
variable runs from one to infinity by consecutive integers by letting 


n= 2m —1. 


Then ifn = 1, m = 1; if » = 3, m = 2, ete. and 
— 44S _sin Qm = 1p a 
Al) 7 eT (23) 


We sce from this equation that a square wave can be represented 
as the sum of an infinite number of sinusoidal components. It 
-would appear logical, therefore, that if we took only a finite number 
of sine waves we should get an approximate value of f(x)—the 
larger the number of terms the better the approximation. This 
question will be investigated in Section 27. 
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20. Change of interval 
We have assumed f(s) to be periodic of petiod 2m a8 indicated 
in Figure 23. Hence we may write the Fourier series expansion 
of lx) as 
Se) = +4 + J (G cosnx + 6, sin nx), 
a 


to) 





Fig. 23 


an = © ff Ae) coo me de, 


a= 1 fA) sin ne de 


In certain applications it is convenient to use the interval 
(a, « + 2m) where « is any real number, rather than the interval 
(—17). We shall prove in this case that the coefficients a, and b, 
are given by 


where 


i 


1 fotze 

H [* ffx) cos nx de, 

ei ai 
b x i flee) sin nx de. 

Break the above integrals into three integrals, then, for example, 


Jo" 12) 008 nx de = [7 se) 008 ae ae + J Ae) cos ns as 


+ JP Re) cos me dx. (75) 


In the last integral on the right make the change of variable 
942m 
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‘Then 
JE" He) 008 nse de = [fly + 2m) cos my + 2n) dy 


= JX Ay) 008 ny dy, 


since the function f is periodic of period 2. On noting that the 
second integral of Equation (75) is precisely ma, and that 


fr==ff 


LEA) con ne ae es ag 


A similar proof holds for b,. 
Commonly used values of a are O and —n. In the case a = 0 
we obtain. 


we obtain the result that 


‘ages 4 * Aix) cos nx dx, 


bn = 2 ff) sin ms ae, 


while in thecase « = —we obtain the standard formulas of Equa- 
tions (71). 


21. Example of change of interval 

To illustrate the gain in simplicity that can be achieved by the 
use of the formulas of the last section, consider the problem of 
expanding the function f(x) given graphically in Figure 24 and 
analytically by the equations 


fe) = Let Qn, et t<z<h (1H) 
100 
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Of course we could define f(x) in the full Fourier interval (—z,) 
by the two equations 


fe) =A e+ Qn-N) <a <d, 
fe) = L@-0 l<z<a, 


but Equation (76) is a simpler representation since it involves 
only one analytical expression. If we let 


fle) = B+ 3 (e conne + by sin nx), 


we may write by Equations (74) with « = —27 + 1, 


Ef Amma 


I 


lp 
= Sheet (e+ 2 — 1) con ne ds, 


rn 
mw Juaet1 


Als) sin mx de 
= fi De + 28 — 0 sin ns de, 


With the aid of a table of integrals (see [10]) we readily compute a, 
and by : 








_ 1 (couns ssinas) met dat) || 
oe et ed Lae tte) Line 
= 2m 20, 
or 
4% = 2, 
and 
1 (sinne x coane) |! Qe —1 (_ cos nz)! ! 
pyc (Se eons) _ oon ne 
GFL, CAA) 
~ 20m wT 28h noe 
7 
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Hence the Fourier series expansion of the function f(x) of Figure 24 
is given by 
2 S (sinn « 
fe) = 1+ 2S (SP cons = 08. sin nx). (1) 


22. Even and odd functions 
A function f(x) is said to be even if 


2) = fiz), (78) 
and is said to be odd if 
f(-*) = fe). (79) 
For example, cos x is an even function, since 
cos (—x) = cos x, 
and sin x is an odd function, since 
sin (—2) — —sin =, 


In the cases in which the function we wish to expand into a 
Fourier series is even or odd, we can deduce simpler forms for 
its Fourier coefficients a, and b,. Suppose f(z) to be even, then 


L pe 
a, = = [i Ae) 208 me de 
Lo L pm 
= SJ fe) 00 me de + a Fhe) con me ds. 
In the first integral on the right let x= —z. ‘Then 
L jo 1 fo 
FJ Ae) cos me de = + f° f(x) cos n(—z) a2) 
Lopw 
= = [3a cos nz de, 
since interchanging the limits of integration changes the sign of the 
integral, f(—z) = f(a) and cos nz = cos (—nz). If we change 
the dummy variable from z to x, we may write a, a8 


aces L fie) con ne ae ae 2 fh) con ne ax 


— 2 fi ste) con ne de. 
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Similarly, the integral for b, becomes 
8, = Ef Ae) sin nx de 
ae © f° sta) sin ne ae + 2 fia) sin ne ae. 
As before, let x = —a in the first integral on the right. Then 
Lf? fee) sin me de = 1 f° ps) sin mfx) dls) 


= — fiAe) sin na de 
since j(—s) = f(z) and sin ns — — sin x(—2). We conclude 
therefore that 
b= Af" fle) sins de + 2 fifa) sin me ae =0. 


Hence, if f(x) is an even function, b, = 0,1 = 1, 2, 3, ... and we 
may write its Fourier series expansion as 


Ae) = te + 6, 00 me 
at fiz) even. (0) 
eis 2 fp) cos ne dx 


Note that half the Fourier coefficients (all the },) are zero and that 
the integration of f(x) cos nx is carried out over the half Fourier 
interval (0,7). 

A similar result holds for odd functions. For, if (x) is odd, 


Bien L ft Ae) sin nx de 


LP fey sin me de + Efi Ae) sin ne de 


fre sin n(—2) d(—2) + 2 fips) sin ms ae 


2 fife) sin ns ds, 


ql 


Google 


58 FOURIER SERIES Art. 23 


ied = Lf Ae) coone de 


= 1 se) coons de + + ff) con ne de 


= LAAs) coo m2) ta) + © ft fe) cos ne de 
= 
‘The Fourier series expansion for an odd function is therefore 
given by 
fe) = by sin ne 
a fle) odd. (81) 
theme 2 Ji Ae) sin ne ds 
In the example of Section 19, the function was odd. Hence all the 
cosine terms were zero, a fact we verified by actually computing 
the a, and 6,. Note that for an odd function no constant term a, 
is present. 


23. Examples of expansion of even and odd functions 
Let fla) be a function defined in the full Fourier interval as 
fs) =ot et neces (82) 
(see Figure 25). Clearly f(x) is an even function, since 
fia) = ot — (Cat = tat = fle). 


ta 





(Art. 23 POURIER SERIES 
Hence, by Equations (80), 


a, = 2 f(x — 35) conned 





= 2 Ags (sins) PE _ [2e coon 

72 Cr) --aes 
-- n#0, 

and 

ao = 2 fit — ay ae = Sot, 


(tx! — 2) sin =] =| 
” o) 


‘The Fourier series expansion of f(x) is therefore given by 


fG) 


con ne. (83) 





Fig. 26 


As an example of an odd function, consider the function f(x) 
sketched in Figure 26 and analytically defined by 


fie) = 20 a << 
fle) = -1 -$<< 
fe) = +1 O<e< 
Ae) = +2 jes< 


2 
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By Equations (81) the b, of f(z) are given by 


2 pr. 2¢ , 
by = = PM sin ne ds + = J" ,2 sin nx de 


oS ak 





or 








=f if 81,35 
am 
b, = 0 ifn = 4,8,12,..., 
ee” a a = 26 Aice, 
=m 
since 
cos 
if n is odd, and 
oon F = (-1" 


if mis even (see Section 28). ‘The Fourier series expansion of f(x) 
is therefore given by 


fey = & -4 3% ww 





24. Half Fourier intervals 


In many applications we are interested in a function f(x) defined 
only in the half Fourier interval (0.7). To expand such a 
function all that need be done is to prolong it into the other half 
interval (—7,0) and then expand this function—now defined in 
the full Fourier interval (—n,n)—by the usual Fourier series 
techniques. This is possible because no matter how we extend 
the function f(x) into the interval (—7,0), it will still represent 
the desired function in (On), as well as in (20,3n), (dm,5m) etc. 

For example, consider the function f(x) = x defined in the half 
Fourier interval (Oj). We may prolong this function to the 


Google 


Art. 24 FOURIER SERIES 61 


half interval (—7,0) in an infinity of ways. In Figure 27a we 
have extended the function as an odd function, in Figure 27b we 
have extended it as an even function, and finally in Figure 27c we 
have prolonged f(:) to (—n,0) so that it is neither an odd nor an 
even function. ‘The series expansions of f(z) as given by Fig- 
ures 27a, 27b, or 27c all represent the function f(x) = x in the 





Fig. 27 


interval (0,n). If we have a choice in the prolongation of a 
function, it is advisable to prolong it in such a way as to reduce as 
far as possible the effort of calculating the coefficients of the Fourier 
series. This would indicate that an even or odd prolongation is 
preferable, since it eliminates the necessity of calculating half the 
Fourier coefficients. Sometimes the cosine rather than the sine, 
expansion, or vice versa, is required by the nature of the problem, 
as will be seen in the examples of the next chapter. Also, even if 
we have 2 free choice of even or odd prolongations, one may be 
more desirable than the other, as will be illustrated in the next 
section. 
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25. Examples of even and odd prolongations 

Let us consider the problem of expanding the function y = x 
given in the half Fourier interval (0,7) into a cosine series, that is, 
of prolonging y as an even function (see Figure 27). In this case 
we know by Equations (80) that 


_2p _ 2 [esinns , cosne 
a= 2 fps smnede = 2 [ 2 te 








lo 


2 
=A lay- 2 #0, 
or 
4 =0 f= 24,605 
a 
and, for » = 0, 
a = 2 [lads =o. 
Hence 
7_4 ne 
s-t-4 3% Se Osecn (85) 





Let us now expand f(x) = », given in the half Fourier interval 
(On) into a sine series (see Figure 272), that is, let us prolong /(=) 
as an odd function. For this case we know by Equations (81), 


ben 2 fprsnmeae = 22 oe) 


= 2 (-28") =-Fo 


7 








lo 


and hence 


x= aZer ans, —n<ercm. (86) 





On comparing this expansion with Equation (85), we see that x as 
given by the expansion of Equation (85) converges more rapidly 
than the expansion of Equation (86), since the terms in Equa- 
tion (85) are decreasing as the square of n, and moreover half the 
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cosine terms are missing. Hence, if we are free to prolong x in 
the interval (—7,0) as an even or odd function, the even pro- 
Jongation is to be preferred in this case. 


26. Points of discontinuity 
‘The question arises as to what value the series will yield at 
points where a function is not continuous. It may be proved 
(see Section 40) that the series yields the average of the values at 
the left and right of the points of discontinuity. Precisely, if 
is a point of discontinuity, 
— AE) + G+) 
Ae) = z » 





where a 
KE-) = Jim fia), 
ret 
and 
AG+) = lim fiz). 


se 


In the example of Section 19 (the square wave), the function is 
discontinuous at the points x = 0,7, +2r, ..., and 


lent) = +h for neven, 
fmt) = —h for n odd, 
fmm) = —h for neven, 
fen) = +h for odd. 


Hence the sum of the series at «= nm equals 
Stor) = 2h F* oo, 
Similarly, in the example of Section 21 (Figure 24), 
Al +2en) = 1, 2 =0, +1, +2. 
and in the example of Section 23 (Figure 26), 
2) _3 
f=0 (5) =F 
7 2 
4-3) =-3 fm=0 
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21. Trigonometric approximations 
Let f(z) be a periodic function with period 2m. The question 
sometime arises as to what is the ‘best approximation ” to f(x) 
obtainable by using only a finite number of sines and cosines. By 
“best approximation” to f(x) we usually mean a function A(x) 
which makes the integral 
T= fl [fe) — A@pae (7) 
as small as possible. ‘The integral J is called “the mean square 
error” and A(s) approximates f(x) “‘in the mean square sense.” 
It will now be shown that the best approximation A(2) to f(x) 
is one in which the coefficients of the trigonometric terms are pre- 


cisely the Fourier coefficients of f(z). Toward this end consider 
the finite series 


Aw) = B+ se cos nx + b, sin nx) (88) 
in which the a, and 5, are to be determined such that 
[Ut — aeypex 
shall be a minimum. Denote the above integral by 


Hay Ay,---Ansbyby--sbn)s 


that is, J is 2 function of the undetermined coefficients a, and by 
only. By the calculus the derivatives of J with respect to the a, 
and b, must equal zero for J to be a minimum. Hence 


o=- 2-28 ~ Ate) 54 de 


2 fue — A(z)] cos kx dx, 


or 


JT Ax) cos ke de = Si. (feowke) de 


N, 
+ [2 em cov me + by sin nx) cos kx dx. 
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By the orthogonality conditions of Equations (68) the term on the 
right is za,; (k = 0, 1, 2, ..., N) and hence 
a= 1 ft sta) con bs dx, k= 0,1, 4. 


In a similar fashion, on differentiating J with respect to b, and 
setting 41/2, equal to zero we obtain 


b= Xf sta) sin ts as, k= 1,2, 





vo Ne 





It is thus seen that the a, and b, of Equation (88) are precisely the 
Fourier coefficients. For example, the best approximation to the 
square wave (with h = 7/4) of Section 19, using one term, is 

Ae) = 
Using two terms, the best approximation is 





ea 


Sz) = sins + Fein 3x, 


and using N terms, the best approximation is 
at 


nls) = Sy 








The graphs of the functions f(x), j,(2), f(x) are plotted in 
Figure 28. 
28. Trigonometric identities 

In the Fourier series expansions of the previous sections we made 
use of certain trigonometric identities like cos mm = (—1)" for 
n=0, 1,2, ... and cos (7/2) =(—1)"" if nis even. Itis con- 
venient to have for reference the following table of similar formulas. 







cos nn | (-1" | +1 +1 +1 
na 

sin o o ° 

con 3 | (ay: =i 41 

sin Vystar cpm | 0 
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29. Fourier series expansion in (—L, L) 


If flz) is @ periodic function of period 2m we have seen that it 
has the Fourier series expansion 


fe) = te + ¥ coon + 6, sin nz), 


where 


a= J, Al2) one as ax, 


b= Ef A) sin mx a. 


In many engineering problems, the interval of interest ie not 
restricted to 2 but may be of arbitrary length, say 2L. To 
obtain a Fourier series expansion for a function f(x), which is 
periodic of period 2L and is given in the fundamental interval 
(EL, L) we introduce the change of variable 


e=fe (89) 


fe) - (2) = 40 


7 


and let 


Here d(z) is a new function of the variable x which has period 2, 
for we see by Equation (89) that when x = —L, 2 = —x, and 
when + = L, «=. Since 4(2) is periodic of period 2n we can 
expand it as was done in previous sections, obtaining 


He) = te + E (acorns + by sin ne), 
where 
om = 5 ff $2) coo ns de, 


8, = © fT ge) sin ne ds. 


If we now return to the original variable x, we have 


He) = I0) = Jon + Zee Be th sin x) (90 
x 
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where 
an = 5 fife Te d( Ts) = 7 Ji se) cos Te de 

oy 


Aft, fe) sin «($) £ fife) sin Se de. 


We have thus obtsined the expansion for a function f(z) of 
arbitrary period 2L, and as mentioned in Section 20, we may also 
write 


1 pean 
a = TIT Me) 008 Fe de, 
1 pete 
wad tena Poa 
where « is any real number. A commonly used value of « is 
zero. In this case, the Fourier coefficients become 
1 pu 
em rh Se) cos Fx de, 
(92) 
1 Fates 
bn = Jp fe) i 





30. Example of a Fourier series in (—L, L) 
As an example illustrating the formulas of the previous section, 
consider the function f(s) defined in the interval (—I,1) as 
fe)= 2 -l<xr<k 


If we repeat this function periodically assuming (—1,1) to 
be the full Fourier interval, we obtain the function sketched in 
Figure 29. Since 2L = 1 — (—I) = 2, L = 1, we have by 
Equations (91), 


a, = 





x cos mmx dx, by = 
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We readily compute the a, and by: 

a, = a [Batetst — 6) cos nx + (n*n’s? — 6nnx) sin onal" 

1 


= 0. Wee Dy Zany 

a = |, ade =o 

bg = ee (Betts — 6) sin rnz — (nina — Genz) conn)" 
= Wa) py, 





Hence, by Equation (90), 


fe) = 3 x Goa) (— 1)" sin anx. 


31. Even and odd functions in (—L, L) 


If f(z) is an even function of period 2L, we can write for it an 
expansion analogous to that of Equation (80) by the techniques 
used in Section 22. We leave it to the reader to verify that in this 
case 


fle) = Lay + Zace Ts 
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2 
4% = Efe e)em Fady n= 0,1,2... (3) 


Correspondingly if f(x) is an odd function of period 2L, the 
expansion analogous to that of Equations (81) becomes 
fle) = Da sin Fe 
(94) 
2 . 
= Z Sofa) sin Fax de, n= 1,2,3,... 


32. Examples of even and odd functions in (—L, L) 
Consider the odd periodic function f(x) 
fe) = + 4<x<4 


illustrated in Figure 30. Since the full Fourier interval 2L is 8, 
its Fourier coefficients are given by Equation (94) as 
La. 


ume 3 dx 


TT | 
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and 
4g in Om 
fe) = — 2 gy ain Ps (95) 


As an example of an even function consider the function illus- 
trated in Figure 31, whose period 2L = 8. By Equations (93) 
its Fourier coefficients are given by 

2 cud 2/4 mm 
a = G yc Pete — | [joo Pade 


2 (.. =m | ar 
Z (oF tan 7) 


Zsin U+Cy), 2 0. 


i 


i 


tha) 





a, = if n is even, 


il 


0 
a £2 ywteminn ifn ie odd. 


(See Section 28). Also 





‘The Fourier series expansion of the function of Figure 31 is thus 


ts 22 A eresim cose, (86) 
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33. Prolongation of functions in (—L, L) 

If f(x) is given in the half 
Fourier interval (0,L), its expan- 
sion may be obtained by pro- 
longing it into (—L,0) as an 
even function, an odd function, 
or neither (see Section 24). 
Consider, for example, the func- 
tion f(z), 

fe)=x OSxeS1, 

Ax) 1sxs2 
shown in Figure 32, which is defined in the half Fourier interval 
(0,2). If we extend this function as an odd function, it will appear 
as in Figure 33. Its Fourier coefficients are given by 


too 





b, = { sin Sx ds +ff sin "Fx de 





= (— Zico Be + isin Se anh ateae 
= cone + aaasin-Te) | 
6, 2 is 
ee 2 if n is even, 
4 ‘igs Sete 
Oy = eg (et — 2a it nis odd. 
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Therefore 
fs) = — 





- dure 


Hyun init gk 
Ae(-e sin 





- 2 > + (-1)" sin (7) 





[We note, of course, that 
2 a: 
Cie 


if n is even, but we prefer to keep the (—I)" factor so that the 
first and third series may be readily combined.] 


lx) 





Fig. 34 


The function defined by Figure 32 could also be extended into 
(2,0) as an even function. In this case, it would appear as 
illustrated in Figure 34 above. The reader may obtain the 
corresponding cosine series expansion for this function (sce 
Problem 19). 


6 
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34. Phase angles 

An alternative form of Fourier series may be considered when 
both sine and cosine terms are present. Suppose we have a 
function f(z) of period 2L given in the interval (—L,L), whose 
Fourier series expansion is given by 


fe) — tag t & (ao Fe + tysin Fe), (8). 


where : 
oi nn 
on = 7 [i Ae) coe Tia ae, 
1 
t= Efi soon Boe 
The sum of two corresponding terms of Equation (98), 
a, con Fx + b, sin Fx, (99) 
is also a sinusoidal function and may therefore be written 
acon a + bysin Ee = egsin (+ 4,). (100) 


To determine c, and qq, expand the right-hand side of the above 
‘equation : 





eatin (Fx toy) — cusin x con fy + cus gy cor Fs 
and equate coefficients of 


er i nr 
sin“Fx and cos “TF 





E 
on both sides : 
Oy = fu 8iN by, by = Cy C08 dye (101) 
Hence on squaring and adding, 
q = (102) 
while taking the ratio, 
dn = are tan ()- (103) 
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‘The Fourier expansion of Equation (98) thus becomes 
1 (0 
fe) = 3% + Zevsin (Fe +4). (104) 
We could also write 
nm am ca 
4, 008 Te + by sine = com (x44); 

and we determine as above that 


Gy, = 6,00 by by = —ty 8 by (105) 
and hence 


= VOFFTR = wewn(—2). (106 


Our Fourier series expansion therefore becomes 
fe) = Fa + Yevcos (Fx +h): (107) 


Care must be exercised in computing ¢, from Equation (103) 
[or (106)] since there are two possible determinations of ¢, which 
satisfy the equation. The value which must be selected is the 
one which satisfies both of Equations (101) for (105). 


35. Example of a Fourier series involving a phase angle 


Consider the function of Figure 24 (Section 21). Its Fourier 
series expansion was determined to be [Equation (77)] 





fe) = 14 2 (+ St coene aia sina). (108) 


To write this series in terms of a phase angle [Equation (104)] we 
have first from Equation (102) that 


«= TEST 
and from Equation (103) 
dn = arc tan (- 





ann) = are tan (—tan 8) = 9—n. 
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The phase angle form of Equation (108) is thus [from Equa- 
tion (104)] 


Ms 


fey 142 


Hein (ne + =n) 


aly 
Ms 


aL sin n(x — 1). (109) 


‘ 


36. Complex Fourier series 

In many engineering applications of Fourier series it is more 
convenient to expand a function into series of exponentials rather 
than into series involving sines and/or cosines. Since the trig- 
onometric and exponential functions are intimately connected by 
Euler's formulas ((24], Section 5.12), 





we can use these relations in Equation (90) to obtain the complex 
Fourier series expansion of the function f(x). If we make this 
substitution and for simplicity let 


= *, a=0, +1, +2, ..., (110) 


the Fourier series expansion of Equation (90) becomes 


Ha) = hay + > a(S 4 =") + = oS 





= font FR te steer + FS lo ttre (uy 


on = Gy — jon (a2) 
and with the Fourier coefficients a, and 6, of Equation (91), we 
have 
- le ue 
6 = iby = Tf, Ae) (008 ene —j sin ent) ds 
1 


=f fife) elntde, (113) 
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Similarly, a, + jb, = é,, the complex conjugate of ¢., 

tn = On tite = Tf Saltoeds. 
But from Equation (110), 

a = GF 2 wy, 


and if we replace n by —n in Equation (113), 


t Ee fisertatds = t f falelmnds = iy (114) 





1 1< 1g. 
He) = a + ES eee + LS ene 


= =(% + Seer 4 dew). (is) 
a a 
Let n = —m in the second summation on the right; then 


Sepia S gms 
x * 


“Sue 





‘Thus Equation (115) may be written 


(a + Dawe + BS cet 
a 





1p 
a= 1 ft teas = oy = caine, 
since 
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so that finally 


fe) = 3 Ee ewie, 


where (116) 
= Li foetveds, 1 = 0, £1, 42, 





= 


Equations (116) are known as the complex Fourier series expansion 
(or exponential Fourier series) of f(x). The set of coefficients {c,} 
is often called the Fourier spectrum or simply the spectrum of f(z). 

‘The complex form of the Fourier series is much more compact 
than the trigonometric form. It is useful, for example, in deriving 
the Fourier integral (see Chapter IV). One disadvantage of the 
exponential form is that it gives no simplification in the case of 
even and odd functions. 





37. Example of a complex Fourier series 


Let us expand the function f(x) = e* given in the full Fourier 
interval (—1, 1) (Figure 35) into a complex Fourier series. With 
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Equations (116) give 
a= fi erchee dy = | 
- TF 1 


di-na) _ ¢-(-a) 





1 jun 
But 
em = cosmm + jsinnm = (—1)" 
and 
‘ ) (abn (t+ jam) (¢ — 4) 
4 T+ nt 7 
Since , 
ene 
z= 1.1752, 





the complex Fourier series of e* becomes, by Equation (116), 


fe) = 1.17982 & (RCE cs, (i) 


38. Double trigonometric series 


In many applications it is found necessary to expand a function 
of two variables into a trigonometric series. Consider a function 
‘flz.y) periodic of period 2n in both x and y. Then, if we hold y 
fixed, we can expand f(x,y) into an ordinary Fourier series, in 
which, however, the coefficients will be functions of y, viz., 


fay) = ¥ [Ady con ne + Byly) si ns). 
% 


If we assume that the functions 4,(y) and B,(y) can be expanded in 
a Fourier series, 


Ax() = (nm 608 my + Pym sin my), 
= 


By) =F (Yam £08 my + tom sin my), 
% 
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and substitute these coefficients in the expression for f(x,y) we 
obtain terms of the form 

cos nx cos my, cos nx sin my, 

sin nx cos my, sin nx sin my, 


with constant coefficients. Hence, in general we may write 


HED) = YY Cam 608 me C08 MY + bpm COS x sin my 
oie 
+ Cam Sin 2x COs my + dam sin nx sin my). (118) 


In order to calculate the coefficients Gan, Bam Cams ANd dam WE USE 
the integration procedure employed in ordinary Fourier series, but 
we apply it twice. Multiplying each side of the above equation 
by sin ax, where a is a nonzero integer, and integrating with 
respect to x from —zn to 7, we obtain, because of the orthogonality 
conditions, 


J Seep) sin ax de = Stn cos my + Dan sinmy. (119) 
ao a 





Multiplying this equation by sin By and integrating with respect 
toy from —z to x, we obtain again, by the orthogonality conditions, 


Te [fr ‘fley) sin ox | sin By dy = md, 


or 
AS, [ten sin ax sin By dx dy, (120) 


In a similar fashion, 





(ag EJ", Jf ea) os ax cos By de dy ) 


1 


4 de [fe cos ax sin By dx dy | (121) 
cs = 4 JL, J fear) sin a coe Ay de dy. | 


Equations (120) and (121) are valid for «, B = 1, 2,3... The 
remaining constants dym, so, Doms Pnor Com» Cnor dom dug ate CO- 
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efficients of terms involving only one variable. Wesee immediately 
that of these constants 

a (122) 
since all these terms involve the sine of zero. To evaluate Cy, 
multiply Equation (119) by cosQy = 1 and integrate from —r 
to m with respect to y: 


fe 3 flsy) sin ox de dy = 2r%q, 





from which 

legume a J, Jo, Aen) sin ox de dy. (123) 
Ina similar fashion, 

buy = sta J, | se) sin By de dy, (124) 


To calculate ayo, n 3 0, multiply Equation (118) by cos ax and 
integrate from —z to n with respect to x: 


[7 fox) cos ax dx = Baan 608 my + 1 Ban sin my, 
oe mnt inno 


and then multiply by cos Qy = 1 and integrate from —n to 7 
with respect to y. 


fe pes flay) cos ax de dy = 2nago, 








from which 

Oa = 53 f. J fsy) 008 cox de dy. (125) 
Ina similar fashion we verify that 

day = she |", [_fe9) cos By de dy. 126) 


Finally, dog is calculated by integrating both sides of Equation (118) 
with respect to both x and y: 


je flsy) dx = 20 z qq COS my + 2m Zim sin my, 


J, J fee) ds dy = Arran, 
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from which 
oe = pe [Ji flood de by. (azn 
It is not our intention to study multiple Fourier series in as 
an exhaustive a fashion as we did Fourier series in one variable. 
Hence for the present we shall consider only a special case, namely 
when f(x,y) is an odd function of both x and y. In this case 
the coefficients 
For Bons Baty Gum Doms bums Coo» Cum 


vanish, and the double Fourier series reduces to 


Sey) = > s ym Sin ma sin my 
a = 
dam = a J”, J” fos) sin me sin my dx dy 


4 pan . . 
= SM fine) sin nx sin my dx dy. 


39. Example of a double Fourier series 

Let f(xy) = k (a constant) in the rectangle 0< x <7, 
0 <y < = and be an odd function of both x and y. ‘The function 
flzy) is sketched in Figure 36 and is defined analytically by the 
equations 





fay) =k O<x<7, O<y<z, 
fey) =k —a<r<0, —9<y <0, 
fey) = —k O<x<% —=4<y <0, 
fly) = —k —4<x<, O<y<n. 


‘Then by Equations (128), 
day = Mf" ff sin nx sin my ds dy 


lk, 
=san «(if n and m areodd, 
=0 otherwise, 
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and the double Fourier series of f(x,y) — h is 





Fig. 36 


40. Mathematical conditions on Fourier series expansions 
Not all periodic functions can be expanded into a Fourier series, 
since certain mathematical conditions must be satisfied to guarantee 
that the sum of the series equals the function. The classical 
sufficient conditions for a function to be Fourier expandable are 
known as the Dirichlet conditions and are listed here in a form 
general enough for most engineering purposes. 
If f(z) is a single-valued function of x, periodic of period 2x 
such that in the interval (—z,7), 
(f(s) never becomes infinite, 
(ii) f(z) has at most a finite number of maxima and minima, 
(iii) f(x) has at most a finite number of discontinuities, then 
f(s) is Fourier expandable. 
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Certain conditions, slightly weaker than the Dirichlet conditions 
lead to the following precise theorem whose proof may be found 
in (6). 

Theorem. Let f(z) be a periodic function of period 27. Sup- 
pose that the interval (—7,n) can be broken up into a finite 
number of subintervals in each of which f(x) is continuous. Then 
the Fourier series for f(s), 


Jay + & os soe ne + by sin ne) 
where 
a, = Lf" fle) cosme de, 2 = 0,1,2,..4 


he tf tnt, as Th, 
converges to the value 
Fhheh) + fe] 


at every point at which f(x) has a right- and left-hand derivative. 

If x is a point of continuity, then f(x+) = f(z—) = fix), and 
the series converges to f(z). 

In summary, we may say that while not all functions possess 
Fourier series expansions, almost any function that is likely to be 
met in a engineering problem is Fourier expandable, and its Fourier 
series will yield the value of the function at every point of con- 
tinuity and yield the average value at points of discontinuity 
(see Section 26). 


PROBLEMS 
1. Establish the second and third orthogonality conditions of 
Equations (68). 
2. Expand the function f(x), 
fe)=4 —-4<x<0, 
fe) =1 O<s<7, 
given in the full Fourier interval (—x, 7) into a Fourier series. 


Google 


FOURIER SERIES 85 


3. Obtain the Fourier series expansion for the function defined in 
(—mm) by Figure 37. 


£0) 





Fig. 37 

4. Write the Fourier series expansion for the function f(x), 

fe) =, —-4¥ ccm 

5. Expand the function fx), 

fe) -—0 weed, 
fe) = O<e<n, 
in a Fourier series. 

6. Find the value of the Fourier series of Problems 2, 3, 4, and 5 
at the points x = —x, 0, 1, 2x. 

7. Determine whether the following functions (considered 28 real 
functions of the real variable x) are odd, even, or neither. 

(a) sec x. 

(b) #. 

(©) sinh =. 

(d) tan x3. 

(©) Ae), f an arbitrary function. 

(0 (sin x), $ an arbitrary function. 

8. Prolong the function f(x) given in Figure 38 into the half Fourier 
interval so that it becomes an odd function, and obtain its Fourier series 
expansion. 

9. Extend the function of Figure 38 into the interval (—z,0) such 


that its Fourier series expansion contains only cosine terms. Obtain 
this expansion. 
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10. Extend the function 
fle) = 2sin3e O<x<n, 
into the half Fourier interval (—7,0) such that it becomes an odd 
function, and obtain the corresponding Fourier series expansion. 
11, Expand the function f(x), 
fe) = cote 4 ccm, 
into a Fourier series. 
12, We are interested in the values of the function 
fa) = 242 —4 
in the interval (—n,n). Obtain a Fourier series expansion that will 
represent f(x) in this interval. 
13. Let f(z) be a periodic function of period 2x defined in the 
full Fourier interval by the equation 
fe) =x -m<s<n 
Calculate the best trigonometric approximation to f(s) using 
(a) One sine term, 
(b) Two sine terms. 
(c) Three sine terms. 
(d) Four sine terms. 
Plot the results of (a), (b), (c), and (d) in the interval (0,7). 
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14. A function ffx) is defined in the full Fourier interval (0, 2c) 
by the equations 
yots O<#<« 
yah -te-9 e<ce<d 
Obtain its Fourier series expansion. 
15, Expand the function 2) given in Figure 39 below into a Fourier 
series, 


0) 





Fig. 9 

16. A voltage ot) = Ey sin cst is passed through a half-wave rectifier 
which clips the negative portions of the voltage. Obtain a Fourier series 
expansion for the output voltage of the rectifier. 

17. If the voltage et) = Eysinwt is passed through 2 full-wave 
rectifier, obtain its Fourier series expansion. 

18. Establish the formulas of Section 31 of the text. 

19. Expand the function given in Figure 34, Section 33 of the text 
as an even function. 


20. Extend the function f(z) given in Figure 40 #0 that it becomes 
an odd function in the full Fourier interval (—2,2). 
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21. Prolong the function of Figure 40 so that its Fourier series 
expansion contains only cosine terms. 


10) 





Fig. 40 
22. Evaluate the Fourier series of Problems 20 and 21 at the points 
= =0, 1, 2,3. 


23. Write the Fourier series corresponding to the function 


Ae) 





O<x<c, 


fe) = —-*e—o cece, 


using sines only. 


24. Write the Fourier series corresponding to the function of 
Problem 23, using cosines only. 


25. Obtain the Fourier series expansion for the periodic function 
shown in Figure 41. 


#() 
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26. Expand the function 
oO <r <d, 
= 3 O<2<2 








into a sine Fourier series involving a phase angle. 
27. Expand the function 
fle) =x  3<* <0, 





fiz) = 0 Oe <3, 
into a cosine Fourier series involving a phase angle. 

28. Expand the function f(x) = e* given in the interval 0 <* <2 
into a complex Fourier series. 

29, Extend the function f(x) of Problem 28 such that it is periodic 
‘of period 4 and an even function, and obtain its complex Fourier series 
expansion. 

30. Write the exponential Fourier series for the function f(z) of 


Problem 28 after prolonging it in such a way that it is periodic of period 4 
and an odd funetion, 


31. Expand the function given in Figure 42 below in a complex 
Fourier series. 





Fig. 42 
32. Write the complex Fourier series expansion for the function 
fe)=-x  -1<x50, 
fie) = = Oss<l. 
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33. Write the double Fourier series for a function f(x,y) of two 
variables in such a manner that the same formula can be used for dg, 
Gems Guts 204 Gym; the same formula for bem and byq; the same formula 
FOF qq 20d Cm: 

34. Write the double Fourier series for the function fiz,y) defined 
analytically in the full Fourier interval by the equation 

Ney) = <a <my, 
—a<ycn 
35. Let filmy) be defined as 
fey) = yO <a<m 
O<y<7, 
in the quarter Fourier rectangle. Extend it as an odd function to the 
full Fourier rectangle and obtain its Fourier series expansion. 


36. Obtain the Fourier series expansion (analogous to Equations 
120-127 of the text) if : 

(2) flx.y) is an even function in both x and y. 

() fix.y) is an even function of x and an odd function of y. 

(©) flz,y) is an odd function of x and an even function of y. 

(d) fey) is an odd function of x. 

(© fiz.y) is an even function of z. 

(8) flxy) is an odd function of y. 

(g) fits) i an even function of y. 


37. Extend the function of Problem 35 so that it becomes an even 
function in both x and y, and obtain its double Fourier series expansion. 


38. Extend the function of Problem 35 so that it becomes an even 
function of x and an odd function of y, and obtain its Fourier series 
expansion. 

39. Extend the function of Problem 35 so that it becomes an odd 
function of x and an even function of y, and obtain its Fourier series 
expansion. 


40. Expand the function in Figure 43 into a double Fourier series. 
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Fig. 43 
41. Given the function f(x,y) defined in the quarter Fourier rectangle 
as 
flay) = ty O<e <2, 
O<y<zn 
Extend this function to the full Fourier rectangle (—1 <x <7, 
—1 <y<-) such that it becomes an odd function of = and y, and 
‘obtain its Fourier series expansion. 

42. Extend the function of Problem 41 so that it becomes a function 

(a) Even in x and odd in y, 

(b) Even in y and odd in x, 

(c) Even in * and even in y, 
and obtain its Fourier series expansion. 

43. Generalize the double Fourier series expansion from the square 
(—4 <x<n, —1 <y <7) to an arbitrary rectangle (1 <= <a, 
—b<y <b). 

44. Obtain the analog of the complex Fourier series for functions 
of two variables. 

45. Deduce the form of the Fourier series and the equations for 
the coefficients of a function f(x,y, 2) which is periodic of period 2r 
in x, y, and 2, and an odd function in all three variables. 
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46. Determine whether the following functions can be expanded into 
a Fourier series or not and give the reasons why. 


@ sat defined in —4¥ <x <x. 


(&) sind defined in 1 <r <2 


On defined in 2 << <2. 


(d) log + defined in 1 <x <4. 
(©) Alz)= 1 — 2-9 for 1-2-8 << 12-9, 2 = 0, 1,2, 
where f(z) is defined in 0 < = <1. 
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Chapter III 
SEPARATION OF VARIABLES 


41. Introduction 

In Chapter I we derived various partial differential equations 
of interest in engineering and physics. In Chapter II the element- 
ary theory of Fourier series was developed. In the present chapter 
we shall combine these results; that is, we shall utilize Fourier 
series to solve the differential equations derived in Chapter I. 
This will be accomplished by the powerful but elementary method 
of “ separation of variables.” 

It is not immediately evident, but it will be shown later (Chap- 
ter VI) that linear partial differential equations of the second 
order fall into essentially three distinct types, of which the following 
are typical examples : 





‘The first might represent the equation of the vibrating string, in 
which case a*= Tg/w. The second might represent the 
temperature distribution in a plate, in which the temperature is 
independent of time. The third might represent the voltage in a 
transmission line of negligible inductance and conductance, in 
which case 1/K = RC is the time constant. We shall show 
that all three of the above equations with suitable boundary 
conditions, can be solved by the method of separation of variables. 

Since most of the equations derived in Chapter I were of the 
second order in two independent variables we have stressed this 
type in the present chapter. However, separation of variables 
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also applies to equations whose order is higher than the second 
(Gection 51) and to equations in more than two independent 
variables (Section 52). 


42. Solution of the vibrating string problem 
‘As a first example we shall consider the problem of the vibrating 
string (sce Figure 44). This figure represents a tightly stretched 
y 


string fixed at the end points 0, L, and subjected to a uniform 
tension of T pounds. We have seen earlier that the motion of 
the string ia governed by the equation 

au au 

R-eS (125) 
where at = Tg/w, with w the weight per unit length of the string, 
and g the gravitational constant, Typical boundary conditions 
were also determined to be 


1 40) = 0, 
2 Lf) = 0, 
3. Fue = 0, (130) 


4. xx,0) = fz) #0, 


where we have assumed the initial position of the string to be 
A). [If flz) = 0 the only solution to the boundary value 
problem is u(z,t) =0.] 

‘The main assumption in separation of variables -consiste in 
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stating that the solution of Equation (129) can be written in the 
form of a function of « times a function of t, 


Mat) = X(x) TI). (131) 
‘This assumption is justified by the fact that it yields solutions to a 
variety of engineering problems. 

If we substitute Equation (131) into Equation (129), we obtain 
we 
x” aT 
where the primes on the functions X and 7 represent different- 
iation with respect to the only variable present. Now the right- 
hand side of the above equation is independent of x and the 
left-hand side is independent of ft. Since they are equal, their 
common value cannot be a function of ¢ or x, and must therefore 
be a constant, say A: 

x” rr 

y= A i a 
KY" — dX =0, MT —dMT=0. 
‘The partial differential equation of Equation (129) has thus been 
reduced to two ordinary differential equations. 

It will be convenient to determine whether A is a positive 
negative, or zero constant before we solve the above equations. 
(See Problem 20 at the end of this chapter.) If A>0, the 
solutions of X” — AX =0 will be 


evs and etVis, 

that is, exponentials, This assumption contradicts the physical 
fact that the motion of the string is oscillatory in the y-direction. 
‘The assumption that A = 0 leads to the solutions, 

X@) = Ax+B, Tt) = Ct+D, 
from which by Ecuation (131), 

ux,t) = (Ax + B) (Ct + D), 

where A, B, C, D are arbitrary constants. Such a solution is 


clearly not oscillatory and furthermore could not be made to 
satisfy the boundary conditions. Our remaining alternative is 


XT’ = @X"T or 


Hence 
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to assume 4——#1<0, where f is a real number. In this 
case the general solution of 

7 X" + PX = 0 

is 


X(x) = A cose + B sin Bs. 

To determine the arbitrary constant A and B, notice that, since 
from the first boundary condition 

0 = 4,4) = X(0) T() 
and T(t) is not identically zero, X(0) = 0. Similarly the second 
condition 

0 = wL,t) = XL) TH) 
implies X(L) — 0. From X(0) — 0 we conclude that 

X00) = Acos0+ Bsin0 = A=0, 


d hence that 
ae Beare X(s) = Bain Be. 


From the second condition 

X(L) = BsinBL = 0. 
But if B= 0, X(x) =0 and hence u(x,t) = 0, since 
u(x,t) = X(x) T(t). This contradicts the assumption that 
u(z,0) = fz) # 0. We conclude therefore that 


in BL = 0, 
from which a 
Ble, B=, wa 23 


Hence for every integral value of n we get a value of B and a dif- 
ferent solution, 
X,(2) = B, sin 7 


* There is no necessity for letting be negative. For suppose we had a term 


XG) = Cn sin 
in addition to bs 
X,(2) = By sin =, 
‘Then we could write their sum, say X,, where n is a positive integer as 


Ks) = Xqlo) + Xgl) = (Bq — Coq) sin 
—which is identical in form with X,(x). 
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Consider now the equation on 7(2), 
T” + PAT =0. 
Its general solution is 
T(t) = C cos Bat + D sin fat, 
and its derivative is 
T(t) = —BaC sin fat + BaD cos Bat. 
From the third boundary condition 


J ue0) = 





we conclude that 
X(x) TO) = 0, 


and hence that T(0) = 0. From this fact 
T'(0) = —BeC sin0 + BaD cos0 = BaD = 
or D = 0. Therefore 
T(t) = C cos Bat. 


Since B = nn/L, we have a solution T,(t) for each integral value 
of n: 





Tilt) = Cy cos 74, 


and for every integral value of n, 1 = 1, 2, 3, .... 
ee 


X,(*) Trt) = by sin Fx cos r (bn = BaCn) 


is a solution of the partial differential equation of Equation (129), 
which satisfies the first three boundary conditions of Equation (130) 
whatever the value of the constant b,. Since Equation (129) is 
linear, any linear combination of solutions is a solution and. hence 
the series : 
t) = > 5, sin = ae, 
u(x,t) x sin “7x cos “7 

is a function which formally satisfies the differential equation and 
three of the boundary conditions. There remains one more 
boundary condition to be satisfied, 


20) = fla). 
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To satisfy this condition the 4, must be determined such that 
#0) = fiz) = 2 bsin To (132) 


It is thus seen that the 5, are the sine Fourier coefficients of f(x), 
and hence are given by 
2 . 
b, = Zp fea) sin Fe de. 
‘The solution of the problem of the vibrating string with the given 
initial conditions is therefore 


u(s,t) = zz [J5-=) sin $x ee] sin Fx cos "74. (133) 


43. Alternate form for the solution of the vibrating string 
A suitable trigonometric identity transforms Equation (133) 
into a very elegant form. Remembering that 


sin A coe B= J [sin (A —B) + cin (A + 2), 
we may write 
sin FE coe Fe = 5 [ain FE ee — a + ain E(w +m], 
and Equation (133) becomes 
(3) = FE rein Zea + FZ sn (+42), (134) 
But 
fe) = $ basin Fs, 


and hence if we let + = x — at in the first sum of Equation (134), 
and x = x + at in the second sum, we may write 


ay = A= + fe bat) (35 
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In this formula we must remember that f(x + at) represents the 
periodic prolongation of f(x). Since at has the dimensions of a 
displacement, the function f(x — at) represents the function f(x) 
displaced at units to the right (see Figure 45). Similarly f(x + at) 





Fig. 45 


represents the function f(z) translated at units to the left, and we 
can therefore think of f(x + at) as representing a traveling wave 
of shape f(x) with velocity Fa. This is the basis for calling the 
partial differential equation 

a au 

G-*Bq-0 
(and its generalizations), the eave equation. 


44. An example of the wave equation 
We shall consider a numerical example of integration of the 
‘wave equation. Assume that the initial displacement of string 
is given by 3 
fo = 5” 
where ZL = 10 and a* = 10,000. [Wenote that f(0) = f(10) = 0] 


‘The Fourier series solution of this problem is given by Equa- 
tion (133), where the Fourier coefficients 6, are given by 


by x(10 — x) sin THe ee 


~ w/e 
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We leave it to the reader to verify that 


a =0 if nis even, 
b= <ty ifn is odd, 
= oe , 
and hence that 
wei Ay Sein Mx con iOest 
a ce nw 10 . 





Art. 44 


From what we have seen in the previous section we may write 


u(z,t) in the closed form 


ste) = AE = 1000 + fr + 1009 


The values of 1000 u(s,t) obtained by taking values of x at 
the points x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 = L and letting 
t = 0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06 seconds appear in the 


table below. 











[oj 0 0 
i838 6 2 0 
| 1s 12 4 0 
| 2% 7 6 0 
| 23 20 8 0 
; a 9 0 
B 20 8 o 
20 7 6 CY) 

1s 12 4 0 

; 8 | 6 2 0 
1 0 | o | 0 0 0 

i 






obbbbbdhle 


‘In computing these values one must remember that /(x) has been 
expanded into a sine series and hence f(—2) = —f(z), and that 


with a = 100 and 
t= w (m= 0,1, 2,...) 
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the tabular values are equal to 
Fife +m) + fle —n) 


‘The results are plotted in Figure 46. Note that the period of the 
vibration is 0.2 second, that is, the value of ¢ for which 


na 















































724 = 10mt = 2n. 
L 
10008 (1) 
4 4 
2a 
ts0ho2 
20 
‘al 1 -0f03 
12 
teoloa 
° 
4 
1-005 
s ze 3 a 
-4 
A ta 
Fig. 46 


45. ‘Temperature distribution in a slab 


As a second example of the application of the method of 
separation of variables we shall consider the two-dimensional 
Laplacian Vu = 0. As was mentioned in Section 41, the equa- 
tion 

au au 3 

Vu = Fat BO (136) 

might represent, for example, the steady-state temperature 
distribution in a plate insulated on its faces. Consider a rectangular 
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plate of sides a and 6 (Figure 47) on three sides of which the 
temperature is sssumed to be zero, while the temperature on the 
remaining side is a specified function of x, namely f(z). We are 

interested in the temperature 


1 u(zy) in the plate which satisfies 
Equation (136) and the boundary 
conditions 

(ob) 
1. u(y) = 0, 

. 2 Kay) = 0, 

(137) 
” 3. wed) = 0, 
oO a 
Fig. 47 4. wx,0) = fiz). 


As before, assume a solution of the form 
sey) = X(@) ¥(y) 
and substitute it in Equation (136), obtaining 


x"Y + XY" =0 
or 





‘The left-hand side of this equation is independent of y, and the 
right-hand side is independent of x; therefore their value is a 
constant, say A. Thus 
xX” AX =0, ¥Y" + AV =0. 
We shall assume 1 = —f* < 0 and show that a negative value 
of ) allows the satisfaction of all the boundary conditions. With 
) = —#F the general solution of 
X" + PX =0 
is 
X(s) = A cos Bx + Bein Br. 
From the first boundary condition 


0 = w(0,y) = X(0) YU), 
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and since Y(y) * 0, we find that 
XQ) =0, o XO) =A=0. 
Similarly the second boundary condition 
0 = may) = X(@) YO) 
implies X(a) = 0, or 
X(a) = BsinBa = 0. 


Since B = 0 would imply u(x,y) = 0, we must have sin a = 0, 
or 
Ba=m, B= =, a=1,2,... 


With these values of 8, 
X,(2) = By sin Ze, 
‘The solution of the equation for ¥(y) yields 
¥() = C cosh fy -+ Dsinh fy 
or, recalling that f is a function of m, 
Y,(y) = C, coh “Ly + D, sinh y, 
It is convenient to write Y,(y) in the equivalent form 
¥,) = By sinh = (y —y) 
and to notice that, since from the third boundary condition 
0 = w(=h) = X(s) ¥) 
and X(2) # 0, ¥(b) must vanish : 
¥,(6) = B, sinh “( —y) = 0. 
Hence 


Yo = 5, 
and the functions 


ws(s:y) = 6, sin 2 sink zy —b) 
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satisfy the partial differential equation of Equation (136) and the 
first three boundary conditions of Equation (137), whatever the 
constants 4,. Since Equation (136) is linear, any linear com- 
bination of the u,, such as 


ay) = =e sin te sinh Gy —s) 


is also a formal solution of Equation (136) satisfying three of the 
boundary conditions. In order that u(x,y) satisfy the fourth 
boundary condition, we must choose the b, such that 


u(z0) = fla) = 3 6, sinh (-2) sine 
“ 


and hence —4, sinh (amb/a) are the sine Fourier coefficients of 
fa): 


, sinh 2 = a fi fle) sin ae dz. 
‘The 5, are therefore given by 
~2 a gi Pt 
&n = a sinh (wnb]a) Je) sin "Te de, 
and the final solution of the problem reads 
1 a ana 
lear (enbja) S72) Oe 4s] 


sin “2x sinh zy —)). (138) 


u(y) = -23 


46. Example of the heat equation 
As an application of the solution of the previous section, assume 
fle) = T (a constant). 


Then 
St ea 
bn = a sinh (nnb]a) Josin Faas 
and 
b, = 0 Pa ieeress, 
AT 
' = evan tmas © iB Ode. 
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Hence by Equation (138), 
— 47 Fp _ sin (wn/a)x sinh [(w7/a) (y —5)} 
i (il) 
If the plate is square (a = 5) and we wish, for example, to find the 
temperature at the center, x = a/2, y = a/2, we obtain 


5 cient 
2, weak Py" 





etd. 


It can be shown that the sum of this series is 7/8 and hence that 
af ey a. 7. 
22) 0 4 
47. Voltage in a transmission line 


The voltage ¢(x,t) in a transmission line (of the submarine 
cable type), grounded at x = 0 and x = L, and with an 
initial voltage distribution f(x), was shown in Section 9 to satisfy 
the equation a ee ; 

mi Kae K =F (139) 


and the boundary conditions (see Figure 48) 





1. e(0,t) = 0, 
2. elt) = 0, (140) 
3. e{x,0) = f(z). 
(ot) (La) 
) t = 
Fig. 48 


If we assume a solution of the form 


ext) = X(x) TD), 
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Equation (139) yields 
Kx’'T = xP oo 72a Tad 
x KT 
Hence 
X" aX =0, PM — KT =0. 


If A> 0, the solutions of X” — AX = 0 are exponential in 
character, which is physically unreasonable. Hence we sssume 


A= 
and the equations above become 
X"+~PX=0, T+ PKT=0. 
‘The solution to the first of these equations is 
X(x) = A cos Bx + Bain Bx. 
From the first boundary condition, 
(0,1) = X(0) T(t) = 0, 





we conclude that 

X@) = 0, 
and hence 

XQ) = Acs 0 + Bsin0=A=0. 
From the second boundary condition, 
Lj) = XQ) TH = 0, 

we obtain 

XL) = 0 
and hence 


X(L) = BsinBl = 0, 
from which, as before, 
Alam, B=, n= 1,23... 
and 
X,(2) = Bain =. 
The solution of 
P+ PKT=0 i T= Core 
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T(t) = Cyan, 
by means of which, due to the linearity of Equation (139) we 
obtain the function 
st) = dasa (Fs) eomnes 
satisfying the partial differential equation and the first two bound- 
ary conditions. This function satisfies the third boundary con- 
dition if 
x0) = fix) = Bsn Fo 
Hence the , are the sine Fourier coefficients of f(z) : 
naif Sie) sin Te de. 


‘The complete sotution of the transmission line problem thus 
becomes 


zt) = ze Lf fe) sin “7 ae] ain ( Fn) inser, 


48. Example of the transmission line equation 


If we assume 
fie) = E (a constant), 
then 
a 2 fe Bain Me de 
= 0 ifmiseven 
= © itn is cae, 
and 
45 it 
a3 a (lL 
ton FR (Ele 
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‘The voltage approaches zero as t + © for any x due to the 
negative exponential. Also, for fixed ¢, the voltage is sinusoidal in 
nature. Both these phenomena agree with our physical intuition. 
While we are interested only in e(,¢) for x between 0 and L, since 
that is the length of the transmission line, the solution is valid for 
any value of x. The value of (x,t) at the point x» plus or minus 
integral multiples of 2 is the same as e( This is due, of 
course, to the fact that f(x) — E is actually a periodic function 
which coincides with E in the fundamental interval (0,L). 








49. Nonhomogeneous boundary conditions 
In each of the three previous examples, only one of the boundary 
conditions was different from zero, but this is not an essential 
restriction on the method of separation of variables, as will be 
shown by means of the following example. 
Consider the equation of the vibrating string, 
ae 1 Ow 
Bt at OF 
when the string is fixed at the end points but has an initial velocity 
(2) a8 well as an initial displacement f(x). The boundary con- 
ditions become in this case 


=, (141) 


1. ut) = 0 

2 Lt) = 0 

z $ ate) = a(x) (142) 
4. u(x,0) = fle). 


We sssume as in Section 42 a solution of the form 
uxt) = (2) TH) 
and obtain the two ordinary differential equations 
X" 4 BX =0, 7 4 par = 
‘The solution of the first of these equations is 
X(x) = A cos Bx + Bein Bx, 
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and from the first two boundary conditions of Equations (142) we 
conclude, as before, that A = 0 and 


a 
b= 2 = 

Since the third boundary condition does not eliminate one of the 
arbitrary constants appearing in the solution of T” + fta'T = 0, 
we have now the solutions 

g(t) = 8, sin cos Ft + 6, sin Fe sin “7s, 
which satisfy the partial differential equation of Equation (141) 
and the first two boundary conditions of Equations (142) whatever 
the values of 5, and cy. Hence the linear combination 
(at) = x ( sin Fx con Ft + cq sin Te sin =.) (143) 
also has this property. Computing the time derivative of this 
function, 


a na nm. nea re 
a Met) = & (40 Fo sin sin + 





and substituting it in the third boundary condition we obtain 
Sue) = we) = oe sin es, 


This equation shows that the constants (ma/L)c, are the sine 
Fourier coefficients of g(x), so that 


Ps 


Sn 





= ise) sin Fe de, 


With these values for the cq, u(2) as given by Equation (143) will 
satisfy the partial differential equation of Equation (141) as well as 
the first three boundary conditions of Equations (142). Simil- 
arly, substituting Equation (143) into the last condition, 


W(x,0) = (x) = ¥ b, sin “x, 
a 
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and we see that the b, arc the sine Fourier coefficients of f(x) : 
b= 2 fi payin eae 
Hence the function 
wees) = S| [Bf tenn Eee sin on 
+ [xe fee nin Pex a] sin Se in Pel (144) 


satisfies entirely the given differential system. 
As a numerical example, let f(z) and g(z) be given by the graphs 
of Figure 49 below, where it is assumed that L = 1. Then 


1004, = <2 fil?) sin wes de + ef, 


seo 
= it nin oad 


=0 if nis even, 


(2 —1) ain we de 


42 


100 Fh 





1004(x) 





Fig. 
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and 


; 


= 2 2 (2) sin nme de + 2 [1 (1 2) sin wm de 


= Oe it nin od 
=0 if nm is even. 
By virtue of these formulas, Equation (144) gives 


“0 = ae 2 gee 
sin mmx (= mat — a sin nat). (145) 


Another method of treating nonhomogeneous boundary con- 
ditions ia that of superposition of solutions. Suppose, for con- 
creteness, we have the partial differential equation of Equa- 
tion (141). Let 4(x,) be the solution to this equation with the 
boundary conditions 


1 gx = 0, 
2 gibt = 0, 

3 ZHe0 =0, 
4 a0) = fe. 


This problem has been solved completely in Section 42. Also, 
let ¥A,t) be the solution to Equation (141) with the boundary 
conditions 


1 Yo9 = 0 
2 Kb) = 0, 
2 FeO) = ate), 
4 x0) = 0. 


This differential system can also be solved by the methods of 
Section 42. Now, clearly the function 


x,t) = Hat) + et) 


Google 


12 SEPARATION OF VARIABLES Art. 50 


satisfies Equation (141) and the boundary conditions of Equa- 
tion (142). First u(2,¢) satisfies the differential equation, since 
both ¢ and y are solutions. Secondly, 

1. uO) = 40.) + YO.) = 0, 
since both 4(0,t) and 4(0,t) are zero by hypothesis. 

& u(L.t) = AL.) + WLt) = 0, 
since again ¢(Z,t) and yL,t) are zero by hypothesis. Also, 

3. Sao) = 420 + 2 ue0) = a2), 
since d¢(x0)/at = 0 and oy{x,0)/et = g(x). Finally, 


4. -ux,0) = f(x,0) + (x0) = f(x) 
since $(x,0) = f(x) and ¥(x,0) = 0. 


50. Temperature distribution in a semicircular plate 

‘The problems of the previous sections were all governed by 
partial differential equations with constant coefficients, and hence 
yielded ordinary differential equations with constant coefficients. 
In this section we shall consider an example which leads to an 
equation with variable coefficients (the so-called Cauchy equation). 
Other variable coefficient equations will be studied in Chapter V. 


y 
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We wish to find steady-state temperature distribution in a 
semicircular plate of radius a, insulated on both faces, with its 
curved boundary kept at a constant temperature U, and its bound- 
ing diameter kept at zero temperature (see Figure 50). The partial 
differential equation governing this phenomenon was found in 
Section 8 to be the Laplacian equation 

fu | ou 
ae + oF = 





Since the cartesian form of Laplace’s equation is not well adapted 
to deal with the problem in hand, we shall transform it into its 
polar form. 

‘The relation between the cartesian coordinates x, y, and the 
polar coordinates r, 6 is 


x=ros0, y= rsind, 








or, inversely, 
r= Vey y¥, = are tan 2 
We must compute 2%u/ax? and éu/ay*. Towards this end we find 
am 1 et 
wre 
#0 _ 2xy 
xt ~ yt" 


and since 
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-3 ("3 


and Vu = 0 is therefore 
a ou Ou 
ott em 
in polar coordinates. 

The temperature » considered as a function of r and @ (see 
Figure 50) satisfies Equation (146) and the boundary conditions 
La) = Us 

2. fr,0) 
3 arin) = 0. 
‘The method of separation of variables suggests assuming a 
solution of Equation (146) of the form 
Mr) = RG) 8). 
Substituting this into Equation (146) reduces it to 
AR’ + rR'O + RO” = 0. 
Dividing by RO and separating the variables, we obtain 
e e 


Rr 
Rt -> (147) 


=0 (146) 






a 
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‘Since the right-hand side is a function of @ only and the left-hand 
side a function of r only, their common value must be a constant, 
say A. Since for A < 0, the solutions of @”+A® = 0 are 
exponential in character, which is not physically reasonable, we 
assume A = f*>0. Equation (147) then yields the two 
ordinary differential equations 
PR’ + 1R — BR =0, 
6” + pO = 0. 

‘The general solution of the second is 

0) = A cos) + B sin 68. (149) 
The first equation is a so-called Cauchy equation. In order to 
solve this equation we make the transformation 


r=e, 


(14s) 


and the first of Equations (148) reduces to 
@R 
Ga — PR =0. 
‘The general solution of this equation is 
R = Cel + De, 
or sincer = e’, 
R = Cr? + Dr?. (150) 
‘The arbitrary constants appearing in Equations (149) and (150) are 
determined by the boundary conditions. From the second 
boundary condition, * 
1,0) = R(r) @(0) = 0, 
we conclude that @(0) = 0 and hence that A = 0. From the 
third boundary condition, 
1,7) = R(r) On) = 0, 
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we conclude that (nr) = 0, or 
O(n) = Bsin Br = 0. 


Since B #0 [otherwise u(r,8) = 0] we must have sin Bx =0, 
from which 


Puan of Bm n=1,2,3.0. 
Hence we find the solutions 
0,(6) = By sin 28. 


Examining Equation (150) we see that as r+ 0, the term r-? +00 
and hence makes u(7,6) arbitrarily large for r very small. This 
contradicts the physical problem, since at r = 0, u(r) = 0. 
Therefore D must be equal to zero. The functions 


ug(7,0) = byr® sin nf 
satisfy the partial differential equation and the last two boundary 
conditions. As before, any linear combination 


lr.8) = Sb sin. nd (si) 
=” 


also has this property whatever the values of the constants 5,. 
The first boundary condition is also satisfied if 


(a0) = Up = 3 ba" sin 06, 


a 

that is, the a"b, are the sine Fourier coefficients of Uy, 
ob, = 2f U, sin nd db 

4U, 

- 





if n is odd, 


=6 if nm is even. 
Hence 
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if n is odd, and is zero otherwise. The series of Equation (151) 
then becomes 


u(r,8) = Me z a (zy sin 76, (152) 


net! 





which is the desired solution. 


51, Transverse vibrations of a beam 


As an example of an equation of order higher than the second, 
consider the transverse vibrations of a simply supported beam 
which has zero initial velocity and whose initial displacement is 
fz) = Aa{L — 2), where A is a constant and L the length of the 
beam. The equation governing this phenomenon was found to 
be (Section 6) 


a + 2% =0 (# = Ee), (153) 
say with the boundary conditions : 
1. (04) = 0, 
2. Lt) = 0, 


oe 
3. WO.) = 0, 

a (154 
4 ga Lt) = 





5. & x20) = 0, 
6. 90) = Axl — x). 
Assume a solution of the form 
Wat) = Xe) TO, 
which transforms Equation (153) into the equation 


xe 7 


XT" + a@X*T=0, or X¥ = ar 
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‘The constant value of the two sides of this equation must be 
assumed positive from physical considerations, namely, to make 
T(t) oscillatory. If we call it 6*, we obtain the two ordinary dif- 
ferential equations 

xe —PX=0, T+ oBT =0. 
‘The general solution of the X equation is 
X(z) = AsinhBe + Bosh Be + C sin Be + D cos Br. (155) 
‘The first boundary condition, 


HOt) = X(0) T(t) = 0, 
implies that X(0) = 0 and hence that 
X0) — B+ D=0. (1560) 


‘The second boundary condition, 
Hs) = XL) T) = 0, 

implies that X(L) = 0, and hence that 

X(L) = AsinhBL + BooshBL + CsinBL + DowBl =0. (1566) 


To utilize the third and fourth boundary conditions, we calculate 
X"(z) from Equation (155) : 


X°(a) = AB? sinh Be + BBY cosh Be — CB sin Bs — DB* con Be. 
and find that the third boundary condition, 
F709 = XO) TH = 0, 
requires that X"(0) vanish. Hence 
X"0) = BB — DB = 0. (1566) 
‘The fourth boundary condition, 
Ato = x OM = 0, 
requires that X’(L) vanish, and hence 
X"(L) = ABP sinh BL + BB? cosh BL 
— CP sin BL — DBP cooBL = 0. (1560) 
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Equations (156a,b,c,) are four linear algebraic equations on the 
constants A, B,C, D, and B. From Equations (156a,c) we have 
B+D=0, B-D=0, 
and hence 
B=D=0. 

From Equations (156b,4) : 

Asinh PL + Coin PL = 0, 

Asinh AL — CsinfL = 0, 
and hence 

AsinhfL = 0, CsinfL = 

The first condition gives A = 0 [since sinh BL = 0 implies 
BL = 0, that is, B = 0 and hence X(z) = 0]. If C is chosen 
equal to zero, then X(x) is identically zero. Hence we assume that 
sin BL = 0, or 

AL = nr, 2 = 1,2,3,... 
and obtain the solutions 


X,(2) = C, sin a 





The ordinary differential equation 
T’ + aT =0 
has the general solution 
T(t) = Ecos oft + F sin of. 
‘The fifth boundary condition, 
$00) = x) TO = 
requires that (0) — 0, that is, that 
TO) = —ofME sin af"0 + af'F cos of) = af'F = 0. 
Hence F = 0 and the solutions 7,(t) become 





T(t) = Eco Et, 
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‘Therefore any function yq(%t), 

Inlet) = X,(x) Tilt) = b, sin t con a BT, 
where 5, = C,E, satisfies the partial differential equation of 
Equation (153) and the first five boundary conditions of Equa- 
tion (154), whatever the values of the constants 6,. The linear 
combination of y,, 


yet) = z 6, sin Fe cosa , 


also has this property, and will satisfy the sixth boundary con- 
dition if 


#0) = AL —x) = > 2, sin Fs. 
a 





Hence the 5, are the sine Fourier coefficients of Ax(L — 





27 _ ter 
b, = 7 ff Al —») sin Fe ae 
salt 
=a if nis odd, 
=0 if n is even. 


The solution to the boundary value problem of Equations (153) 
and (154) is therefore 





a Jsin x cosa. (157) 


52. The vibrating membrane 
As an example of a problem involving three independent vatia- 
bles, consider the equation of the vibrating membrane (Section 4) 


re ee 


Bsa (158) 
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We shall assume that the membrane is rectangular in shape (see 
Figure 51) and satisfies the following boundary conditions : 





1 xy.) = 0, 

2 w(a,y,t) 

3. u(z,0,t) 

4. w(sb,t) = 0, 

5. 2 atay.0) = 0, 


6. u(xy,0) = Tay(x — a) (y —5). 


Fig. 51 
‘The first four boundary conditions state that the membrane is 
held fixed on the boundary for all-time. The fifth states that the 
membrane has no initial velocity, while the sixth gives the initial 
displacement. 

‘As customary in the method of separation of variables, assume a 
solution of the form 

u(x,9,t) = X(x) ¥(9) Te). 
Substituting it in Equation (158) and dividing by XYT, we obtain 
SE a (159) 


‘The right-hand side of this equation is independent of x, and the 
left-hand side is independent of ¢ and y. Their common value 
must be a constant A. Further, from physical considerations, 
A must be negative, say A = —f*, for otherwise X(x) would vary 


9 
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exponentially and the boundary conditions could not be satisfied. 
Equation (159) thus gives the two separate equations : 


xX” + px =0 (160) 
and aps 
yrapte (161) 


In Equation (161) the right-hand side is independent of y, while 
the left-hand side is independent of t; hence their common value 
must be a constant, say ». As before, we let « = —a%, and 
hence Equation (161) gives the two separate equations 
"4+ e8¥ =0 (162) 
and 
T’ + Aat + PT = 0. (163) 
‘The partial differential equation of Equation (158) has thus been 
replaced by three ordinary differential equations, Equations (160), 
(162), and (163). The solutions of these three equations are, 


respectively, 
X(2) = Aconps + Bein Bs, 





¥(y) = Ccosay + Dainay, (164) 
Nit) = Ecosyt + Fainyt, 
where 
YP = det +P). 


From the first boundary condition, 
™0,y.t) = X0) YO) TH) = 0, 


we conclude that X(0) = 0 and hence that A = 0. From the 


second, 
Mays) = X() YO) TW) = 0, 
we conclude that X(a) = 0 and hence that Bsinfa = 0. If 
# is not to be identically zero, we must have sin fa = 0, or 
fa =m, n= 1,23... (165) 
Similarly, the third and fourth boundary conditions imply that 
C = 0 and that 
ob = mn, m= 1,2,3,... (166) 
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Finally, from the fifth boundary condition, 
J x90) = Xe) YO) TO) = 0, 
we conclude that 
TO) = —yEsin0 + yFcos0 = +F =0 

and hence that F = 0. 

Indicating the value of BDE corresponding to the integers n and 
min Equations (165) and (166) by da, we may write 

er an \* ‘wer \? 

sele98) = Anta Ein yon [er V() + (F))]- 
‘These functions satisfy the partial differential equation and the 
first five boundary conditions for any choice of the constants dy_, 


and for any integral values of n and m. In order to satisfy the 
sixth boundary condition we assume for w the linear combination 


Heyt) = ZY Henle 
Letting t = 0 in this series, we see that 
u(xy,0) = s 3 am tin So vin Fy 


and hence the a,,, are the sine Fourier coefficients of the function 
u(x,y0). By Section 38, 


a Sh feina Tempe 
= AE fate 0) sin de fay — 0) sin Ey ay 
4T 2a a 
= Fe ay — Bom - 0 
— Te it a and mare odd, 


=0 otherwise. 
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Our final solution is therefore 


wey = ATS 








a8... 
me, 


sin ae sin “Fy co0yants (167) 


me (+B 


53. Nonhomogeneous equations 

‘The partial differential equations considered in the previous 
examples were all homogeneous. In this section we shall consider 
a certain class of nonhomogeneous equations, namely, those in 
which the function on the right-hand side of the equation is a 
function of one variable only (the so-called forcing function). 
Many physical problems give rise to equations that fit into this 
category, of which a few will be mentioned. 


1. Radioactive Decay. We have seen in Section 7 that the 
temperature distribution in a bar of length L satisfies the one- 
dimensional heat equation 

me 

Ox* 
If now we assume that heat is also gained by radioactive decay 
of the material in the bar (see Problem 20, Chapter I) an additional 
term of the form Ae-* must be added, and the above equation 
modified to read 


where 








a 


ou au 25 
Fe tS = -Aem, (168) 


‘Typical boundary conditions associated with this equation might 
be i 





1. 0(0,t) = 0, 
2. u(L,t) = 0, (169) 
3. u(,0) = fl), 


where f(z) is a given function of x. 
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2. Decomposition. Consider again the one-dimensional flow 
of heat and assume that heat is also gained by the decomposition of 
the material in the bar (see Problem 21, Chapter I), This heat is 
generally assumed to be generated at a constant rate A, and the 
heat equation therefore assumes the form 


a au 
Ba -G--4 (170) 
Appropriate boundary conditions might be 
1. 4(0,t) = 0, 
2 WL) = 0, 0) 
3. Wx,0) = 0. 


3. Forced Vibrations. Consider the equation of the vibrating 


string, 

ing, las 
of” 

and assume that an external force normal to the string is applied 
which varies, say linearly, with the distance from the end x = 0 
(see Problem 22, Chapter I). In this case the equation assumes the 
form 





au 
on 


Ge ~ OG = Ae (172) 
Boundary conditions for such a problem might be 
1. x0,t) = 0, 





2 
3. u(z,0) = 0, 


a 
4 By M(s0) = 0. 


(173) 


Note that here, as in the second example, all the boundary con- 
ditions are homogeneous, yet there exists a nontrivial solution, 
since the differential equation is nonhomogeneous. [If both the 
equation and the boundary conditions are homogeneous, the only 
solution is the trivial one, u(x,) = 0.) 

The technique used in solving equations of the type given 
above consists in introducing a new function of the form 
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w(xt) + 9x), where w(x,t) satisfies a homogeneous partial 
differential equation and (x) is determined from the solution of an 
ordinary differential equation. We shall illustrate this method 
in the next section. 


54. The beat equation with radioactive decay 

To start the solution of Equation (168) we introduce a new 
function w{z,t) defined by the equation 

u(x,t) = (xt) + 2), (174) 

where (x) is a function of x only, as yet undetermined. We then 
try to determine (x) 20 that (x,t) satisfies the homogeneous 
equation obtained from Equation (168) by setting A = 0. 
Substituting Equation (174) into Equation (168) we obtain 


a + ¥"(%) - oe = —Ae-*, 
Hence, if 9) is such that 
¥'@) = Aes, (175) 


t=(x,¢) will satisfy the homogeneous equation 





Equation (175) is readily integrated and gives 
Wa) = —4ewt Bet, (176) 


where B and C are arbitrary constants. Now any value may be 
used for B and C, and (x,t) will still satisfy a homogeneous 
equation. At first we might be tempted to let B and C be zero. 
A little reflection will indicate, that, since we have the arbitrariness 
of B and C at our disposal, it might be well to determine them 
such that the boundary conditions assume as simple a form as 
possible. From the firet of Equations (169) and Equation (174) 
we see that 
u(0.t) = w(t) + (0) = 0, 


(0,1) = —K0). 


or 
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Similarly, from the second boundary condition, we obtain 
MLt) = W(L.t) + WL) = 0, 


oe 
Ls) = AKL). 
‘The constants B and C will be so chosen that 
HO) = 0 = HL). 
‘Thus we obtain from Equation (176) 


A A 
W)=-S+ C= © CHS 
Also, 
WL) = —4qt4 r+ 4-0, 
or 
1(A4 A 
8-7 (4-+-4). am 


With these values of B and C, the function e(s,t) satisfies the 
homogeneous equation 





ao eo 
at 
and the boundary conditions 
1 0) = 0, 
2. WL t) = 0, 


3. wfx,0) = fiz) — Wx) 
= fe) + 4 [fet — 0-1]. 
‘Once this system is solved by the techniques developed in previous 


sections, the original function u(s,t) is readily determined from the 
equation 


He) = mat) + Ys) 
= ost) + 4 [- + pct N+ i}. (178) 
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55. The general case 

A rather wide clase of homogeneous linear partial differential 
equations of the second order in two independent variables is 
solvable by the method of separation of variables—namely, the 
class of equations which can be put into the form : 


Ne) Fe + 60) Fe + 
+ HO) E+ Ue) + hoe = 0. (179) 


‘This family includes as special cases many of the equations con- 
sidered thus far. If we assume a solution of the form 

uy) = X(#) ¥(y) 
and substitute it in Equation (179), we obtain 


AX"Y + G:X¥" + fX'V + XY’ + (f+ d)X¥ = 0, 


YAK" + faX’ + feX] + ATHY” + WY’ + $Y] = 0. 
Dividing through by XY, this equation reduces by the usual 
arguments to the two ordinary differential equations 

SX" + fAX' + fk = AX, 

OY" + bY’ + bY = —AY, 
where A is an undetermined constant. The solution of these 
ordinary differential equations is, in general, not a trivial problem, 
but the separation of the variables has been achieved. 


or 





56. A trial function 

Another elementary method frequently useful in solving partial 
differential equations is to assume a solution of a particular form 
involving a number of arbitrary constants and/or functions and then. 
try to determine these constants and/or functions so as to satisfy 
the given differential system. We shall illustrate this method 
by means of an example. Consider first the familiar equation of 
the vibrating string, ou ow 


=e (180) 
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with the boundary conditions 


L u(0,t) = 0, 

2 Lt) = 0, 

3. x0) = Dein Fx, 
4 a) =0 


A trial function involving sines and cosines seems appropriate 
since only second derivatives are involved in Equation (180); 


(@ sin x)/dx* — —sin x, (d¥ cos x)/ds* — — cos x; 
and furthermore the boundary conditions involve a trigonometric 
function. Let us therefore assume a solution of the form 

u(x,t) = A sin ax cos pt, (181) 
where A, «, and f are constants as yet unspecified. Substitution 
of Equation (181) into Equation (180) yields 

—a'a2A sin ax cos Bt = BA sin ax cos ft, 

and shows that Equation (181) is a solution of Equation (180) if 


att — f, 
that is, if 
uz) = A sin ax cos aat, 

where A and a are still arbitrary. ‘This function can be made to 
satisfy the boundary conditions by choosing A and a properly. 
In fact, the first boundary condition is automatically satisfied, 

u(0,t) = A sin 0 cos aat = 0. 
The second condition gives 

UL,t) = Asin al cos aut = 0, 


or 


aa neh. 


The third boundary condition yields the equation 


na nt 
0) = Asin "Te — Dain Fs, 
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which becomes an identity if 
n=1 ad A= 





Hence 


u(x,t) = Dein prs Ft (182) 


is the desired solution, since it satisfies entirely the differential 
system (the fourth boundary condition being automatically 
satisfied). In general, if a more complicated trial solution than 
Equation (181) had been assumed, the additional terms would 
have dropped out. 


57. A trial function method based on Fourier series 


A variation of the trial function method consists in assuming 
the solution of a partial differential system can be written in 
the form of a Fourier series. This method might be called the 
“‘ direct Fourier series method ” and is analogous to the methods 
used in solving ordinary differential equations by power series. 
To illustrate this approach we shall consider the heat equation in 
two dimensions. If we have a two dimensional plate, the equation 
satisfied by the temperature u was found to be (Section 8) 

mu tu le 

eta tar (183) 
We shall assume that the plate (see Figure 47) is insulated on its 
lateral surfaces, has its four edges kept at zero temperature, and 
has a temperature T (a constant) at t = 0. The boundary 
conditions for this problem may therefore be written 


1. Oy.) = 0, 
2 way!) 
3. a(s,0,t) 

4. u(s,b,t) = 0, 

5. u(ay,0) = T. 
Let us prolong the function u(x,y,0) as an odd function in both 
the and y variables so that the temperature distribution u(+,y,t) 





(184) 
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will also be an odd function of x and y for all time t. That is, 
we shall assume 1(2,y,t) can be written in the form 


arnt) = Ze ZS sin resin Ey, (185) 
‘ 
where the coefficients f,., of the double Fourier series are functions 
of time. We see that if u(x,y,t) is so defined, it automatically 
identically satisfies the first four boundary conditions of Equa- 
tions (184). 
If we substitute a typical term of Equation (185), 


Fam(t) sin = sin *y 
into Equation (183), we have 
mT ‘mn \* a. mer 
(Fyn a toy 





1 ._ aw. mer 
= alm sin a sin Hy. 
If this equation is to be identically satisfied we must have 
: 
45) (F) 
SF nm(t) + Bam fam(t) = 0, (186) 
sap 
ne= (I+ (FY 
The solution of Equation (186) is 


Fam(8) = Ane eoom's, 





le 
bm = hfe 
6 


where 


and hence 
o(,9,t) = SS Aneto sin x sin My (187) 


is a function which satisfies the partial differential equation of 
Equation (183) and the first four boundary conditions of Equa- 
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tion (184). If u is to satisfy the fifth boundary condition, we must 
have 

oe fr. ater 

(2,90) = T= a> Aum sin x sin Sy, 


Hence the Aya are the Fourier sine coefficients of the double 
Fourier series of u(x,y,0) : 








4 fe . ‘. 
Aan = 55 J [pT sin “Poe sin Fy ae ay 
AST ifm and m are odd, 





=0 otherwise. 
The solution to the heat problem of Equations (183) and (184) is 
therefore 


ways) = OP 





Lig Be gig, Me att 
<— sin 2x sin FE yeetonmee 
as , 2 6 . 


na= (5) + (F) 


In conclusion we might say that while the trial function method 
by Fourier series is generally simpler than the method of separation 
of variables, it is of limited application, due to the fact that the 
assumed solution must automatically satisfy one or more boundary 
conditions. 


(188) 
where 


PROBLEMS 
1, Verify the values of the coefficients of the Fourier series in the 
example of the vibrating string, Section 44. 


2. Figure 52 shows the initial displacement of a perfectly flexible 
string of weight 1 Ib per ft, fixed at the end points and under a tension 
of T lb. Attime t = 0, the string is released from this position, 


(@) Write the expression for the deflection u(x, t) 


Google 


(b) Assuming T = 
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10 Ib, wm = 0.0322 Ib/ft, g = 32.2 fi/sect, L = 


10 ft, a = 1/3, plot the graph of u(x,t) at t = n/100, where » = 1, 
2, 3, 4, 5, dividing the interval (0, L) into ten equal parts. 


2 
ub, or 3-uls,0} 





Fig. 52 


3. Solve the problem of the vibrating string for the boundary 


conditions 


1 
2 


3. 
4 


0(0,1) = 0, 
Lt) = 0, 
2 u(x,0) = 0, 


4x0) = A(L sinh x — x sinh L). 


(Use the constants 7, w, g, L of Problem 2.) 


4. Find the displacement u(x,t) if the 





displacement of a 


vibeating string ia zero, but the initial velocity is given by the function 
of Figure 52, Assume a = 0.5. 


5. Solve the problem of the vibrating string for the boundary 


conditions 


1. 0,t) = 0, 
2 ML) = 0, 
3 a u(x,0) = AxL — »), 
4. (x0) = 0. 
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6. Solve the problem of the vibrating string for the following boundary 
conditions : 


L (0,1) = 0, 

2 aL = 0, 

3 Fux) = ae 

4 a0) = Gy OO< <d 





» <r 
Assume a* = 2500 and L = 1. 109 

7. Solve the problem of the vibrating string if the string (of length L) 
is fixed at the end points, and has an initial velocity given by Figure 52 
(« = }) and an initial displacement given by Figure 53. 





the change of variables 
zs=€tn yot—on 
and obtain an explicit integral representation for u(x, y) subject to the 
boundary conditions 
(20) = fiz), Z ute) = ate). 


9. Deduce the formula derived in Problem 8 by a consideration of 
the Fourier series of Equation (144), Section 49. 


10. Solve Problem 3 of Chapter I. 
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11. We have seen in Problem 11 of Chapter I that the angular 
displacement @ of a shaft vibrating torsionally satisfies the equation 
20 


we 
Ox** 
where a = Gg/8, G being the shear modulus, 5 the density, and 
g the gravitational constant. Find the twist 6(,t) of a steel ahaft 10 ft 
long, of specific gravity 8,and shear modulus 12 x 10* psi if it is fixed 
at both ends and given an initial twist (=, 0) = “°C —=—*) and initial 
velocity zero. 

12. Attempt the solution of Problem 17 of Chapter I for the longitu- 
dinal displacement of the bar. 


=@é 








13. Verify the values of the Fourier coefficients derived in solving 
the heat equation of Section 46. 


14. Find the steady-state temperature distribution in a rectangular 
slab of sides a and 5, insulated on the lateral surfaces and satisfying 
the boundary conditions 


1. Oy) = 0, 
2 way) = 0, 
3. ux) = 0, 
4. wfx0) = Tafa — x). 





15. Solve the steady-state heat equation for a slab in the form of 
a rectangle (see Figure 47) if the side along the x-axis has 2 temperature 
distribution of the form of Figure 53 (L = a), the remaining three sides 
are kept at zero temperature, and the lateral faces of the slab are insulated. 


16. A rectangular plate of dimension a units along the s-axis and 
5 units along the y-axis has two of its edges x =0 and x =a insulated. 
The edge y = 0 is kept at a constant temperature 7, and the edge 
y = b is kept at zero temperature. Find the steady-state temperature 
distribution in the plate. Assume that the lateral faces are insulated. 


17. Find the steady-state temperature distribution «x, y) in 2 semi- 
infinite strip whose faces are insulated, whose edges x = 0 and x =a 
are kept at zero temperature, and whose base y = 0 is kept at a constant 
temperature T, 
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18. Solve Problem 17 for the case in which the edges instead of being 
kept at zero temperature, are insulated. 


19. Solve Problem 17 for the boundary conditions in which the edge 
x = 0 is kept at zero temperature, the edge x = a is insulated, and 
the base is kept at a constant temperature T. 


20. The constant A in Section 45 was assumed negative. Try carrying 
out the solution if A= a? > 0. 


21, A bar of length Z has an initial temperature distribution as 
shown in Figure 54. From t = 0 on, the ends of the bar are kept at 
zero temperature. Find the equation for the temperature u(z, 2) of 
the bar. Assume the surface of the bar is insulated. 


u(xo) 





L 
2 
Fig. 54 
22, Find the temperature distribution s(+,t) in the bar of Problem 21 
if it has an initial temperature distribution as given in Figure 53. 
23. Solve the heat equation for one-dimensional transient flow with 
the boundary conditions 
1. uO) = 0, 
2. wL,t) = 0, 
3. u(z0) = aL — 2). 
24. Solve the heat equation 





subject to the boundary conditions 
1. (0) = A (a constant), 
2. (Lt) = B (a constant), 
3. u(#0) = a(L — 2). 
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25. Solve the one-dimensional heat equation 
eu yu 
Saak 
Ox* a 


with the adiabatic boundary conditions 
a 
1 gM) = 0, 





a 
2 Zul.) 
3. u(x0) = x. 

26. Solve the submarine cable equations for the current and voltage 


in a cable of length / if at x =0 and x =/ the cable is short-circuited, 
while the initial current distribution i(,0) is =(!—). 


27. A high-frequency transmission line of length 1 is grounded at 
x = 1 and open-circuited at # — 0. If the initial current and voltage 
distributions are i(e,0) =I, sin (m/l) # and ex,0)—=K, (sinh x — 
x sinh J), respectively, determine the current and voltage in the line at 
any time t. 

28. A semicircular plate of radius A has its circumference kept at 
a temperature u(A, 6) = T0r—6), while the bounding diameter is 
kept st zero temperature. Find the temperature distribution 1(r, 8). 
Assume the lateral surfaces of the disk to be insulated. 

29. A ring-shaped plate (eee Figure 55) of inner.radius @ and outer 
radius 6 is insulated on its lateral surfaces. Find the steady-state 
temperature u(r, 6) if the initial temperature distribution on the inner 
circle is A6(2n—@), and the initial temperature distribution on the 
outer circle is BO* (2n —6). 

30. Investigate the solution of Problem 29 as the inner radius @ 
approaches zero. 


31. A beam of length L is simply supported at x = 0, x = L, 
and is acted upon by two bending moments Mg equal and opposite 
in sign at these end points, The deflection due to these moments is 

M 
sal — 9). 


Attime t =0, the moments are suddenly removed. Write the equation 
for the transverse motion y(x, t) of the beam. 


10 
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Fig. 55 


32, A beam of length L is simply supported. If its initial velocity 
is zero and its initial displacement is D,+ sin m2/L,find the equation of 
its transverse vibrations. 

33. Solve the example of Section 57 by the method of separation 
of variables. 

34. Solve Problem 49, using the method of separation of variables. 

35. Solve Problem 50, using the method of separation of variables. 

36. Find the steady-state temperature u(x, y,3) in a rectangular 
parallelepiped (see Figure 56), if the face in the xy-plane is kept at 
‘a temperature given by 

(2,7,0) = Ugey(a — 2) (b —y), 
and the remaining faces are kept at zero temperature. 

37. Solve the heat equation with radioactive decay, 


ou ou 
OG tA = 
with the boundary conditions 
1 40) = 0, 
2 WL.) = 0, 
3. x0) = T (a constant). 
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(0,80) 
y 
Fig. 56 


38. Solve Problem 37 with completely homogeneous boundary 
conditions. 


39. Solve the heat equation with decomposition, 


ga Rt 48 
with the boundary conditions 

1 u(0,t) = 0, 

2 (Lt) 

3. wfx,0) 





40. Determine the general form of the solution of the following 
boundary value problem : 


Baek + fe); 
1 Ot) = a, 
2. wLt) = 5, 
3. (x0) = fs). 


41, Investigate the solution of Problem 40 where f(x) is replaced by 
an arbitrary function of time, $42). 


42. Solve the equation for the forced vibrations of a string, 


au au 
Rao + as, 
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with the boundary conditions 


1. 10,2) = 0, 

2. w(L,t) = 0, 

3. u(z,0) = 0, 

4. 2u(m0) = 0. 
ot 


43, Solve Problem 19 of Chapter I for the temperature in the thin 
bar, Assume that the ends of the bar are kept at the temperature of the 
surroundings, 7, and that the initial temperature u(x, 0) = Te-*, 

44, Solve the equation of the vibrating string (of length L) in which 
the initial velocity is given by 

Vy sin? F2, 

and the initial displacement is zero. If the string is 3 ft long, weighs 
O.1 Ib, is subjected to a constant tension of 5 Ib, and V» = 1.5 fps, 
calculate the displacement of the mid-point of the string at a time 0.01 sec 
after the string was given its initial velocity. 

45. Solve the equation of the vibrating string in which the initial 
displacement is given by 

D, (sin Tx + 4 sin 2), 

and the initial velocity is zero. 

46. Solve the heat equation 





aon 
Ox Ot 
for the boundary conditions 
1 ut) = 0, 
2. w(t) = 0, 


3. ulx,0) = T sin ee 


47, Solve Problem 28 of Chapter I by assuming 2 trial function 
in the form of a polynomial. 


48. Solve the illustrative example of Section 52 by the trial function 
method with a Fourier series expansion. 
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49. Let a rectangular membrane, fixed on its boundary (Figure 51) 
have an initial displacement and an initial velocity of 
u(x,9,0) = Asty (a — x) (6 —y), 
@ 
Hy 9.0) = 0, 


respectively. Find the coefficients dam of the double Fourier series 
expansion for u(, y, t), using the trial function method with a Fourier 
series expansion. Plot the displacement of the center of the membrane 
asa function of time for ¢ = 3,6 = 2, A =2. 


50. Solve Problem 49 with the boundary conditions 
u(x,y,0) = 0, 
é@ 
$e,3.0) = Ant (a —2) (6-9). 
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Chapter IV 
THE FOURIER INTEGRAL 


58. A heat flow problem 

Consider an infinitely long thin bar insulated on its lateral surface 
(sce Figure 57). We wish to determine the temperature u(:,t) 
if the initial temperature is given by 


ux0) = fz) me <x<o, (189) 
The equation satisfied by u(x,t) was shown in Section 7 to be 
ou au 

iS (190) 


and this problem appears similar to those treated in Chapter III. 
‘There is, however, one important difference : the bar is infinite 
in extent. If we try to solve this problem by the method of 
separation of variables we shall see that Fourier series are not 
adequate to yield a complete solution. 


Tt) 






ulxo)m #00 


Fig. 57 


142 
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Assume a solution of Equation (190) of the form 


ust) = Xs) TE, 
which when substituted in' Equation (190) yields 
r xf 


XT = AX"T, or a —, 


where f is a real constant. The ordinary differential equations 
resulting from this equation are 


X"+—pxX=0, "4+ aPT=0, 


and their general solutions are 
X(z) = A’ cosBx + BY sin Bx, 
T(t) = Cees, 

so that 


u(x,t) = X(x) T(t) = Ae" cos Be + Be~*#™ sin Bx, (191) 
where A = A'C, B = BC. 

In the case of a bar of finite length, the values of 8 were 
determined by one of the boundary conditions at the end of the 
bar. Here the bar is infinite in length and there are no such 
boundary conditions. Hence 8 cannot be determined as was done 
in the case of a finite bar. Furthermore, since the initial function 
fla) = ulz,0) extends from —© to +, it is not necessarily periodic 
and hence may not have a Fourier series expansion. 

From the above remarks we see that if this and similar problems 
are to be solved, the methods of Fourier series must be extended. 
We shall undertake this extension in the next section in a systematic 
fashion and thus define the Fourier integral. Before undertaking 
this treatment we can carry the solution of the present problem 
one step further. 

The function u given by Equation (191) satisfies Equation (190) 
for any real value of B. Hence we may write an infinite, continuous 
set of solutions 


ugx,t) = A(B)e-*8* cos Bx + B(B)e-*** sin Bx, (192) 
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where B may take on any real (not necessarily integral) value, 
and A and B may be different for each value of , that is, A and B 
are functions of B. Since Equation (192) satisfies the linear partial 
differential equation of Equation (190) for any values of f, the 
linear combination of an infinity of functions, up, 


wet) = K f° [ABle-#* cos Bs + BB“ sin Be] dB, (193) 
where K is an arbitrary constant, will also satisfy Equation (190).* 


59. The sine and cosine transforms 


Given a function f(t), not necessarily periodic, extending from 
—© to +» (see Figure 58), we wish to obtain for it an expression 
analogous to the Fourier series expansion, but valid between —2 
and +=. This can be done as follows : 


+n) 


Fig. 58 
Define a function fr(t) such that 


fit)=f) -Z<e<F, 
Frit) = 0 otherwise, 
* There is no necessity for extending the range of integration to —°, for calling 
A‘(f) and B’(8) the arbitrary functions A(—f) and —B(—8), respectively, 
ACT cos (—B)x = A’(B)e<*#* cos As, 


asd BUA)" sin (fz = BY f)er# sin fax; 
that is, these terms are of the same form as the terms appearing in the integral 
of Equation (193). 
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(see Figure 59). If the graph of fr is repeated every T units 
we obtain a periodic function of period 7 (see Figure 60), and 
presumably, if 7 is increased without limit, the Fourier series 


fr(t) 


Fig. 59 


of fr(t) will approach that of f(t). Since the function fr(t) has 
been prolonged as a periodic function, its Fourier series expansion 
is given by Equations (90) and (91) where 2L = T: 


1 < 2m . 20m 
t) = Fay + (« +h, 2m) 
FA) = a0 +B (mp on Ft + by sin 2 


and 


TN 


z TnaliX(t) cos = 2m at, 





% [Efe sin 2 sat, 


fr(t) 





Fig. 60 
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Letting 
A, 


i f: 
In = 9m Be = by = 


felt) = * 2 ZS (Ax coswxt + By sin ot), (194) 
a 


where 
2 
A, = te fit) 008 wnt dt, 
(195) 
B= fe Alt) sin aigt de. 
Let us now set 
2 2 Qn 
Ao, = aye — ay = EY _ 2m 


in Equation (194), which then becomes 
felt) = 84 2S A, con ant Aer, + 1S B, sin unt Aer 
Wael 7 pel 


But since w, = 0, and hence cos wot = 1, sin wet = 0, we may write 
Ay _ ApAey 





Ag cos cet 
te Otte 





and adding and subtracting 4q/T to the last equation, 
fr) = — AB EX Aco toxt toy + By sin eagt Keo,). (196) 


It seems reasonable to surmise that as T —+ , the above infinite 
sums approach integrals. We shall try to make this limiting process 
appear plausible. The integral 


E Aw) cos wt deo 


can be closely approximated using a finite sum by the very 
definition of the definite integral (see [24], Section 1.9). Suppose 
we divide the a-axis from 0 to T into N pieces, each of length 
2n/T. Clearly, the larger T, the greater N, and the smaller the 
length of each interval. Let 


ee SOAs 
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‘Then, approximately, 
xa 
fi A) coe wot deo = Si (erp) 000 cogt Aurgy 
w= 


where Au = cn4;—Wn = 2a)/T. The larger T, the better the 
approximation. 
With this justification, and noting that 
lim —42 = 0, 


me T 


Equation (196) assumes the form 
ft) = + Jo (Alo) 008 ot + Blo) sin ef] deo, (197) 
and by an even simpler argument 


A(w) = f Ai) cos wt dt, 
(198) 
Biw) = { Ae) sin wt dt. 


Equation (197) is one form of the so-called Fourier integral. ‘The 
functions A(w) and B(w) are called the cosine and sine transforms 
or Fourier transforms or spectrum of the function f(t). Conversely, 
(f(t) is often spoken of as the inverse Fourier transform of A(w), B(w). 

Note that while in a Fourier series we get contributions only 
from integral values n of wT/2n, in a Fourier integral, we get 
contributions from every real value of w. 

‘We leave as an exercise for the reader the proof of the following 
‘statement analogous to those for Fourier series : If f(t) is an even 
function, B(w) = 0, and if f(t) is an odd function, A(w) = 0. 


60. Continuation of the solution of Section 58 


We shall use the results of Section 59 to continue the solution 
of the heat flow problem started in Section 58. Returning to 
Equation (193), and letting K = 1/m, u(x,t) becomes 


eet) = 1 Jo ABP cos Be + BBI*M* sin Bs] df, (199) 
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so that the initial condition 1(x,0) = f(x) takes the form 
w.0) = fle) = 1 J [A(B) cos Bx + BIB) sin Bx) dB. 


This equation shows that the unknown functions A(B) and B(8) 
are simply the cosine and sine transforms of f(x), and hence, 
by Equation (198), are given by 


AG) = [2 A) cos pe as, 


B®) = we sin px dx. 


By actually substituting the cosine and sine transforms given above 
into Equation (199), 


1 


wen) — be | [fe A&) £08 BE «i et cos fx 


of [fn JG) sin BE «] crap F 


20 Ds (cos BL cos Bx + sin BE sin Bae“ a 
on changing the order of integration. Recalling the trigonometric 
identity 

cos BE cos Bx + sin Bf sin Be = cos A(f — x), 
u(s,t) assumes the form 


tat) = Lf 10 [JF eee coo ak 2) a8] af. (20) 


This is the solution of our problem whatever f(x). However, 
the second integral 


Tm Ae eet" cos BIL — x) dp (201) 


is a definite integral involving only known functions, and if it 
can be evaluated, Equation (200) can be put into a simpler form. 
We shall compute this integral in the following section, 
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61. A digression on integration 


We shall evaluate certain definite integrals that will be useful 
in the completion of the solution of the previous section. Consider 
the integral 


i= feta 
and its square 


P= (Sc eat as) (I° eo dy), 


which may be written as the double integral 
pea is Se enolate dx dy, (202) 


If we introduce in the above integral a change of variable to 
polar coordinates, 
x= 1co8, y=rsind 


so that as x and y vary from 0 to ~ (the first quadrant), r varies 
from 0 to ©, 6 varies from 0) to 7/2 and dx dy =r dr dé, Equation 
(202) assumes the form 


pea tc rr) ie re dr, 


The value of this integral is 





and hence 


j= f ett de 





‘We can now proceed with the evaluation of the integral 
Kab) = is "=" cos bx dx, 


which is a function of the two parameters a and b. If we 
differentiate I with respect to b and integrate by parts, 
a 


F = — Jo see sin bade = sy (> sin be dle“) 
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= pellsin beje#2] r — abe Je bere coe be de 


b 
= 0-301 
Hence 
a by 
&~ ia" 
Integration of this ordinary differential equation yields 
Heb) = Acre, 
where A is a constant of integration. A is evaluated by noting that 
1(@0) = A= £ ot! de = J = ae 
Hence 
Kap) = YE evn, 
or 
Je eee 0s bo de = Yeon, (203) 


62. Completion of the solution of Section 58 

We are now in a position to complete the solution of the problem 
started at the beginning of this chapter. Returning to Equation 
(201), let 





B=x, ct=a, ¢ 5, 


in Equation (203), which becomes 





fo 20 com BC — 2) 4B = Meee (-“*)- 


Substitution of this quantity in Equation (200) gives finally 
wet) = sbeefo Aetna, 


which is the solution of the problem for any f(x). 
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63. A numerical example 


If the function f(x) of Section 58 is defined in the following 
fashion : 
fle) = 0 #<0, 
As) = T (aconstant)  >0, 


the solution of the last section becomes 
u(x,t) = aor [e etn dt. 
If we introduce the change of variable 
=((—a)eyi sothat df = 2X Vids, 


and when {= 0, 2 = —2/2cy/i, and when [= ©, z= © then 
u(x,t) becomes 


me) = Fe Pa 


(a= Take = 


(= {0 saeetian neni) 


Now 


and u(x,t) assumes the form 


wet) = Je fp ett ds + Je et ae. 
Since, as was shown in Section 61, 
gies F 
wen = Dy e Ee, de, 


The function a 
fh 
oe a 


is known as the probability integral or error function and is written 


Google 


152 ‘THE FOURIER INTEGRAL Art. 64 


‘The properties of this function are well known and its values 
are tabulated, for example, in [22] and [10]. Hence we may write 
our solution in the compact form 


wet) = 2 (1 + ety 





64. Vibrations of an infinite string 
As a second example of the use of the Fourier transform we 


shall determine the displacements u(x,t) of an infinite string. The 
function u(s,t) satisfies the partial differential equation (Section 2), 


Pu 1 Oe 
Be hoe ee 
and we shall assume the boundary conditions 
é 
2 ue0) = 0 
me0) = fe) -e << 208) 


Proceeding as in Section 42, we try a solution of the form X(x) T(t) 
and find, by using the first boundary condition that u,, 


uglst) = Xx) Tt) = (AWB) cos Be + BCA) sin Bs} cos aft, 


is a solution of Equation (205) and satisfies the first boundary 
condition of Equations (206) for any real value of B. Hence 


w(sst) = 1 J? (AG) cos fe + BG) sin Be] cos aft dB 207) 


also has this property. If u(x,t) as given above is to equal f(x) 
att=0, 
u(e,0) = fz) = a JF (AG) cos Bs + BOB) sin Bx} dB, (208) 


and A(f) and B(f) must be the cosine and sine transforms of 
f(z). By Equations (198), therefore, 


A() = [7 Ai) cos BS a, 
C (209 
Bp) = fA) sin pe av. 
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‘Substituting these values for A(8) and B(8) in Equation (207) 
and interchanging the order of integration, 

ats) = 1 f° AO f° (con con Be + sin BU sin Br) 


cos oft dB dt. (210) 
By use of the trigonometric identities 


cos BE cos Ax + sin BE sin Bx = cos A(t — x) 


and 
con (E —) con ft = [os {l= + af) + cosA(C —2 — at), 
Equation (210) becomes 
ate!) = ie fof; Je ME — 2 + a) dp at 
ZL AO fo Somat — 2 —atyap.dt, 1) 


which is the desired solution. 

We can also determine a form of the solution analogous to 
that obtained in Section 43, Substituting (8) and B(B) as given 
by Equations (209) into Equation (208), we obtain 


Fe) = 2 §° AO J? (cont cos Bx + sin BL sin Bx) dB at 
= LIL AO JF coe ae — 2) ap at. 


If we replace x by x + at in the above equation, 
fle at) = 1 f* 4G) [? coo AE — F a) df at, 
and comparing this with Equation (211) we see that 
tei fe + at) ae fle — at) 


65. The Fourier transform 


In this section we shall develop a form of the Fourier integral 
based on the complex form of the Fourier series. ‘This form is 
‘of great importance and has a wide range of applicability even 
n 
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outside the field of partial differential equations.* Problems 
solvable by means of the sine and cosine transforms can, of course, 
be solved by the complex form of the Fourier transform; but 
this form allows the use of a method totally different from 
separation of variables for the solution of partial differential 
equations. 

The complex form of the Fourier series of a function fr(t) 
(Figure 60) is given (Equation 116) in the notation of this chapter, 
by 

Fijon) — [i fale tmn' dt, 
(212) 


St) = eS Fijag) ett 


‘The discrete set of numbers F(jw,) is often called the (discrete) 
spectrum of the periodic function f;(t). If we introduce the 
notation 
2am 2n 
Wy = A Ay = Wnt — Oe = 


as in Section 59, the second of Equations (212) becomes 


ft) = yD Alie,) ee Bary 


and by the same argument used in Section 59, 28 T+ «, Equations 
(212) become 
Foja) = fe) etm dt, 
13) 
$0) = of, Ria) om de. 


The function F{jw) is called the spectrum or Fourier transform 


‘* For example, linear systems, network theory, amplitude modulation, in- 
formation theory. In these cases t in usually interpreted aa “‘ time” and w a2 


“angular frequency.” 
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of f(t), and f(t) is called the inverse Fourier transform of Fljw). 
F(jo) and f(t) are also called Fourier transform pairs.* 

While mathematical rigor has not been stressed in this book,** 
it is essential to note that if F(jw) as given by the first of 
Equations (213) is to be meaningful, the integral must exist, and 
for this to be true, 

lim f(t) = 0. 
iP. 
6. Examples of Fourier transforms 

In this section we shall calculate the transforms of a few functions 
of practical interest. These results, as well as others, are listed 
in Table I. For an extensive table of Fourier transforms, see [4]. 


fn) 


a Tt tro 





Fig. 61 


Consider first a “train of pulses,” that is, the function given 
in Figure 61 which is periodic of period T, and compute its complex 
Fourier series by means of Equations (212) : 


F(ja,) = [opt tet at = h [tere at 


a (os — one) = 2 it one, 
Jun Ws 

* If fit) is an even function, then F(jw) is a real,even function of w. For replacing 
4 by —t in Equation (213), 
jo) = > Aer ay = E fdede dt = F-je) = FGe), 
and FU—a) = je) = PGs) = Fue). 

** A more detailed discussion of questions involving mathematical rigor appears: 
in Section 70. 





Google KiiERE OO RCSAN 


156 ‘THE FOURIER INTEGRAL Art. 66, 


We can actually plot F(jw,) (see Figure 62), since this function 
is real. The ordinates of Figure 62 constitute the discrete spectrum 
of the function f(t}. We also note that the envelope of these 
discrete ordinates is the continuous function 2h (sin wa)/w. 


Fim) 





Fig. 62 

Suppose now that instead of having a train of pulses, as in 
Figure 61, we have a single pulse as in Figure 63. (This pulse 
may be thought of physically as a single dot in a Morse code 
telegraphy system or as a distributed load on an infinite beam 
resting on an elastic foundation.) This function is not periodic, 
but if we evaluate its Fourier transform F(jw) by means of 
Equations (213), 


Rij) = f°, Aovriet at = h ft tet dt = 2h 
mo 


sin wa 
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wwe see that the Fourier transform of the single pulse is the envelope 
of the lines of the discrete Fourier transform of the train of pulses. 
‘The function F(jw) is plotted in Figure 64. 


F (ju) 2h S82. 





Fig. 64 


If we consider the limit of the pulse of Figure 63 obtained 
by letting A» © and a@—+0 in such a manner that the area of 
the pulse, 2ak, remains constant, say equal to K, we obtain the 
so called Dirac deita function, &(t). (In the terminology of the 
civil engineer &(x) represents a concentrated load, and in the 
terminology of the electrical engineer 8(t) represents a unit impulse 
function.) ‘The Fourier transform A(ju) of 8(t) can be obtained 
as the limit of the Fourier transform of the finite pulse, 


Agia) = gg AiO 8 Sit 


K. 








Hence the Fourier transform of a unit impulse function 8(¢) 
(K=1)is 1. 


As the next example we shall calculate the transform of the 
so-called “‘ unit step function” (Figure 65) which is defined by 
the equations 


ut) =0 1 <0, 
wf) = 1 t>0. 
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From Equations (213) 
Uj) = Flee) = [P eset de = 


u(t) 


cme, 
Jwo 


Fig. 65 
But 
lim, ef" = lim (cos wt — j sin wt) 

does not exist. We surmount this difficulty by the following 
device. Consider the function et with «>0. The function 
u(t) e-% (Figure 66) is such that 

u(t) = 0 8<0, 

ut) =e 2 > 0, 





and as a approaches zero, ¢~* u(t) approaches 11). ‘The transform 
of the function e-** u(t) is easy to calculate by ee (213): 


Flues] = [J meta = 1. 


If we define the transform of u(t) as the limit of the transform 
‘of e** u(t) as a approaches zero, we find 
1 1 


Ufa) = Jim, Fes] = im, 55 = jo: 


+ We shall use the notation J to indicate the Fourier transform of a function 
and write F(jw) = Ff]. Correspondingly, we write f(t) — J! [Fijw)) and 
read : f(?) is the inverse Fourier transform of F(ja).” 
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‘Transforms of other elementary functions may be found by 
using the process above, that is, by introducing a factor e-*!(a > 0), 
computing the transform, and then letting «—> 0. 


u(t) eat 





Fig. 66 


67. Properties of Fourier transforms 

An important property of the Fourier transform is the so-called 
“‘ translation property ” which enables us to calculate the transform 
of ft —7) (where 7 is a constant) in terms of the transform of 
fit). By definition 

Fut —2)) = {s fit — re dt 

where r is a constant. In the above integral introduce the change 
of variable: t= €+ 7. Then 


Fue — on) = f% AG) etetern de = eter [* seyetot ae 
= eRe), 
where F( jw) is the Fourier transform of f(t). Hence 
FU = 2) = © FLA). 
If we interpret ¢ as “ time” and 7 as a “ time delay,” we see 
that a delay in time is equivalent to multiplication of the transform 
by an exponential. 
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‘Taste I 














wo Fu! 
1. f@) | Fw) 
2 fe—r) | FUa)e to" 
3. jwF Go) 
4. Ft) Gay"FGw) 
5. Finite pulse of Fig. 63 ah “= 
6. 3) 1 
7. &¢—7) eter 
1 
8. u(t) s 
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10. u(te*, a>0 shi 
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13. u(t) Gor =-Z 
14. u(t) sin at gare" 

jw 
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Another important transform formula gives the expression for 
the transform of the derivative of a function f(t). Write 


Furol = f2, HO emt at. 
If we integrate the above integral by parts, 


+o 
FIO) =f + Ga) J Aiea 
and assume, as was remarked in Section 65, that 
Hf = 0 
we find that 
FU) = joFGe), 
where F(jw) is the transform of j(t). We see, then, that 
differentiation with respect to t is equivalent to multiplication 


of the transform by jw. If f'(t) also vanishes at --~, we can compute 
the Fourier transform of f(t) in the same fashion, 

FUP) = (Giw)*F(jo) = —w*F(juo), 
and similarly for higher derivatives. 

In all the examples given above, transforms were computed 
starting with a function of time (). It is clear that we could instead 
have started with the transform (a function of w) and obtained 
the time function by use of the inverse Fourier transform formulas 
Equation (213). However, in general, this is not an easy process, 
and frequently necessitates the evaluation of integrals by contour 
integration in the complex plane—a topic beyond the scope of 
this book. 


68. Transform of the heat equation 


In this section we shall apply the Fourier transform to the 
solution of the heat flow problem solved in Section 62. We recall 
that the boundary value problem was given by 

oo ao 

oO oe 


(20) = fz) 0 <<a. ae 
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We shall take the transform of the above differential system, 
solve for the transform of the unknown function v, and apply 
the inverse Fourier transform of Equation (213). ‘The transform 
of o(x,t) will be taken with respect to x, letting ¢ be a parameter. 
Then, by definition, indicating by F, the Fourier transform taken 
with respect to x, oe 


Floto) = VUE!) = [2 oes tae 


(where ¢ is the transformed variable corresponding to x). From 


‘Table I, 
Fel Zs e9] = Gerruen = -erUE0, is) 
while 


Fe Le oe] = ia 2 o(e,te te de 


-2 [ Ji, neve +] - £10. 216) 


Operating by the Fourier transform on both sides of Equation 
(214) and using Equations (215) and (216), we obtain the trans- 
formed equation 


4 MGE») = 88 VUE. 17 
‘The transformed boundary condition is 
VUEO) = Flks0] = FAP = FG. ale) 
Equations (217) and (218) form an ordinary differential system 
on V. The general solution of Equation (217) is 
Ves) = AGH te, 
where A, the “‘ constant of integration ” is a function of ¢. From 
the boundary condition of Equation (218), 


FUE) = VEO) = Alife? = ALE), 
and hence 


VE) = FGF (219) 
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‘The desired solution v(x,t) is the inverse Fourier transform (with 
respect to *) of V(j,t), and by Equations (213), 


oe) = 2 IP, Mipestiee at, (20) 


which is the required result. 
If, for example, as in Section 63 we let f(x) = Tu(2), then 


PUD) = Fats) = F 
and 
Hes) = ao fh ge eetete a. iy 


‘The inverse transform of (1/jge~**" can be found in [4]. In 
fact, Formula 727 of [4] yields (in our notation) 


Fg a] = 2 (stave -})- 


1 (1 ne) =i xn) 
Filge =a (I +e ZF): 
Hence we conclude as in Section 63 that 


o(z,t) = F(! +erf an): (222) 





or 


L 
69. Transform of the wave equation 


‘As a second example of this technique we shall solve the 
illustrative problem of Section 64. The boundary value problem is 


Fo a 
of OF 
(6,0) = fle) 2 <3 <2, (23) 
2 x0)—0. 


‘Transforming this differential system with respect to x we have, 
letting V(jé,t) = F.{o(2,t)], 


Fe[ Zoe] = vErvuE, 224) 
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Fe[eeemo] = fo Zenner ae 

= a Viet), (225) 

Felels0)] = VEO) = FUE) = FUO, 226) 

F.[Z ote] = ZFdo%0) = 2 rue. — aan 


From Equations (224), (225), (226), (227) we obtain the ordinary 
differential system 


Fe YU69 = oUEVUE, @2sy 
VE) = Fue, ei 
Bes 
2 1UE0 = 0, 


on the transformed function V. The general solution of Equation 


is 
Ce) VGEt) = AGE + Bis H*. 
The constants of integration A and B (with respect to ¢) can 
be evaluated from Equations (229) : 


FAs) = VUE) = AGE) + BU, 
0 = 4 vGE0) = jag Ait) — jab BLE. 
We conclude from these two algebraic equations that 
AGO = BUD = L7G 
and hence that 
VUE) = J PUB (ot + eH, 230) 


The desired function o(s,t) is the inverse Fourier transform of 
V(jé,t), namely, 
otet) = Lf viene ae 
2n jo 


= LP rune tenmeras. aay 
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This equation gives the solution of the problem for any f(+). 
By means of the translation property (see Table I), 


Fe FUE eer] = fle + at), 
FAG] = fle — at). 


(et) = SREB + FORGO He} 


= Fite + at) + fle — an, 232) 
which is the familiar equation for the traveling wave. 


hence 


70. Mathematical conditions on the Fourier integral 
‘The Fourier transform is defined only for functions which are 
absolutely integrable, that is for functions f(x) such that 


[Ue ld <e. (233) 
The Fourier transform F(w) of such an f(x) is defined as 
Fla) = JP Aeyrioe de, 4) 


and can be shown to be a continuous function of w. 
On the other hand, Equation (234) does not imply that f(x) 
is the inverse Fourier transform of F(w), that is, that 


1 fe 
fe) = 3 SE Fw) deo, (235) 
unless certain additional conditions are satisfied. A sufficient 
condition for Equation (235) to be valid at a point x = xp is 
that the right- and left-hand derivatives of f(x) exist at x = x,* 
in which case 


A feet) + fa = af, Meyennde. 236) 


* By definition, the left-hand derivative of f(s) at the point x = 2 is 
Pee) = tims), 
8 
and the right-hand daivative of fl) atthe point = ayia 
Pvt) = lin fo) 
(eee Section 26), Bh 
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If, moreover, (x) is continuous at x = so then flx)+) = fl%—), 
and Equation (236) becomes identical with Equation (235). 
Equation (235) is valid for all + (and not just at a point as 
described above) under various sets of hypotheses. Perhaps the 
simplest of these is given by the following theorem. 
Turormm. If f(z) and F(w) are absolutely integrable, /(2) is 
continuous, and 
Flu) = [%, fleet de, 
then 
fle) = 2 fn Flue! dw 
holds for all x. 
For a further discussion of the subject including a wide variety 
of theorems and proofs, see [2]. 


PROBLEMS 
1. Show that if f(t) is an even function (see Section 59), B(w) = 0, 
and if f(t) is an odd function, A(w) = 0. 


2. Find the cosine and sine Fourier transforms A(w) and B(w) of 
the function f(é) given in Figure 67. 


(6) 





Fig. 67 
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3. Find A(w) and B(a:) for the function of Figure 68. 
He) 





Fig. 68 


4. Find A(w) and B(w) for the function of Figure 69. 
1 





Fig. 69 


5. Find A(w) and B(w) for the function of Figure 70. 
ue) 
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6. Find Aw) and B(w) for the function f(é) : 


MH =1-F -1<€ <1, 
1 =0 otherwise. 

7. Plot the functions A(w) and B(w) for Problems 2, 3, and 4 so 
that they can be used for any values of the parameters = and m, that 
is, use nondimensional parameters. 

8 Find the Fourier transform FYjw) of each of the functions of 
Problems 2, 3, 4, 5, and 6, and plot the magnitude and phase of F(jw) 
for each. 

9. The lateral surface of a semi-infinite bar (one-dimensional case) 
is insulated, and the end x = 0 is kept at zero temperature. Find the 
temperature distribution u(x,t) if the initial temperature u(z,0) is a constant 
T. Plot this temperature distribution between x = 0 and x = 10 at 
t = 4 sec and t = 10 sec, Assume A/c} = 0.25, where k = thermal 
conductivity, ¢ = specific heat, 5 = density. 

10, Why is it necessary to assume u(x,0) = — T for negative values 
of x in Problem 9? 

Il. Find the temperature u(x,1) in an infinite bar in which the initial 
temperature distribution is 

ux,0) = Tel, a > 0. 
Asoume the lateral surface of the bar is insulated. 


12. The lateral surface of a semi-infinite bar is insulated and at its 

end x = 0 the temperature gradient 

a 

Hut) 
is zero. Find the temperature distribution u(x,t) if u(x,0) = f(x), 
O<z< 

13. Consider a thin slab of infinite extent whose lateral surfaces are 
insulated. Let (x, y) be the initial temperature distribution in the plate. 
Find the temperature u(x, y, t). 

14, Find the temperature distribution in a semi-infinite slab 
(Figure 71) if the initial temperature distribution is u(x, y,0) = f(z, 9), 
O<e< a, —& <y< ©, and edge of the plate, x = 0 is kept at 
zero temperature. Assume that the lateral surfaces are insulated. 
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/ 


Fig. 71 


15. Solve Problem 14 if, instead of having 1(0, y, t) = 0, this condition 

is replaced by the adiabatic condition 
a 
Fetus) = 0. 

16. Consider a‘ quarter-infinite ” plate (Figure 72) insulated on both 
lateral surfaces. Let the initial temperature distribution be u(x, y, 0) 
$29), J<*< =, O<y< om. Find u(s,,,#) if the edges » — 0 
and yJ= 0 are kept at zero temperature. 

y 


° 


1 
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17. Solve Problem 16 if the conditions (0, y,#) = 0 = u(x, 0, ) 
are replaced by adiabactic conditions. 

18, Solve Problem 16 if the boundary conditions are 

1. xy,0) = Hey), 
2. (0,y,t) = 0, 
3 MeO) =10. 

19. A. three-dimensional infinite solid has an initial temperature 
distribution u(x, y, #,0) = (x, 9,4), — © <*#<—, —@<y<™, 
— = <2< ©, Find u(s,y, 9,2). 

20. A semi-infinite string, fixed at the origin (Figure 73) has an initial 


displacement and velocity given by fix) and g(x), respectively. Find 
the transverse displacement w(x, t) of the string. 


y 


io 
Fig. 73 
21, Solve Problem 20 if the boundary conditions are u(x, 0) — fle) 
and a 
#0) = 0. 
22. Solve Problem 20 with the boundary conditions 
1. 40,1) = 0, 
2. Wx,0) = 0, 





3. 2 u(x0) = ale). 


23. Consider a membrane, infinite in extent, whose initial dis- 
placement is given by u(z, 7,0) = (x,y) and whose initial velocity is 
zero. Find the transverse displacements u(x, y, £) of this membrane. 


24. Solve Problem 23 for the boundary conditions 
1. a(xy,0) = 0, 


2) Zu2y0) = xy). 
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25. A semi-infinite membrane (Figure 71) is held fast along the y-axis. 
Determine the displacement u(x, y, t) if the initial velocity is zero and the 
initial displacement is u(s, y,0) = o(s,9), O<* <@, —@<y <a 


26. Solve Problem 25 if the initial displacement is zero and the initial 
velocity a 
ae 70) 

is We, 9). 

27. Figure 72 represents a “ quarter-infinite” membrane held fixed 
along the x- and y-axes. Find the displacement u(z, y, t) if u(x, y,0) = 
3,9) and a 

a 4=y0) = 0. 

28. Find the displacement uz, y, t) of a quarter-infinite vibrating 

membrane if ¢(x, y) = 0 and 
a 
290) = Kay) 


29. A semi-infinite shaft, fixed at the origin is vibrating torsionally. 
+Find the equation of motion for O(+,) if the initial angle of twist 
(2, 0) = F(x) and initial velocity is zero (see Problem 11 of Chapter III). 


30. An infinite shaft vibrating torsionally has an initial velocity 
a 
Fae0) = Ge) 
and zero initial angle of twist. Find the angle of twist (2, #) at any time f. 







31. A high-frequency transmission line, infinite in extent, has an 
initial voltage distribution e(, 0) = E(x) and an initial current distribution 
i{z, 0) = I(x). Find the currents and voltages in the line at any point = 
and at any time ¢. 

32. The high-frequency line of Problem 31 is semi-infinite in extent, 
with zero voltage at the sending end (origin). Find es, t) and i(s, 1) 
in a line subject to the same boundary conditions as in Problem 31. 

33. A submarine cable, semi-infinite in extent has zero voltage at 
the sending end and an initial voltage distribution ¢(x,0) = B(s), 
0 <s<-«. Find the voltage e(x, t) at any time f. 

34. A submarine cable, infinite in extent, has an initial current 
distribution ix, 0) = I(x). Find i(s, 1). 
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35. The function 
Si(s) = [84 ae 


is known as the “ sine-integral manbee Show that lim Si(s) = 22. 
[Hint : Consider the integral 
f ao a 5a 
where o is a parameter.) 
36. Prove : 





1 
; wala = log (1 +0), 2 = 0. 
37. Show that 


[leon ae 





38. Prove that 
* loge if #21 
fF 0g (1 —2¢ cos x + 1%) dx = fF eed 

39. Establish the formula 


f e* sin x 
0 


7 
ya = 5. 


40. Prove that 
FIRE —2)] = eter 
by considering 8(t —r) as the limiting case of the finite pulse of Figure 74. 
f(t) 





Fig. 74 
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41. Prove that 
Fimo) = Gap 
where u(t) is the unit step function. 
42. Find the Fourier transforms of the functions 
u(t) sin at, — u(t) cos bt, 
where (¢) is the unit step function. 
43. Show that the Fourier transform of e~*!*, o > 0 is given by 
2a 
efor 
44. Calculate the discrete or continuous spectrum of the following 
functions and plot the results. 
(a) The periodic function sin t. 
(b) The periodic function f(t) of period 4Nn defined in the interval 
(0, 4Nm) as 
Slt) = sint 0 <t <2Nn, 
ft) =0 2Na <t < 4Nn. 
(c) The nonperiodic function g(t) defined as follows : 
a(t) = sint 0<i<e, 
at) = 0 —m <ct1<0 
(a) The nonperiodic function A(t) defined as follows : 
Ki) = sint = —2Nw <t<o, 
Ht) = 0 2 <1 <—2Nn. 
In parts (b) and (d) assume that N is a large positive integer. Which 
results can be considered as limiting cases of the others? Discuss these 
results. 


45. By a direct application of the inverse Fourier transform formula 
[Equations (213)] calculate the inverse Fourier uansform of F(jw) = 


(=) 
Solve the following problems by the methods of Sections 68, 69; that 


is, calculate the Fourier transform of the unknown function and solve 
the resulting transformed equation for the unknown Fourier transform. 
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. Problem 11. 
. Problem 12. 


Problem 20. 


). Problem 21. 


Problem 22. 


. Problem 29. 
}. Problem 30. 
}. Problem 31. 


Problem 32. 
Problem 33. 
Problem 34. 
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Chapter V 


LEGENDRE, BESSEL, 
AND MATHIEU FUNCTIONS 


71. Some important ordinary differential equations 

In previous examples when we used the method of separation 
of variables to solve a partial differential equation, the resulting 
ordinary differential equations were generally linear with constant 
coefficients. In this chapter we shall consider certain partial 
differential equations of engineering importance which give rise 
to variable coefficient ordinary differential equations of great 
interest, namely, Legendre’s equations, Bessel’s equation, and 
Mathieu's equation. 

‘We shall discuss first certain essential preliminaries such as 
orthogonal functions and the gamma function, and then solve 
Legendre’s, Bessel’s, and Mathieu’s equations and study some 
‘of their more elementary properties. At this stage of our 
development we shall be in a position to solve completely the 
original boundary value problems which led to these ordinary 
differential equations. 


‘72. Legendre’s equation 

Consider a sphere of radius A with a certain potential 
(gravitational, electric, or magnetic) distribution on the surface; 
we wish to find the potential distribution at interior points of 
the sphere. The potential u(x,y,2) satisfies the three-dimensional 
Laplacian V4 = 0, 


where, in Cartesian coordinates, 
a lala 
og OO 
Vegtgt ar 

175 
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However, since our problem concerns a sphere it is convenient 
to transform V® into spherical coordinates, just as in Section 50 
where we changed the two-dimensional Laplacian V? into polar 
coordinates because our problem involved a semicircular slab. The 
spherical coordinate system we shall usc is illustrated in Figure 75 





Fig. 75 


where in terms of r, 6, and ¢, the Cartesian coordinates x, y, 
and 2 are given by 

x= rosésin®, y=rsingsind, x =r cod. 
From these equations we may calculate 


ie ae cle 


in terms of r, 0, and ¢, and write 

1a eu 1 a Ou’ 1 au 
7a (-%) + ans oo (80%) + into OF 
‘The details of the partial differentiations follow the same pattern 
as in Section 50 and will be left as an exercise for the reader. 


Vu = 
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In particular, if the potential is symmetric about the x-axis, 


ou 
Bre 


and the above equation assumes the form 


‘ a ou 2 (.. 2 ae 
sino 3 (n 2) + & (vin 055) 
‘The boundary condition, giving the initial potential on the surface 
of the sphere can now be written 
48) = f@), (238) 
where f(9) is the given potential distribution. 

Equations (237) and (238) form a differential system which we 
shall attempt to solve by the method of separation of variables. 
Substituting a solution of the form 

(7,8) = R()8(0) 
into Equation (237), we obtain 





0. (237) 


0 sino 4 (ary + REO sin) =o, 
or dividing by RO sin 0, 
Td 5 
Gen 8 ao sin 6) = -feer,. 
By the usual argument, the two sides of this equation must be 
equal to a constant, say A, so that 


(8' sin 8) — 20 = 0, (2392) 


als 


1 
sind 


Ma 


gy OR) +R — 0. (2398) 


The second of these equations becomes, on performing the 
indicated differentiations, 


AR" + WR + R= 0 
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which we recognize as Cauchy’s equation (see Section 50). If 
we introduce the change of variable 


x = 1088, 
a iad _ _ gd 
a6 ~ dds O° 


in Equation (2392), this equation assumes the form 


t 


+4 vint 2) —20 =0, 


or 
a 4 - 
Z[« -9 2] —r0 =0. 





‘This is one of the fundamental equations of mathematical physics 
and is called Legendre’s equation. If, as is customary, we let 

A= =nln + 1), 
Legendre’s equation assumes the classical form 


Z[d-9 2] + e+ ne =0. (240) 


‘The properties of this equation are studied in Sections 77-79. 


73. Besse!’s equation 
Suppose we have a thin circular disk insulated on its lateral 
surfaces, whose temperature depends only on the distance from 
the center (and not on the direction). We wish to determine the 
transient temperature distribution in the disk for a given set of 
boundary conditions. The partial differential equation governing 
this phenomenon is 
ue Lae 1 ee 
a 


oa +e ba mo mH 
If we call a the radius of the disk, typical boundary conditions 
that could be associated with Equation (241) are 
u(a,t) = 0, ufr,0) = fir). 


* This equation was derived in Problem 18 of Chapter I directly from physical 
considerations. 
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These indicate that the temperature on the boundary of the disk 
is permanently zero, while the initial temperature distribution is 
a given function of r, say f{r). 
As in previous examples, we shall substitute in Equation (241) 
a solution of the form 
wrt) = Rr) Te) 


and obtain by separating the variables, 


LR 1? 
Perea ae * 





where A is a constant. If A is negative, A——#t, the solution 

of the equation in 7(t) will contain a negative exponential, as 

one expects from physical considerations, and the above equation 
gives rise to the two equatiors, - 

T+ APT =0, (2420) 

r+ te + pR=0. (2426) 


"The second of these equations is another fundamental equation 
of mathematical physics and is known as Bessel’s equation. If 
we let s = Br, Equation (242b) becomes 


0, 


| 
* 
| 
| 
+ 
~ 
i} 


or 
@R aR 
ets tar = 0. (243) 


‘This is the classical form of the so-called Bessel’s equation of the 
zero order. The more general equation 


@R dR 
Sets + or = 0, (244) 


which arises in other boundary value problems, is known as Bessel’s 
equation of the nth order. Clearly Equation (243) is the particular 
case of Equation (244) corresponding to n = 0. 

We shall solve this equation and make a study of its properties 
in later sections (Sections 81-83). 
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74, Mathieu's equation 
Let us consider an elliptic membrane, fixed on the boundary 

and subjected to an initial displacement and an initial velocity. 
We wish to determine the transverse displacements of the 
membrane. The partial differential equation (in Cartesian 
coordinates) governing this phenomenon is (see Section 4), 

CE ater a) 

at T Bt ~ at at 
It will be convenient at this point to climinate the time variable. 
Since essentially we shall assume a solution with the separated 
variables, and since the free vibrations of the membrane vary 
sinusoidally with time, we let 

u(x,y,t) = o(x,y) (Acos wt + B sin wt). 


Substituting this equation into Equation (245), we obtain 


(245) 


Om aw 2 
ae t pit Gre =o (246) 


which is an equation in the two variables x and y. 

Because of the shape of the membrane it is convenient to 
transform Equation (246) into elliptic coordinates (&,), by the 
transformation 

* = ycoshfcosm, y = ysinh &siny, (241) 

where y is a constant. The variable 7 varies from 0 to 2m, and 
the variable € from 0 to». 

‘We shall discuss this transformation in some detail. If we 
let € be a constant, Equation (247) defines a locus which is an 
ellipse, since 





7 ; 
yeoke = °° — yainkg = 9% 

and é ‘i 
Geo FF * Gein BF = om 
‘The semimajor axis of this ellipse is y cosh g, and the semiminor 
axis is y sinh (Figure 76). ‘The distance between the foci 


2Viq cosh EF — (sink FF = 2y. 
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Hence we see that the curves £ = constant form a family of 
confocal ellipses. When ¢ 0 we obtain the degenerate ellipse 
which is the straight line between (—y,0) and (7,0) (Figure 76). 
As > ©, the ellipse approaches a circle of infinite radius with 
center at the origin. 





Fig. 76 
If, on the other hand, 7 = constant, Equations (247) yield 
eo Fg 
peng NG = sinh 
pee ee 
(y cosy? (y sin 
‘These curves represent hyperbolas with principal axes coincident 
with the x-axis. The distance between their foci is 
2V(y cos 9) + (y sin) = 2y. 
Hence the hyperbolas are confocal and, moreover, have the same 
foci as the ellipses, £ = constant. 


= cosht £ — sinht é = 1. 
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If » = 0, we find from Equations (247) that 
x= yoohi, y=0, 


and the corresponding hyperbola reduces to the straight line 
extending along the x-axis from x= y to x= +o. Similarly, 
ify =7, 





x=—youhf, y=0, 


and the hyperbola reduces to the straight line extending along 
the x-axis from x = —y tox = —=. Also, if 7 = 7/2, 
2=0, y=yesinh é 

and the corresponding locus is the upper half of the y-axis. For 
1 = 3n/2 the locus is x= 0, y = —y sinh , or that portion of 
the y-axis extending from y = 0 to y = —=. Ina similar fashion 
we may determine the remaining portions of the hyperbols 
sketched in Figure 76. 

We leave as an exercise for the reader the problem of expressing 
Equation (246) in elliptic coordinates. ‘The result is 








SE + Ga + Mlcoeh 2F — om 2mm = 0, (249) 
where oy 
te em) 


Suppose then that our elliptic membrane has semimajor axis « 

and semiminor axis B. Its boundary equation is 

ae 

atg@ol 
(Figure 77). The constant y of Equation (247) may be chosen 
so that the membrane is one of the confocal ellipses of Equation 
(247). For this purpose notice that the distance between the 
foci of the membrane ellipse is 


ave —B, 


y= Vor — Be 
—-p. 251) 
Since a and are given, this equation determines 7. 


and setting 
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In order to solve Equation (249) we substitute in it a solution 


of the form 
w(En) = S(€) HC). 
and separate the variables, obtaining 
E _H 


“a + Mewh 2g = 57 


z + Bt cos 29. 


: 
(@,0) 


Fig. 77 
Calling the constant value of the two sides of this equation, 
we obtain the two ordinary differential equations 

5" —(A—Fooh 298 = 0 

H” + (A—# cos nH = 0. 
Equations (252) represent another fundamental type of ordinary 
differential equations known as Mathieu equations. In later sections 
we shall investigate the properties of their solutions and then 
complete the solution of the boundary value problem for the 
elliptic membrane. 


(252) 


75. Orthogonal functions 
At this point we shall consider certain concepts necessary to 
a more complete understanding and appreciation of the special 
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equations derived in this chapter. We shall revert to Fourier 
series and consider some of their more abstract properties. In 
this connection it is convenient to relabel the functions cos nz, 
n—=0,1,2,.. and sin nz, n= 1,2,... as a single infinitude of 
functions, say, as indicated below : 
02) = 1. 
Ale) 
xz) 
$e) 
dz) 





Fm) 
Panl*) 





sin nx. 


Calling the inner product of two functions f(x) and g(x), the integral 
(a) = [°Ae) ate) de, 


the integrals appearing in the Fourier coefficients of f(z) may 
be indicated in this notation by the symbol 


(iba) = J Me) dala) a. 
‘The Fourier series expansion for an arbitrary function f(x) in 
the interval (—1,7) then takes the form 
1 1S 
Sle) = zAldoWbls) + — 2 (fbadbal®) 
and the orthogonality conditions for the sine and cosine functions 
become 


(Gnd) = 78am %,m integers not both zero, (ofo) = 2m. 


We shall show that Legendre, Bessel, and Mathieu functions 
(solutions of Legendre's, Bessel’s, and Mathieu's equations, 
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respectively) also satisfy orthogonality conditions, and that an 
arbitrary function f(x) can be expanded in an infinite series of 
Legendre, Bessel, or Mathieu functions. 


We shall first consider some elementary properties of orthogonal 
(functions in general. A set of functions 


Gil2) Galt), Pal), 
is said to be orthogonal in the interval [a,b] if 


J'BAedbMe) de = 0, x85, 


J * dda) dx v4 0. 
If, say, 
[FGF a = 0,7 0, 
we can normalize the ¢((x), that is, find a new set of functions, 


$2) = b4x)/V/a which are orthonormal in the interval [2,6], that 
is, such that 


SL PeWde) de = By. 
In the inner product notation, the set of functions {i,(z)} is 
orthonormal if 
Wold = Bu 


Let us attempt to expand the arbitrary function f(x) in terms 
of the orthonormal functions y,(2), that is, let us write 


fe) = x abl), 


where the a, are undetermined constants. By the condition 
(Hoshs) = Bes we have 


(Fh) = (Zagat) = LayGbush) = EanPan = aa 
Hence we may write 


fia) = ¥ ibedile. (233) 
% 
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So far our work has been purely formal, Does the above equation 
actually represent an identity? Not always, since, for example, 
the functions sin x, sin 2x, sin 3x, ... are orthogonal in (—n,7) and 
yet we cannot expand cos x in terms of this orthogonal set. For, 
if we attempt to, we see that 


(cos x, sin Rx) = 0 
for all k, and we have the absurd result that 
cos x = 0. 


Whenever Equation (253) does hold for all values of x, the 
orthonormal set functions {y,(2)} is ssid to be complete over 
the interval [a,b]. Or, more briefly, the {¥(x)} form a complete 
orthonormal set over [a,b]. 

The trigonometric functions sin kx and cos kx, the complex 
exponential functions, the Bessel functions, the Legendre functions, 
and the Mathieu functions are complete orthogonal sets, and 
therein lies their usefulness : we can expand an arbitrary function 
in terms of Bessel, Legendre, or Mathieu functions just as we 
did for sines and cosines. On the other hand, the sine functions 
are not complete in the interval (—z,r), as shown by the above 
example; whereas they are complete over (0,7). 





76. The gamma function 


As we extend our knowledge of mathematics, certain functions 
other than the elementary ones appear in our work. For example, 
we have previously met the error function, erf (Section 63) 
and the sine integral function, Si(2), (Problem 35, Chapter IV), 
and we shall soon be compelled to study Legendre, Bessel, and 
Mathieu functions. 

In the definition of Bessel and Legendre functions it is 
convenient to introduce another new function which can be 
considered a generalization of factorial. We recall, that if » is 
a positive integer, 


al = nln — 1) (n — 2) 24 
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while the symbol «! is at present meaningless when « is not a 
positive integer. A function which is defined for all values of 
the variable a and equals «! when a is a positive integer, is the 
gamma function, T(q). It is defined by the equation 


Te) = fp tet ds, > 0. (254) 
It is easy to show that if a is a positive integer, 
T(@) = (@—1). 
Integrating by parts the integral of Equation (254) : 


Tha) = —ertee 2 + [Ee — tatters de 


0 + @=1) fo ate de 


= (@—-IMe—1)," a> h (255) 
Applying this result repeatedly, we see that if « is an integer, 
Te —1) = « —2F(@ —2), 

T@ —2) = @—30a —3), 


TQ) = 20(2), 















TQ) = IT), 
* For any «, 
lim tet = 0. 
For oe 
Be 
- + 
eet tee 5 
Now choote N so large that V > a—|. ‘Then for large =, 
Gee ee 
Ca ee ate 
mt *wFDI M * WFD! 


and since x —a + 1 > 0 for all 2 N, every term in the denominator approaches 
infinity as x approaches infinity, and hence x*~te~* approaches zero. 
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from which 

Ta) = (@ —1) (« —2) (« —3) » 21T(). 
‘The value of I(1) can be found by direct integration : 


T) = i werd = [° erde = 1, 


T@) =(@-D! «a =12 


Equation (255), however, holds for any positive value of the 
argument and shows that if [(a) is known for 0<a<l, the 
gamma function can be determined for any positive value of a. 
For if B is any positive number, we may always write 


Bp=mta, 





and hence 





where m is an integer andOQ<a<l. Ifa=1, then 
16) = 
since f is an integer. If a # 1, then by Equation (255), 
Tp) = 6-1) @ —2)... 6 —mF@), 


which is readily computed when Ta) is known. For the values 
of (a) see [14] and Figure 78. 

So far we have defined T(a) only for a > 0. If « is negative 
we define F(a) by means of Equation (255) as 


FeLe Met 1) 


HN +1)<a<—-N (256) 
for N = 0, 1,2, 


1-32) - AP) _ ae [232] - 63 a2) 5O5 [2 a 


no) = 


For example, 








amenTAeHTO” 


© 2) (2-2)(12) a (12) [ 
0.6891. 
As an immediate corollary of the above definition, 


T@) =(«—IT(e—1)  « #0,—1,-2... 
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When a = 0 (or a negative integer), the gamma function becomes 
infinite (see Figure 78). This can be shown analytically by writing 
Equation (256) with « = « > 0, 
— Me+)) 
no = Met 


and then taking the limit as e+ 0. Ase 0, [e+ 1) + (1) =1 
and 

T(0)— lim T%) tim, e+) me 
For negative integral values of N, Equation (256), applied 
repeatedly, gives 


TiN) = =, TO) 


1 
CNYCNF ICN +2). (1) 


W 


qy 
GF ro) = +. 


71. Solution of Legendre’s equation 
We have called Legendre's equation in Section 72 the equation 
[Equation (240)] : 

(1 —a8)y" — Day! + nie + Dy = 0. (257) 

‘This ordinary differential equation may be solved by the method 

of Frobenius (see [25], Chapter VII). Substituting y = =” into 
Equation (257), we obtain 

{on — 1-4 — (om — Le — Qe +n + Ie 
= [m(m — 1)]* + [-mt — m + n(n + 1 (258) 


and 
m(m —1) =O implies m = 0,1, 
(259) 
m+ m—n(n+1) = 0 implies m = n,—-(n + 1). 
‘The smaller of the exponents of x in Equation (258) is m—2; 
hence by means of the indicial equation m(m— 1) = 0 we may 
‘obtain solutions of Legendre’s equation in ascending powers of x. 


Google 


Art.77 LEGENDRE, BESSEL, MATHIEU FUNCTIONS 191 


Since one of the roots of m(m — 1) = 0 is zero, we may assume 
a solution of the form 


He) = eo + eget + eget +. = Set (260) 
a 
If we substitute this in Equation (257), there results 


S anak — ews — 3 22k — Nenstt — 2 3 hens 
% % & 


+ n(n +1) Y cust! = 
& 
‘or collecting powers of x: 


lee + 292k + Dees — 2424 — Nien 
— Ake + (m+ Ip) — 0. (261) 
If Equation (260) is to satisfy Legendre’s equation, the coefficients 
of all powers of s in Equation (261) must vanish, and hence 
ey + n(n + leq = 0 coefficient of 2°, 
4+ 3g —2- Mey —2-2ey + n(n + Ney = coefficient of x*, 
6+ Seq —4+3cq — 2 4ey + nln + Ie coefficient of x4, 











From the equation for the coefficient of x0 we obtain 
ey a mee (262) 


“The second equation gives 
- Bonet O,, — CH DOF), 


and proceeding in this fashion we obtain all the cy, and hence 
the integral of Legendre’s equation : 


We) = ceysl#), 
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n(n +1).4 (8-2 (D+NO+3), 
Sy a 
(2 = 4) (m —2 (n+ 1) (H+ 3H 4 9) 
Ga =O EDO FDO +I 04. 063) 


Since the difference of the roots of the equation m(m — 1) =O 
is not a multiple of the difference of the exponents of x of Equation 
(258), we obtain another linearly independent solution 


. x const, 
é 


corresponding to the root m— 1. [If the second root had been 
a multiple of 2, this would be the same solution as Equation (260) 
since it would then involve only even powers.) Substituting the 
above series into Equation (257) and equating coefficient of powers 
of x to zero we obtain 


Hee) = eye + cae + og + 





WE) = arl=), 


where 


piejax =D EH + = )O-DOFACTD 


— (2=5) (2 —3) (21) (n+ 2) (9 +4) (2+ 6), 
1 


Hee (264) 
Using the test-ratio test, [24], Section 5.10, we can show that 
the power series y,(z) and y,(x) are convergent for |z| <1. Hence 
if |x| < 1, the general solution of Legendre’s equation is given by 
He) = Ay(*) + Byrle), 
where A and B are arbitrary constants and y,(x) and y,(x) are 
given by Equations (263) and (264), respectively. 
Turning now to the second indicial equation m*-+m—n(n-+1)=0, 
we assume first a solution of the form 


We) = age + ayxr® + ager + ... 
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in descending powers of x corresponding to the root m = 
the root m = —n —I leads to a solution of the form 


Ha) = bpd + byes $ bg + w. 


Equating the coefficients of all the powers of x to zero, we 
determine as before two linearly independent solutions of 
Legendre's equation : 
nn —1) a, Mn —1)(n—2) (0 
ne) = — 90_ iy" + “Sin — 1) Qn —3) 
n(n — 1) (n — 2) (wn — 3) (x —4) (0 = 5). 
— Hilde =) Qa —3)@e—3) 
bow (268) 
4 OHNO ‘ 2) 
20n+3 
4 et peed 3) +4), 
BHDIGn + 3) Qn + 5) 


4 O+NO+2) +3) O44) (n+ 5) (0+ 6). 
253I(2n+ 3) (n+ 3) (n+ 7) 
Ay (266) 


‘These two power series may be shown to be convergent for |a| > 1. 
Hence 








Bucs 


yat(s) =m ) gna 








He) = Cz) + Dyr*e) (267) 


is the general solution of Legendre’s equation for |x| > 1 where 
C, D are arbitrary constants, and y,*(2), y:*(x) are given by 
Equations (265) and (266), respectively. 

‘A careful analysis of Equations (263), (264), (265), and (266) 
shows that the functions 7, ¥2, 91%, ya" are related to one another 
when 1 is an integer. In particular it may be shown that 


(a) If nis a nonnegative integer, 
ne) = Ain), 
ni") = hayes) 


if mis odd, where k, and ky are nonzero constants. 


if nis even, and 
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(b) If m is a negative integer, 
n*(2) = y(*) 


IaX(®) = beonl=) 
if n is even, where ky and k, are nonzero constants. 

We see therefore that if m is a nonnegative integer, the series 
defined by y,*(2) is finite, since after a certain term each numerator 
in the coefficients of the powers of x will contain a zero. These 
finite solutions are essentially the so-called Legendre polynomials. 
It is customary to define the Legendre polynomials, P,(x) for 
n= 0,1,2,... by the equation 


— a 

PA) = Fone 

"These polynomials are also called sonal harmonics. ‘The functions 
O.(2) defined by 


if n is odd, and 


*()- (268) 


20) = ge ovr 


are sometimes called the Legendre functions of the second kind or 
sonal harmonics of the second kind.* 

We leave as an exercise for the reader the task of verifying 
that the first few Legendre polynomials are given by : 


Pfx) = 1. 
Px) = =. 
3 1 
Pye) = ae — > 
3 (269) 
Pe) = oe — Fe 
35. 15 3 
ge eh ee 


P(e) = 


* If m is a nonnegative integer the general sclution of Legendre’s equations 


for |z| > 1 is given by 
2) = CrP yl) + CQal@s 

where C; and C; are arbitrary constania. If [s| < 1, then the Pq(x) are still solutions 

of Legendre’ equation, but the Q,(+) must be replaced by 9,(+«) or 74(2), depending 

on whether # is odd or even. 
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78. Orthogonality of the Legendre polynomials 

The usefulness of the Legendre polynomials lies in the fact 
that they form a complete orthogonal set (see Section 75), 80 
that arbitrary functions can be expanded in terms of the 
polynomials. We shall show in this section that the functions 
{P,(2)} are orthogonal in the interval —1 <x <1, and that 


Ji, PalePae) ae 


0 men nm =O), 204 
1 2 ee) 
fi, Pete) de = BET AOL 


To prove the orthogonality conditions, recall that if P,(s) and 
P_(z) sre two Legendre polynomials they satisfy Legendre’s 
qustion, 

zlo -»£r + n(n + DP, = 0, 
and 

[a -29 GPa] + mint WP. = 0, 
respectively. Multiplying the first of these equations of P,, the 
second by P,, and subtracting, we obtain 


da ad 
PEI — =P) — Pa Zl — Pd 
+ mm t 1) — met DPPe =O. QT) 
But upon integrating by parts, 
1 Pe Glll —29P olde = PP. [)— [! 0 29 Pade 





=f) 0 — PAP de, 


* If we define 
d+ 
Pal) = 4. ¢ Pa) 


the set of functions {p,(z)} will form an orthonormal set in {—1,1], and we may write 


Ji, Pcerm(E = Ayu m=O, 12 ore 
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and integration of Equation (271) yields 
= fd PP ds + fh oP Pn de 
+ (mt +m — nt — 2) [!, PaPy de =0. 
If we assume that m # n and hence m* + mn? + n, 
JL Pa@PG) de = 0, men 


To prove the second of the orthogonality conditions, we assume 
known a property of the Legendre functions which is proved 
in the next section, namely, 





1 < a 
Woe 2 PC). (272) 
If we write Equation (272) for another value of A, say 4: 
= DorPale) 
= 


and multiply this by Equation (272) we have, on integrating from 
—lto +1: 








f ats dx 
+t (1 = 2px 
Zere [nore e7) 





‘The integral on the left-hand side of this equation is (see [22], 
Formula 113) 
1 1+ Vie 
sh ae 
Vu 1 Vm 
If we expand the above function in powers of Au, 


1 gg lt Vm oJ eg = 
Jig Vie ~ vag oe + VID og (1 — VAn)] 


De Mt at 
(r+ anor) tat ey +o) 
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and equate powers of Aw in Equation (273) we see that 


[rade = hye 


We have snown that 
Ji, Pa Pa) de = 57 Bam (274) 


‘The orthogonality of the Legendre functions is therefore 
established. The completeness of the {P,(2)} will not be proved. 





79. Generating function for the Legendre polynomials 
Consider the function 
FQ) = 





1 
Vi =e +o 
which arises in the mathematical theory of potential. If we expand 
it in a power series in A, 


Fila) = eR), @75) 


where the coefficients of the powers of ) are functions of z. We 
shall prove by expanding F(2,) by the binomial theorem and 
then equating powers of A in Equation (275) that the R,(x) are 
precisely the Legendre polynomials P,(x). The binomial expansion 
of F{A,x) is given by 











Fda) = 1 — S08 —D) + 5(53) OF — Dat 
1 (1-3-5 
= xta($33) @ - ae ae 


4 GM U3-5- 24 — Mae rrop 





* From the power series expansion we could also conclude that 


f PalOPn(edde = 0, nem 
since the coefficient of 1".™ (nm > m) on the left-hand side of Equation (273) is zero. 
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Given an integer B, every term of this expansion for which 
B2=nZBi2 contains the term . For example, the f’th term 
of the above expansion may be written 


(Gh O28 3 0P 0 de () (22-41) * 


and the coefficient of A in it is 


(3-5-8 Dh apap a 1:35 B-Day 


Similarly the (6 — 1)st term is 
GP 3 5B — Mos § (°') (03) 4, 
and the coefficient of A in this expression is 
gre 3+ -3) Cy ')onxys 
1-3-5... (28 —3) 
=—“ga-a 


The coefficient of 4° in the (8 — 2)nd term may similarly be shown 








(Px osx 
135.9 
= HB — 91 ; 

Continuing in this fashion we find that the total coefficient of 

? in the expansion of F(2,) is 

RAe) = PAs oD [# - oe 


BG —) G-) 8-3), 
+ “mag 1) en —3) 
and comparing this with y*,(x) of Equation (268) shows that 
Res) = Pyz). 





Bl 


+ (f) inthe binomial ceticient: (8) = Wom 


ke 
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80. Completion of the solution of Section 72 
We considered in Section 72 the problem of the potential in 
a sphere which led to the partial differential equation 
sin 9 2 (3) +8 (ia 0%) =0. 276) 
By means of the substitution 
4.8) = RCO) 


we reduced Equation (276) to the two ordinary differential 
equations 
AR’ + WR — nn + 1)R = 0 


£[¢-9S + nla + 1)@=0, = cos 0. 


‘The first we recognized as a Cauchy equation, which reduces to 
Ri + R— nln + I)R=0 
by means of the substitution r = e'. Its general solution is 
er) = Cyent + Cerne 
or, in terms of r, 
Re) = Cy + Cy, 
"The second equation is a Legendre equation whose general solution 


is given by 
(6) = Crys (coe @) + Cys (cor 0). 

[We use the solutions (x), |x| < 1 since |cos 6] <1.] When@ =0 
or =, (the z-axis), |cos 6 = 1 and the only solutions of Legendre’s 
equation which are finite for |cos 6] — 1 are the Legendre 
polynomials P,(cos 0). Also as r—>0, r-(*#1)-+ o, which is 
physically unreasonable. Hence the only solutions which remain 
finite are of the form 


ufr,0) = mP{c08 6). 
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Any linear combination of these u, is a solution of Equation (276), 
and hence the function 


ulr,8) = Sa,mP,(cos 8) 71) 
= 
satisfies the partial differential equation, whatever the values of 


the constants a,. The a, can be determined from the boundary 
condition giving the potential distribution on the surface of the 


sphere, 
WA.) = f(6), 

which we find convenient to write in the form 
(8) = g(cos 8). 


By means of Equation (277) this condition gives 

(4,8) = g(cos 6) = Dy a,A"P,(cos 6). 
Due to the orthogonality of the Legendre polynomials, we have, 
‘on multiplying both sides of this equation by P,(cos @) and 
integrating from —1 to +1,* 


Ji e@)Pals) ds = Ara, | Pats) de = 





mee 


Pg J eavP ale) de. 


The final solution of the boundary value problem is therefore 





Oy = 


ary > [lan + nf etree ax] (4)'Pucoee). 278) 


‘As an example, suppose that the initial potential distribution is 


SO) = g (cos 8) = sin (m cos 0). 
Then 
az) = sin zz, 


* For simplicity we have set cos @ = s. 
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and 
ay = +f! Pele) sinzsds = +f sinwz de = 0, 
a= as P(e) sins ds = Br Ji ssinmede = 3. 
a = shy f', Pas) sinnads = 35, 
an = yyy f' Pale) sinweds = sv S', [3#-3] sins de, 
cok 108 
AS ~ AAP 
Hence, by Equation (277), 
10: e 
00) = 25 P(e — (B-2) Fro +... 


81. Solution of Bessel’s equation 
If we write Bessel’s equation (Section 73) in the form 
ay" + ay’ + (# — nly = 0 (279) 


we find by the method of Frobenius two linearly independent 
solutions 


a af 
ahs -[! ~ atl + wT NO? 





2 
ether y@sy + ~] 
and 


2 # 
ale) = =! aay + 202 =a =m) 


2 
— wad) 2—yG—y + =f 


If nis not an integer, both y,(2) and y,(x) converge for all values 
of x. If m is a positive integer, y,(x) converges but y,(x) does 
not; if m is a negative integer, y,(x) converges but y,(x) does 


4 
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not. Finally, if n = 0, y(x) = y4(2) converges for all values of x. 
We thus see that if n is an integer only one of the above functions. 
is a solution. The other linearly independent solution cor- 
responding to ,(#) may be proved to be (see [25], Chapter VIII) 


wt) [Sey (2%) 
with a similar expression corresponding to y4(x). 
‘The function ) 
we 
IM*) = BED) 
is called the Bessel function of the first kind of order n. The function 
Jats) = safth (281) 


is called the Bessel function of the first kind of order —n, and 
is linearly independent of J,(z) if mis not an integer. We see, 
therefore, that if m is not an integer, the general solution of Besse!’s 
equation may be written 


WE) = Cin(*) + Ca) 


where C, and C, are arbitrary constants. If m is an integer, the 
second linearly independent solution of Bessel’s equation [Equation 
(280)] is known as the Bessel function of the second kind of order n 
and is generally indicated by Y,(x). Hence if m is an integer, 
the general solution of Bessel’s equation may be written 

We) = CaIn(2) + C¥ale)s (282) 


where C, and C, are arbitrary constants, 


82. Properties of the Bessel functions 


We shall obtain in this section some properties of Bessel 
functions necessary to prove their orthogonality (see the next 
section). Other identities and relations are left to the problems 
(ece end of chapter). 
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With the aid of the gamma function, the Bessel functions of 
the first kind of order +n may be written in the compact form 


SC [al 
10) = Bana (3) 83) 


" =o 
didi= Sankveen(3) ; oe 


When n is a positive integer we have (see Section 76) 


Ty-4t+atl)=@ for «= 0,1,...,8—1. 
Hence 


Tle) = Lard are a(t a 


= Se (2) 2 Bu (eye 
Zate- Ta(d ) ~ feb érnm(z) 

on changing the dummy index of summation from ato B = a —n. 

But if m= 0, 1,2, 3,... then J,(x) as given by Equation (283) 


becomes 
&(-1y fe) 
Ie) = Boe P5(3) 











Hence 
J-n®) = (AAP) # = 01,23... (285) 


Another formula to be used in Section 83 is obtained by 
multiplying Equation (283) by +* and then differentiating with 
respect to x: 


d ag (ye aintte 
aMLO =F Laorerepy FE 


-- ete 
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But if we replace n by n —1 in Equation (283), 


_ Bw ta) (x) 
= ones (G) 
Hence 


Ae nel = tha). (286) 


The Bessel function most frequently encountered in the 
applications is J,(z). It is convenient to write out explicitly its 
power series expansion : 


Bi -ip7 wie 
1 = 3 Sar(3) 
='-h4+ eae (287) 


Its graph has the appearance of a damped cosine wave. [See 
Figure 79 for an accurate plot of Jy(z).] The other linearly 
independent solution Yq(x) of Bessel’s equation of order zero is 
sketched in Figure 80. 
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83. Orthogonality of the Bessel functions 


It is clear from the graph of Figure 79 that J,(x) has infinitely 
many positive roots, which we shall denote by A,, Ag, Ay» «.. where 
0<’,<’,<.... We shall now show that the functions 





ve h(44)| i= 1,2,3, 
are orthogonal in the interval {0,<] and more particularly that 
freh(at)1(aZ) a - Suaume,+ ese) 
Before starting the proof of Equation (288) we remark that 
the above statements can be generalized to Bessel functions of 
the first kind of all orders n 2 0, viz., 
(2) Ja() has infinitely many positive roots, say 
DAM AM AM <.... 
(>) The functions 


J(ar) 


are orthogonal in the interval [0,c] with respect to the weighting 
function =, and 


Soxda( eZ) (mE) de = Soe Qey 


To prove the orthogonality of the Jy functions it is convenient 
to introduce the functions 


(2) = VeJdox), os) = Vej{Bs), (289) 


* To aay that the functions | VJ) are orthogonal in the interval (04) 
is equivalent to saying thatthe functions J,( 4) are orthogonal in the interval 
[lic] with respect to the weighting function x. 
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where a and f are two constants. Since Jy(x) satisfies Bessel’s 
equation (Equation 279 with n=), Joax) must satisfy the 
‘equation 

(xy gay Thos) + (08) ge Idan) + (os Ide) 


= 8X jos) + 2 E Jdes) + atstflea) = 0. 
Substituting in this equation the value of Jy(xx) a8 given by 
Equation (289), we obtain 
aetluy” + xMuy + ate{etly) = 0, 
or 


deta” + (dota + Iu = 0. (290) 
It may be similarly proved that 0 satisfies the equation 
4xto” + (4B + 1)o = 0. (291) 


If we now multiply Equation (290) by o(x) and Equation (291) 
by u(x) and subtract, we obtain 


AsX(ou" — uel”) + de%(ctuo — frou) = 0, 


or 
wo” — ou” = uefa? — 6%), 22) 
and integration of this equation between the limits of O and + gives 
Jeo" — ow") de = wo’ — ow) — fF u'o' de + [Pow ae 
= (£4) f gue de, 
or 


GP = BY JF ap) = wade") — ofa '(0)] 


— [A OV(0) — of0yu(0)). 
Recalling the definitions of u and v we see that 

0) = 0- J40) = 0, 

(0) = 0- Jo{0) = 0, 
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and hence that the above equation becomes 
8) fF Sde ABO dl = x Joa) JulBs) + 4 Jes) JdBs) 
= JB) I'dee) — 5 J4B*) Je) 
= Aloo) Jd) — Jd)’ eed]. (293) 
But from Equation (286) with n = 0, 
Aue) = 4), 
and from Equation (285) 
‘ie Ja) = — 10. 
EZ ldex) = —a(ox) 
Z labs) = BIB»). 
Utilizing this result in Equation (293), we obtain 
#8) fF Idea JABD do = afae(Pa) Yon) — Blox G)- (294) 
If we now let x=, «= Ne, B = pc, Equation (294) becomes 
or 49 512) nwt) ac = &[Areone —£10rn09]. 
and if A and p are two distinct roots of Jo(), 
Jyato(r2)to(w2) ax =0 Atp (295) 


‘This proves the first orthogonality condition. If A= are the 
same zero of Jo(x)=0, we cannot divide the equation, as we 
did above, by 4*—y, But by differentiation of Equation (294) 
with respect to a we have 


de fF Ble yTABE) ae + (at BY ft A Ilo yFABD ae 
= +[ MB Nos) + MBO Z Mes) — BS J4es)i8s)]. (295) 
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If we let x — ¢ and a — B= Ne, where d is a root of Jy(x) — 0, 
we see that a?—£*—0 and hence the second integral on the 
left of Equation (296) vanishes. We also notice that 


Ide) = Jor) = 0 
and that by Equations (285) and (286) 


A eee) = eI-slee) = —eh(oe) = —eI,0). 


2 fone(is) w= [3 


or, since A 0, 


Hence 





—44001,0)] 


. frata(Qe) a = Seon 7) 


which complete the orthogonality proof. 

The roots (or zeros) of the function J,{x) are often needed 
in the applications; the first twenty are tabulated in Table II (see 
also Figure 79). 






‘Taste IT 
Zenos oF THE Besse. FUNCTION Jo(x) 
A= 240s / | Au = 33.716 
Ar = 5.520 | Au = 36.917 
As = 8.654 | Aus = 40.058 
Ag = 11.792 | Au = 43.200 
As = 14.931 Aw = 46.341 
Xe = 18.071 | Dae = 49.483 
As = 21.212 | Ary = 52.624 
Ay = 24.352 | Aue = 55.766 
Ay = 27.494 | Aus = 58.907» 
Are = 30.635 | Aro = 62.048 


84. Completion of the solution of Section 73 
We considered in Section 73 the problem of finding the 
temperature in a thin disk insulated on the lateral surfaces and 

with the boundary conditions 
1. xa,t) = 0, 


2. ufr.0) = fir). on 
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‘The equation satisfied by temperature u(r,6) was found to be 
Ow law 1 


wt re kao ie 
and assuming a solution of the form 
Hr) = ROT, 
we arrived at the two ordinary differential equations 
T’ + kB'T = 0, (3002) 


Ru + te + BR = 





(3005) 
By means of the substitution s = fr, the second equation reduces 
to the standard form of Bessel’s equation of order zero : 
®R aR 
eae teG + aR = 0, 


whose solutions are Jy(:) and Y,(+). Hence the solutions of 
Equation (300b) are Jy(r8) and Y,(7B), and the function 

urst) = MDB) + eYLeB)] (301) 
satisfies Equation (299) whatever the values of the constants B, 
b,c. But referring to Figure 80, lim Y,(2)=—s, and sinée 


1(0,*) is certainly finite, we must eliminate the function Yy from 
our solution. Hence 


ur,t) = be¥*Je(rB). 
From the first boundary condition of Equation (298), 
wat) = be-** J(aB) = 0. 
This implies J,(a) = 0 and shows that af must be a root of 


Id#) = 0, 
B=rAm PA §= 12.05 


Due to the linearity of our equation, we conclude that the function 


). 





a¢rt) = BS byerterlere Jobe 
Zz 
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where the 6, are arbitrary constants, is a solution of Equation (299) 
which satisfies the first of the boundary conditions of Equation 
(298). If u(7,#) is to equal f(r) when 1 = 0, we must have, moreover, 

fe) = ¥ budo( WZ). (302) 


‘We now obtain the 6, by means of the orthogonality properties 
of the Bessel functions. Multiplying Equation (302) by 


h(n £) 


and integrating from 0 to a with respect to r, 
f ON dae) im z +f. 1&5) bz) a 


SUM 





from which, 
= 2 ft # 
ba = ayaa Seton (E) & 
‘The solution to the temperature problem thus becomes 
28.1 [pe ‘ Paw 

ste) — BB jay [Sfton( eed) a] 2g )reen s 
Infinite series of the type appearing in Equation (303) are often 
spoken of as “ Fourier-Bessel” series in analogy with ordinary 


Fourier series. 
As a numerical example, assume f(r) = 1. Then 


2 @ r z 7. r 
bs = arpony JoehoM 5) = apy SordMS) 
If we make the change of dummy variable of integration, 
sat 


ony 


2. 1 - 
be = jaca TF Le 2) de = ray 


Google 


212 LUGENDRE, BESSEL, MATHIEU FUNCTIONS Art. 85 


[See Eq. (286)]. The numerical result is therefore 


wea) = 2B payrnnen Hae): 


aot 


‘Values of the Bessel functions are tabulated, for example, in [30]. 


85. Solution of Mathieu's equation 
‘The ordinary differential equation 


¥" + (A—2g 008 2a)y = 0 (304) 
was called a Mathieu equation (see Section 74). If we let 
=f 


in Equation (304), it reduces to 

¥" — (A—2g cosh 2xly = 0, (305) 
which is another form of Mathieu’s equation. It will therefore 
suffice to study Equation (304). While Mathieu’s equation may 
appear simpler than Legendre’s or Bessel’s equation, we note that 
it depends on two parameters, A and g. Hence certain features 
not encountered in these earlier equations will be present in its 
solution. 

We note first that if g = 0, Equation (304) reduces to 


y” +d =0, 





so that with A= m?, cos mx and sin mz are two linearly 
independent solutions. Hence we should like any solution we 
obtain of Equation (304) to reduce to cosmx or sinms when 
q¢=0. We shall scarch for periodic solutions y(x,A,g) (of period 
7 or 2n) of Mathieu’s equation. If »(x,A,g) is one such solution, 
and q must be interrelated, say 


Ds mB tb gg tag? + ag! to... (306) 


Since we want A to equal m* when q = 0 we have chosen a as m*. 
Solutions of the form 


(,A,9) = cos mx + ga(x) + gar) + gras(x) + 
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reduce to cos mx for g = 0. These solutions are generally denoted 
by 

Semn(*9)s 


where “ce” is read “ cosine-elliptic.” Similarly, solutions of 
the form 


sin me + gb,(x) + g%(x) + gba(2) + -.. 

reduce to sin mz when q = 0, and are generally denoted by 
em(%9)> 

where the “se” is an abbreviation for “ sine-elliptic.” Both 

ce,(zq) and seq(2,9) are referred to as Mathieu functions of the 

first kind of order ».* ‘The first few cosine- and sine-elliptic 

functions are listed below. 

cndeg) = 11g om 2 + et eo te 


= Gt (Foo 6 —7 oo 2s) tos 


ce(%q) = 008 = — go 00 Be + deat oo Se — cos 34) due 
(307) 





coleg) = con 2x — Fa( Foo 4 —2) + sgt on Ge 2 


* Since the value of A is different for cey and seq these functions cannot be 
solutions of the same Mathieu equation (except in the case q= 0). ‘The other 
independent solutions of Mathieu's equation, fe,(x,9) and ge,(z,¢) corresponding 
to cey(x,9) and se,(x,q), respectively, are nonperiodic. The general solution of 
Equation (304) can therefore be written either as 


Wx) = A ceq(xa) + B feq(xa) 








ya) = C seg(x2) + D gen(*ia), 
where A, B, C, and D are arbitrary constants, and n is an integer. (For details 
of these other solutions see [20].) It will be seen that as in the case of Legendre 
‘and Bessel functions of the second kind, the fe, and ge, functions do not often 
appear in the solution of elementary boundary value problems. In particular 
‘we shall not need them in our problems, 
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: Ll. 1/1. ” 
se(xg) = sin x — pg sin 3x + 9%( 7 sin Sx + sin 3) — ..., 


seg) = sin 2 — foaten 4 + ge? sin 6 Padiad 


‘The functions ce,(x,g) and se,(x,9) are periodic of period 2m if 
nis odd and periodic of period 7 if m is even. We may therefore 
write 


cen(#9) — 2) Asa) cos 2ax, 
a 


Cen e(9) = Ztaanen cos (2a + 1), 
(308) 


sem(*g) = x B,p™ sin 2B, 


Seansn(%9) = Bro sin (2B + 1)s, 


where the 4’s and B's are known constants depending on the 
parameter g. 
We have seen that Equation (304) can be transformed into 

Equation (305) by replacing x by jx. Hence 

cen(jxg) and —_se,(J2,9) 
are solutions of Equation (305). The solutions of Equation (305) 
are generally written 

Cex) and — Se,(x), 
where 


Con(aig) = cen(ing)» 
(309) 
Sen(x) = —j sea(iz9), 


and are called modified Mathieu functions of the first kind of integral 
order. 
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86. Orthogonality of the Mathieu functions 
‘We shall now show that the functions {ce,(x,q), seq(x,9)} are 
orthogonal in the interval (0,2m). Let 


Jala) = cen(%q) and yal) = ce,(x,9) 
be solutions of the Mathieu equations 
o” + (x —2g 008 2x)y = 0 G10) 
»" + Os —2g 008 2s)y = 0, Gl) 


respectively, with different A's but the same g. Multiply Equation 
(310) by y, and Equation (311) by y, and subtract, This yields 


Weds — W's = On — dae 
Integrate this equation from 0 to { : 


Jolrstedy's00) — salads] de = Oa — Ad fOyvladrle) a. 


If we integrate by parts as was done in the corresponding 
development for the Legendre and Bessel functions, we obtain 


Deve) — 9's. = a —A) [Ennledvde) de. 12) 
Since y, = cen(x,g) and yg = ce,(x,g) are both periodic of period 
zr or 2n depending on whether # and m are odd or even, the left- 
hand side of Equations (312) vanishes for £ = 27, and 

Qa —A) fp” cenl sn) ceq() de = 0. 
If mn, ~A and 
IP" conta) centers) de 
Tf, instead = m and n is odd, say n = 2r + 1, 


Jo ects de — fo" [3S ase oe (an + 1] a 


” 2 ark #0 


W 


and 





0 men 
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by the orthogonality of the cosine functions. Similarly, if » is 
even, say = 2s, 


[Ol ectten) ae = rotate + oS asey = 0. 
% 
Using the same technique, we can readily show that 


Jo s<n(esg) seg) de = 0, em, 


Jo" sehxa) de = ¥ er = 0, 
a 


Je sehen de = oS aR = 0, 
a 
and 
J." cenlte) sea(eg) dz = 0. 


Note that by virtue of Equation (308), this last equation holds 
whether the q appearing in the cosine-elliptic term is the same 
as or different from the q appearing in the sine-elliptic term. 


87. Completion of the solution of Section 74 


In Section 74 we considered (Equation 245) the problem of 
finding the transverse displacements of an elliptic membrane and 
saw that a function u(x,y,¢) of the form 


u(z,y,t) = a(x,y) (A cos wt + B sin wt) (313) 


was a solution of Equation (245) if w satisfied the equation (in 
elliptic coordinates) 


Bw Aw 
Se + SS + wtcosn 2g —cos nye = 0. (4) 


The membrane was assumed to be an ellipse of semimajor axis a 
and semiminor axis f, and the focal distance y (see Equation 247) 
was so chosen that the membrane became one of the confocal 


Google ' 


Art. 87 LRORNDRE, BESSEL, MATHIEU FUNCTIONS 217 
ellipses, that is y2 = a? — f (see Equation 251). It was that 
particular ellipse of the family corresponding to (sce Equation 248) 


£ = & (aconstant), 
where 


tah, = £ dma tanh 2 
‘The constant k in Equations (249) and (250) was given by 
_ ad — ft) 
ae He 15) 
Since the membrane was assumed fixed on the boundary with 


an initial displacement and zero initial velocity, the boundary 
conditions, in the elliptic coordinate system £, 7, become 


1 wont) = 0, 
aa att G16) 
3. 5 MlGn0) = 0. 
‘The solution of Equation (314) led to two Mathieu equations 
[Equations (252)], whose periodic solutions were found to be 
Cea(fg) or  Sen(€,9) 


cen(mg) OT 8 nC 9)» 
respectively, where 


and 


m= B G17) 


‘The solutions of Equation (314) must therefore consist of terms 
of the form 


Cen(E,9) cen(7.9) 
Se(E,9) sen(7.9), 


where Ce, and ce, correspond to the same value of A, while Se, 
and se, also correspond to the same value of A, but not the value 
appearing in Ce, and ce,. Hence the function 
un(Emt) = Cen(E,9) cen(9) [An 008 Wat + By sin cyt] 
+ Sen(Esg) 8ex(7,9) [Cy 008 wryt + D, sin cont] 


and 
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satisfies the partial differential equation of Equation (245). We 
may use the third boundary condition of Equation (316) to eliminate 
B, and D,. For, 

ZB ulGnt) = ColEa) conlnal—tinAy sin wat + eonBy 608 cf] 

+ Seq(E,g) eal mg)[ ergy Bin cage + cyDy C08 cyt] 

and 

0 = J u(Eins0) — Bytom Coal) coals) + Datta Se) seal) 
Since the elliptic functions are not identically zero, 

B, =0= Dy 
and we may write 
(Eins!) = Ay Cos(€a) cose) 208 ent + Cy Seal leg) eey(n) £08 ent. 
From the first boundary condition of Equation (316), 
Hq(Eort) = An Cen(Gor9) cen(n9) 008 ergt 
+ Cy Sen(Eorg) 8€n(79) 008 wryt = 0. 


Since A, and Cy are arbitrary and coa(74) #0, coe wat #0, we 
must have 
Con(Eont) = 0. (318) 


It may be shown that corresponding to each n, there is an infinite 
number of roots, say 
Gas Gu Yass oe 


‘of Equation (318), and for each gam we obtain from Equations (315) 
and (317) a frequency am : 


x, erin d, (319) 
Similarly we obtain from 
Sen(Eng) = 0 (320) 
an infinite number of roots, say 


Gar Gee Tar + 
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and the corresponding Dm: 
Tat = She. eat) 
We therefore see that the linear combination 


“ent = & YA Ceq(6,tem) C60 am) 608 Cyt 


=o me 


+ 3 > Com Seal E Gam) 82x(7Gnm) C08 Ginmt (322) 
wom 


satisfies the differential equation and the first and third boundary 
conditions of Equation (316). It will also satisfy 


(.n.0) = Aén) 
if the A,, and C,,, are so chosen that 


6:00) = En) = 3 then Colton) cn) 


+ Cum Seq(E nm) 2¢n(7Fmm)]- (323) 
‘The orthogonality conditions proved in Section 86 are useful 
for expanding an arbitrary function in terms of sine-elliptic and 
cosine-elliptic functions. However, these orthogonality conditions 
were based on the assumption that the se, and ce, functions 
corresponded to the same value of g. Clearly the functions of 
Equation (323) have different values of g. Hence, in order to 
evaluate the A,, and Cy, we must prove another type of 
orthogonality condition involving Mathieu functions. 
We have seen that 


Cen(E Gam) Cn 2am) (324a) 
and 

Seq(E,dam) S€x(%nm) (3246) 
are solutions of the equation 


> + at 2ym(Coth 2£ — cos 2r)o = 0. 


Google 


220 ‘LEGENDRE, BESSEL, MATHIEU FUNCTIONS Art. 87 


Hence, calling vyn(£im) and v,,(én) two of the functions of 
Equation (324), 


am + oo + Unm(cosh 2f — cos 2n)onm = 0, (325) 





+ + 2plcosh 2£ — cos 2nJor, = 0. (326) 


Multiplying Equation (325) by »,, and Equation (326)' by eam 
and subtracting, we obtain 
Wag a: = agg BiH ag athe = Ha 
ra GEEeom — Pam ap Pre + Ore Gr Pam — Pom Bae Pre 
+ 2Udnm — Gra) (Cosh 2E — c08 27)mtrs = 0, 
or 
a, @ a af, 2 a 
Rid asad tA oan actos 
+ dan — dra) (Cosh 2E — 08 2n\vamtrs = 0- 


Now integrate both sides of the above equation from 0 to 27 
with respect to 7 and from 0 to €, with respect to €. 


Jo (ea fetan— ton 0)" or fo (end nm Pm Ze) 


+ 2am — a1) J2" fo" YanrAcosh 2 — cos 2n) df dy = 0. 


ia0 
dé 
lo 





‘The first integral on the left vanishes at ¢ = 0 (see Equation 308 
and Section 86); the integrand vanishes at ¢ = , since the 
membrane was assumed clamped at =; and the integrand 
of the second integral vanishes since v,, and t._ are periodic 
of period 2m in 7. Hence we conclude that if dam Yen 


ir i VnmBes(Cosh 2¢ — cos 2n) dé dn = 0. (327) 
If n= 1, m=, dam = Yen and 


JR Je emmt&sn) (cosh 2 — cos 2n) dé dm 24 0, (328) 
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since P_,2 > 0, and for any € and any 9, 
cosh 26 — cos 2 = 0. 


By virtue of the orthogonality conditions of Equations (327) 
and (328), multiplying both sides of Equation (323) by 


Ce,(E tre) ee(7HGr1) (cosh 2§ — cos 29) 


and integrating from 0 to £ with respect to £ and from 0 to 2 
with respect to 7 we obtain 


JR Jet A600 Cox€.tn) conan) (cosh 2€ — 008 2n) dE dy 


= Aa Jer fot CeX€sted conte) (cosh 2 — con 2n) df dn. 
Hence 
‘ 
WP” Jot AGsm) CenkEstes) cerknites) (comb 26 — cos 2m) df dy 


IP [ee CorM€an costar) (cosh 26 — eos 2n) a dy 





Similarly, 
JR Jet Abn) Seb) seen (cosh 26 — 000 In) db dy 
WP Jot Ses 


We have thus determined the Ay, and Cy, of Equation (322) 
and completed the solution of the boundary value problem of 
Section 74. 


) se.Mndr) (coh 26 — cos 2n) df dy 





88. A unifying principle 

‘We have considered in this and earlier chapters four sets of 
orthogonal functions : the trigonometric functions, the Legendre 
polynomials, Bessel functions, and Mathieu functions. All these 
functions were solutions of second-order linear ordinary differential 
equations with certain appropriate boundary conditions. It may 
be reasonable to ask whether there is any general theory of such 
equations of which the above four examples are merely special 
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cases. The answer is yes, and is found in the general Sturm- 

Liouville theory of ordinary differential equations (see [12]). This 

theory considers a differential system of the form 
é [no gas] + [ole) + Me(e)]a(e) = 


ca(a) + fau'(a) = 0, 

y(t) + Bu) = 
defined over the interval (a,b), where A is a parameter, (Actually 
the theory embraces a slightly more general differential system 
than the one considered above.) 

Any second-order linear differential equation with a parameter, 
namely, 
Flew" e) + eleml(e) + He) + bCayIa(z) = 0 
may be written in the form of Equation (329) on multiplying 
by the integrating factor 
1 fitter 


Fey 
For then this equation may be written 


[Sve] + Hae + ad ay = 0, 





(329) 


i) 
and on letting 
we) = J, 
ae = SOA, 
ne) = Ahh, 


this becomes Equation (329). 

In general, the completely homogeneous system of Equation 
(329) will not have a solution unidentically zero. However, for 
certain values of the parameter A the system may have solutions 
which are not identically zero. The values of A for which there 
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exist such solutions are called characteristic numbers. The functions 
4(2) which correspond to these values of A are called characteristic 
functions. If 4(x) is a solution of Equations (329) for some A, 
then so is A(x) where A is an arbitrary constant. In general, 
a suitable normalization is applied to the characteristic functions 
which then specifies A. 

For example, the particular differential system 


wiz) + du(x) = 0, 
40) = 0, un) = 0, 


over the interval (0,7) has, in general, no solutions which are 
not identically zero. However, for 


=n 


where n is an integer, Equations (330) do have a solution 
unidentically zero, namely, 


(330) 


sin ax. 
Hence the set of squared integers n* are characteristic numbers 
for the differential system of Equation (330) and the {sin nz} 
functions are the corresponding characteristic functions. 
‘The fundamental theorem associated with Equations (329) may 
be stated as follows : 


TueoreM. Let p(x) p’(a), (2), r(x) be continuous functions 
in the interval a<x<b. Let p(x) >0 in [a,6]. Then the 
Sturm-Liouville system 

[oem + [ale) + (a)]}u(x) = 0, 
cafa) + fu'(a) = 0, 
db) + 8u'(8) = 0, 

has an infinite number of characteristic numbers, 
Ac cdc... 


and corresponding characteristic functions {4,(2)}. Any function 
‘flz) which is continuous on [a,b] and has a derivative everywhere 
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except at a finite number of points in [a,b] can be expanded in 
an infinite series of the characteristic functions {44(x)}; 


fe) = Seats 
where 
a 
6 = fo rayferb(a) ds, 


and this series converges uniformly to f(x) in any subinterval of 
[a]. 

We have seen that the trigonometric functions satisfy the 
equations u + Au = 0 and that A, =n are the characteristic 
numbers, $,(x) = sin mx, the characteristic functions. Also, the 
sine functions are orthogonal in the interval (0,7). 

Legendre’s equation, Equation (240), may be written 


& [0-9] + me) = 0 


Here p(x) = 1— 2%, g(x) = 0, r(2) = 1. We have seen that the 
characteristic numbers 2, are n(x -+ 1) and that the characteristic 
functions P,(x), the Legendre polynomials, are orthogonal with 
respect to the weighting function r(x) = 1. 

Also, Beseel’s equation of order zero, Equation (243), may be 
written in the form 





A towed + daw) = 0. 


We have scen that for the interval [0,1] the characteristic numbers 
Aare A,!, where the A, are the roots of the Bessel function J,(x) 
=0. Also, the functions J(A,x) are orthogonal with respect to 
the weighting function r(x) = x. 

Actually, a slight modification of the Sturm-Liouville theory 
has to be used in the case of Legendre and Bessel functions, 
since the condition p(x) >0 of the theorem is not satisfied in 
these two cases. 
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PROBLEMS 


1, Derive the expression for the Laplacian 
Ou + Au += 
aa + os + ae 

in spherical coordinates. 

2. Derive the expression for 

Vio + Yo = 0 
in elliptic coordinates. 

3. Show that the three-dimensional heat equation 
au eu Ou ou 
gett mr eg 

‘becomes 

eu 1 ow au 1A Ou 
wt ret ott oe Oe 

in cylindrical coordinates (r, 6,’z). Sce Figure 81. 
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4. Determine the constants a4 By, yx 80 that the functions 
As) = ay 
Bee) = Br + Ber, 
Ce) = nn + ye + ys, 
become orthonormal in the interval (0, 1). 


5. Show that the functions 


1 cosx sin cos 2¥ sin 2¥ 


Vin’ va Va Vn 
form an orthonormal set in [— 7,7]. 








6. ‘The beta function, defined by the integral 
Bue) f el —x)1de, wo > 0 
often appears in mathematical physics. Prove that Bu, ©) = Bo, »)- 


7. Prove that the beta function is related to the gamma function 
by the equation 





Y 


[Hint : Make the substitution 1 x = sin™@ in the expression for the 
beta function and the transformation x = rcos*6, y = rsin®@ in 


TW) Te) = fo fP eter yn ter de dy] 
8. Prove by the result of Problem 7 that 
1 
M4) =v 
9. Show that 





where m is a positive integer. 
10. Derive, in detail, the solutions of Legendre’s equation 
(1 — aly” — 2ay’ + n(n + yy = 0 
for |z| > 1 and |x| <1. 
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11. Show that P,(2), the Legendre polynomial of order n, is an odd 
function if m is odd and an even function if n is even. 


12. Prove that 


Pl) = 1, 2 = 01,2... 
13. Establish the formulas 





Pe) = Ge 


Ett pie sa 
15. Establish Rodrigues? formula 
PAs) = sig Sa — 


[Hint : Integrate, m times, the expression for P,(+) (Equation 268) from 
Os] 


16. Show that the Legendre polynomials satisfy the recursion formula 
mPq(%) — (20 + 1)8P,(2) + (+ 1)Paia(#) = 0. 


(Hint: Consider the generating function of Equation (275) for the 
Legendre polynomials and differentiate it with respect to 2] 


17. Establish the recursion formula 
Pian) — P’eal®) = (2 + IPA). 


[Hint : Differentiate the generating function of Equation (275) with 
respect to =] 


18. Prove that 


flay de = 5 ae) — Pentel 
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19, Establish formulas for 
(@) —  #Pae) as, 
(bf em Pels) de, 
(©) #P,G) de. 
20. Generalize the results of Problem 19 to 
Jo ePala) ae. 
21. Prove 
fi, AerPate) de = SUP ft cat — 1 FE) ae 
[Hint : Use Rodrigues’ formula.] 
22, Expand the function f(z) = »* in a series of Legendre polynomials. 
23, Expand the function f(x) 
fe =-1 -lezr <9, 
fa)=t1  Ocrcl, 
in a series of Legendre polynomials. 
24, Expand the function 





fe) = —<s#<9, 
fe)ste Oce<l, 
in a series of Legendre polynomials. 


25. Obtain the Legendre series expansion for the function 


fe =0 -1ex#<G, 
fe) 2  Ocw<l 


26. Let a sphere of unit radius have its lower hemisphere kept at 
a constant temperature U,, and its upper half kept at zero temperature. 
Find the steady-state temperature distribution u(r, 6) inside the sphere. 


27. Solve Problem 26 for the temperature outside the sphere. 
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28. A hemisphere has its curved surface held at a fixed electrostatic 
potential while its bounding plane is kept at zero potential. Find the 
electrostatic potential at any interior point. 


29. A hemisphere has a magnetostatic potential ys which is symmetrical 
about the vertical axis through the center of the sphere, while on the 
bounding plane, s = 0, aJ/26 = 0. Find the magnetic potential at 
any interior point. 


30. Derive in detail the infinite series solution of Bessel’s equation. 


31. Obtain the second linearly independent solution, Y,(+) of Bessel’s 
‘equation of integral order 1. 


32. Prove 
Id-#) = (API) = 1 
33. Establish the formulas 


Jule) = af 2 sins! 


1 fz 
Fouls) = a] ons = 
34, Establish the following recursion formula for the Bessel functions 
of the first kind. 
eJ'2) + Jol2) = *Jn-x(%)- 
35, Establish the identity 
a 
Ate 16@) = nla 
36. Prove the recursion formula 


Jnale) + June) = Ble). 
37. Prove that 


fF EID de = 2h (2). 
38. Calculate é 
Jy Jae) dx. 
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{Hint : Write 
fda) as = — fore] Qa)) dx 
and apply the result of Problem 37.] 
39. Evaluate 
fe 8 Jo(2) de. 
[Hint : Use Problem 37 and integrate by parts.] 


4. Expand the function 
fe)=e Ocxa2 


into a Fourier-Bessel series of Bessel functions of the first kind and 
order zero. 


Al. Expand f(s) : 
ft)=h Oce<l 


in a Fourier-Beosel series of Jo(Ass), where A, are the positive roots of 
Idd) = 0. 


42. A heavy string of length L and weight w Ib/ft hangs from a fixed 
point. If the initial displacement of the string u(x, 0) is a given function 
of x, say f(x), and the initial velocity is zero, find the displacement »(x, t). 
{Hint : The equation of the string is 

ue _ 1 oe 
"ate oe 
where g is the gravitational constant. See [25], Section 11.4.] 


43, Determine the longitudinal vibrations of a bar (Equation 14) 
fixed at the end point, » = L which is given an initial displacement f=) 
and has zero initial velocity. Assume the material of the bar is such that 
EA = kx, aA =k’, where hand } are constants. 


44. A thin flat circular disk hes its lateral surfaces insulated, while 
its circumference is held st zero temperature. Find the temperature 
u(r, t) in the disk if the initial temperature is a function of the radius 
only, 2(r, 0) = 73. 
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45. Find the temperature in a semi-infinite right circular cylinder 
whose base is insulated, whose lateral surface is held at zero temperature, 
and whose initial temperature u(r, 0) = fir) is independent of 8 and z. 


46. Find the temperature u(r, t) in the cylinder of Problem 45, if 
the cylinder has its lateral surfaces also insulated. 

47. Find the transient temperature distribution u(r, t) in an infinite 
right circular cylinder whose lateral surfaces are insulated and whose 
initial temperature is u(r, 0) = fir). 

48. A right circular cylinder of radius a and height A has one base 
and its lateral surface kept at zero temperature. Find the steady-state 
temperature «(r, 2) in the cylinder if the other base is kept at a temperature 
fir) independent of 8. 

49. Solve Problem 48 if the lateral surface of the cylinder is insulated 
instead of being held at zero temperature. 

50. A circular membrane, fixed on the boundary, has an initial 
displacement u(r, 0) = fir) and a zero initial velocity. Find the motion 
u(r, t) of the membrane. 

51, Find the equation of the transverse vibrations of the membrane 
of Problem 50 if the membrane has an initial velocity 


a 
FeO) = ae) 
and a zero initial displacement. 

52. A right circular cylinder of radius a and height h has its entire 
surface kept at zero temperature. Find the temperature u(r, 2, ) at any 
point and at any time if the initial temperature of the cylinder u(r, z, 0) 
is a given function, f(r, 2). 

53. Solve Problem 52 if the lateral surface is insulated, and the other 
boundary condition remains unchanged. 

54, A right semicircular cylinder of radius a and height h has ite 
bases and curved surface kept at zero temperature. Find the temperature 
u(r, t) in the cylinder if the plane surface is insulated. Astume that 
the initial temperature u(r, z, 0) = fr, 2). 
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55. A thin circular disk is completely insulated. Find the temperature 
u(r, 6,4) in the disk if initially u(r, 8,0) = f(r, 8). 

56. A semi-infinite cylinder has its curved surface kept at zero 
temperature, its base insulated, and has an initial temperature f(r, 8) in - 
dependent of z. Find the transient temperature distribution u(r, 0, ). 

57. Solve Problem 56 if the entire surface of the cylinder is insulated. 

58. An infinite right circular cylinder has its entire surface held at 
zero temperature. If its initial temperature u(r, 0,0) = fir, 6), find the 
temperature u(r, 8, t) in the cylinder. 





59. Find the temperature (7, @, 2) in a right circular cylinder of 
height h if the lateral surface and one base are kept at zero temperature, 
while the other base, s = 0, is kept at a temperature u(r, 0,0) = fir, 6). 

60. A semicircular right cylinder of height k has its curved surfaces 
insulated, its plane face and one base held at zero temperature. Find 
the temperature u(r, 6, 2) in the cylinder if the base (x = 0) is kept at 
a given temperature f(r, 9). 

61. A circular membrane clamped on the boundary has an initial 


displacement u(r, 6,0) = j{r,@) and a zero initial velocity. Find the 
transverse vibrations of the membrane. 


62. Find u(y, 6, t) for the membrane of Problem 61 if 


u(7,80) = 0, 2 u(r,8.0) = ar). 
63. Show the equivalence of 
yy” + (A—2q cos 2x)y 
y" — A—24 cosh 22)y 
under the transformation x 


64. Derive the cosine-elliptic periodic solutions 


ceg(ag) and ce(,9) 
of Mathieu's equation. 


i 
° 


and 


" 
° 





65. Derive the sine-elliptic periodic solutions 


sey(q) and —_seq(49) 
of Mathieu's equation. 
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66. Establish, in detail, the orthogonality of the sine-elliptic functions. 

67. Consider an elliptic membrane, clamped on its boundary with 
zero initial displacement and initial velocity a given function of € and 7. 
Find the displacement u(€, , #) of the membrane. 

68. Find the temperature in an infinite right elliptic cylinder whose 
Jateral surfaces are kept at zero temperature and whose initial temperature 
distribution is independent of the z-axis. 


69. An elliptical disk has its lateral surfaces insulated, Find the 
temperature u(€, 7,1) in the disk if its circumference is kept at zero 
temperature, and its initial temperature distribution w(é,, 0) = f(é, 9). 

70. A semi-infinite right elliptical cylinder has its curved surface kept 
at zero temperature. Find the temperature in the cylinder if its initial 
rate of temperature change 2/dt u(£, 7, 0) is a given function g(£, ) and 
the base is insulated. 
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Chapter VI 


PROPERTIES OF SECOND-ORDER PARTIAL 
DIFFERENTIAL EQUATIONS 


89. Introduction 


In keeping with the modern spirit of partial differential equations 
and the practical approach of this book, we have been engaged 
so far in the study of “ boundary value problems,” -that is, we 
have been interested in differential systems. However, we can 
obtain a deeper insight into the structure of partial differential 
equations by studying some of their general properties. Even 
though some of our conclusions may be generalized to more 
complicated equations, we shall consider in this chapter only linear, 
homogeneous, second-order partial differential equations in two 
independent variables with constant coefficients, that is, equations 
of the = 





a* fit Bay + CR te FER EO OM) 
where A, B, C, D, E, and F are constants. It will be noted that 
most of the equations studied in previous chapters belong to this 
class. 

Our first task will consist in classifying equations of the 
above form into three distinct types depending on whether the 
quantity B*—4AC is positive, negative, or zero. The proce- 
dure followed requires a transformation of the *independent 
variables in Equation (331). We shall then see what simplifi- 
cations are possible by specializing the transformations, and in 
this way reduce the equation to a form which is more readily 


analyzed. 





234 
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90. Transformation of variables 
Consider the most general transformation of the independent 
variables x and y of Equation (331), namely, 
&= ey), 1 = ney), (332) 
and compute the derivatives of u, using ¢ and 7 as intermediate 
variables : 








If we substitute these values into Equation (331) we obtain 


" » ow ou Fe oe 
AAR pete H+ PE gy t Fue 0 (333) 
where 
. ag) g aE ag\* 
a= a(S)’ + 283 + (8) 


B 


0g ag ag 3 
usp oot 4 2) £ 2c 


c=a @)’ + ane + c (3)\ 


Is 


Ox, 


2 
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Ho B  p%, pa 

Cay Ex, 
Beay + COR t Pas + Fay 


es Fon ay am. pn, pan 
bmAD + BR + Cr+ Da + en 


D-A amie 


(334) 





F=ik, 
Equation (333) is then Equation (331) expressed in terms of the 
new variables and 


91. Canonical forms 


We shall consider what simplifications can be made in Equation 
(333) under the assumption that B*—4AC is positive, negative, 
or zero. Assume for the moment that B}—44C >0. Then 
A’ and C’ may be written 


= 4(E-a8) E-a8) 


+ =e (*2 a2) (20 
o = 0 (2-43) (az). 
where A, and 2, are the roots of the quadratic algebraic equation 


AM + BX+C=0, 
and hence are real and unequal. 
Suppose we choose ¢(x,y) and 7(x,y) such that 
3 a _ om oy _ 
me R= and By = 


and the transformation of Equation (332) is reversible, that is, 
that we can solve for x and y in terms of ¢ and 7. For example, 
if we let 


$=Arty a= Arty. (335) 
‘Then certainly, 
a a 
% 8-0 9 a27= 
and furthermore we can solve for x and y from Equation (335) : 
rata ge 


b= 1 — Ae 
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since by hypothesis 1, #2, and hence 2,—, 0. With such 
a transformation 4’ = 0, C’ =, and Equation (333) assumes 
the form 
ou 
" ve E> Fu=0. 

Bet UB + Et Puno (36) 
‘We must check whether B’ ee in which case Equation (336) 
would be a first-order equation. From the quadratic equation 
AX + BL+ C=0 we see that 


B=—AQ +A) C= Ady 


{If A were equal to zero, we could interchange the roles of x 
and y. If both A and C were zero, Equation (333) would already 
be in the form of Equation (336).] 

The coefficient B’ can thus be written 





w= 42% — tay (22+22) + aa, 25] 
=~ 4 [ER nS — Bs ae SB 
~ ASH AEB + Se 
~+(C8-98) @-98) (E98) 2-98) 
-4(E-a8) 203) 
since a oe 
But 7 e 
x att Eon Ra, 
for, from Equation (335), 
1,8 
a-Az 
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Hence 
B40, 
and we may divide rieees 336) by B’ to obtain 
ee om =0, (337) 
a + teats. (37) 
where 
ee i ,_F 
tap =p foe 


Equation (337) is the canonical form of Equation (331), cor- 
responding to the condition B* — 4AC > 0. Any linear differential 
equation such as Equation (331) in which B*— 4AC > 0 is said 
to be of hyperbolic type, and can be reduced to the form of 
Equation (337). We shall examine this equation in more detail 
in later sections. 

Suppose now that B*—44C=0. Then the roots of the 
algebraic equation AX + BX+C =O are real and equal, say 
A=A,=A,. In this case the coefficients A’, B’, C’ become 





c 


i] 





If we choose a ¢ and an 7 such that 
% _ 29 = fhe 
Ez Ag xo Ge 10, 


Be 


then 
B=C'=0 and A’ #0. 


For example, the reversible transformation 
n=wty% E=% 
transforms Equation (333) into the equation 
ee - 
ABTD E+ ES + Fu =o, (338) 
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and since A’ #0, we can divide through and obtain 
Fu yOu, Oe 


mtd’S é "« = 0, (339) 

oe + a + a + fu G39) 
where D a 7 

a3 eo rank 


Any equation of the form of Equation (331) with B*—4AC =0 
can be reduced to the canonical form of Equation (339), and 
is called an equation of parabolic type. 

Finally, consider the case where B‘—4AC <0. The roots 
A, and A, of AM* + BA-+ C = 0 are in this case complex conjugate 
numbers, say A and X. If, as in the hyperbolic equation, we choose 
and 7 such that 

0G 4 i 
= - A= 0 Z -iGi=0 
the constants A’ and C’ will vanish, and our equation becomes 
identical with Equation (336). However ¢ and 7 are now complex 
variables. That is, 
& = Hay) + May) 
1 = Hy) — May) = € 
In order to remain in the domain of reals, we shall write our 
equation in terms of the real variables $ and . From the above 
equations, by addition and subtraction, we have 


1 1 
$= Eto b= HE-2. 





Hence 
B-2E+ RR - 3503) 
4G +8) 
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Introducing this second change of variables, Equation (336) 
becomes 
ahh 
4 





HH) -FB-08) 
+3 (e+ig) + Fu =0. 


Since B’ # 0 by the same arguments used in the hyperbolic case, 
we can divide by B’/4, and the above equation becomes 


we mw Tm, 
apt op tag +B Gy tye = 0. (340) 





and »/ are real. In fact, 
+e) — HE 
Se 6 PS ae 





Wo pn 
ya d 


Any equation of the form of Equation (331) with Bt—4AC <0 
can be reduced to the canonical form of Equation (340), and 
is called an equation of elliptic type. 


92. The hyperbolic equation 
We have seen that the general hyperbolic equation can be written 


Au ou a 
Be ttm hy twee (341) 


In order to simplify Equation (341), introduce the new variable w 
defined by the equation 


ulx,y) = wolayle tee, 


Then 
MM prey — boro, 
Me Beier — aur ten, 
Bs = oy ebay _@ s etry _ 5 =e enteay 4 abee 2-09, 
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Substituting these expressions in Equation (341), we obtain the 
equation im 


Be t © ake = 0. 


A solution of this. equation leads to Bessel’s equation (sce 
Problem 13 at the end of this chapter). 


93. The wave equation and the Cauchy problem 


We have considered at some length in previous sections the 
wave equation 


eu lee 
mm aa (ay 
which, if we let -y = at, assumes the form 
Au eu 
mF ao 643) 


Since A = 1, B=0, C=—1, B'—4AC = 4 >0, we conclude 
that this equation is of hyperbolic type. The roots of the cor- 
responding algebraic equation 

AM + BE C= 
are +1; hence, if we let 





€=x+y% a=x-y, 
Equation (343) reduces to 
fu 
8 oy 
This equation is readily integrated, its most general solution being 
«= 4) + Hn), 
where ¢ and y are arbitrary functions. Reverting to the original 
variables x, y we may also write 


w= be +y) + He —y), 


=0. 


and hence 
u(x,t) = o(x + at) + Ye — at) 
is the general solution of the wave equation of Equation (342). 
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Suppose we wish to solve the wave equation (for example, the 
vibrating string equation) with the following initial conditions ; 


u(z,0) = f(2), rc) = a), (344) 


where f(z) and g(x) are arbitrary functions. ‘This classical problem 
in the theory of partial differential equations is known as Cauchy's 
problem, With 

ufx,t) = o(x + af) + Wx — at) 
we have, at t= 0, 


(20) = fle) = He) + He). 


oe 
a 


Also 
= ap'(x + at) — ay(e — at) 
and ‘ 
Jeo) = ete) = of(2) — ap, 
‘Thus the conditions of Equations (344) require that 
w=d) +40 Boge vm 045) 
Integrating the second of the above equations between 0 and x gives 
He) — oe) = 4 flats ae + K. 
If we add this to the first of Equations (345) and divide by 2, 


He) = Spe) + fiat ar + %, (346) 
and if we subtract, and divide by 2, 
2) = he) S X fiat ax - x (47) 


Replace x by * + at in Equation (346) and by x —at in Equation 
(347) and add; then 


He + at) + We —at) = Fie + at) + fle —at)] 
+ fp ae a + z Seas) = 
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But u(x,t) = $(x + at) + Wz —at). Hence the solution of the 
Cauchy problem is 


wet) = Tite tat) + teen) + [ads 48) 
(See Problem 8 of Chapter IIL.) 


94. The homogeneous equation of hyperbolic type 
If, in Equation (331), D = E = F = 0, the equation becomes 
homogeneous* in the second derivatives, namely, 
au a a” 

Amt Bay t+ Cy =o (349) 
If B*—4AC >0, Equation (349) is hyperbolic and upon in- 
troducing the transformation 

$mAe ty arty 
where A, and A, are roots of A* + BA-+ C =0, we find 


4 =0, D=9, 
B= -AQ,—Ay, OE = 0, 
Cc =0, F =0. 
Hence Equation (349) becomes 
OY ai, 
aa 7 
whose general solution is 
4 = $E) + Wn). 


Returning to the original variables x and y, we see that the general 
solution of ose ss i 
Amt Bay + Cao 
where 
BY — 44C > 0, 


* By homogeneous we now mean that all terms appearing in the equation have 
derivatives of the same order. In particular, Equation (349) is homogencour of 
the second order. (It alto happens to be homogeneous in our previous definition.) 
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(ny) = dx +) + Wee +y), 
where ¢ and are arbitrary functions, and A, and A, are the two 
distinct roots of the algebraic equation 


AM + BX + C= 0. 


95. The parabolic equation 
If the original Equation (331) is homogeneous in the second 
derivatives and parabolic, 
au eu eu 
ARt Fae t CRA o (350) 


it reduces to 
Ou 


oe 
under the change of variables 
f=5 n=e+y, 
where A is the only real root of the algebraic equation 


AX + BA+C=0. The general solution of Equation (351) is 
therefore 





(51) 


= b(n) + Hn) 


or, reverting to x and y, the most general solution of Equation (350) 
is 


ulary) = ape + y) + Ye + ys 


where ¢ and y are arbitrary functions, and A is the only (real) 
root of the algebraic equation AN + BA + C =0. 


96. The elliptic equation 


We have seen (Section 91) that the general elliptic equation 
can be written 


au au ou 
Btete +5 +a =0. (352) 
If we let 

u(x) = w(x,y)e torr, 
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the derivatives of u are 


= 2 tes Sage torn, 
=. — fue F 


= 88 ten a 2 eden + $ we tent, 


ost 


it ze Haztoy) Swetecsn, 


aw ow oe 
a B® pteeiy pS deein 4 OP yg testers 
at b + pw z 


QF gig UP ge 


Substituting these expressions in Equation (352) yields 
aw | ok aoe 
aa t oa t («-$-4)e =o. 
This equation has been met before in connection with the study 
of Mathieu functions. 
If Equation (331) is homogeneous in the second derivatives and 
elliptic, that is, if 





au au au 
Aga t Bary + Cha = 0 
and 
Bt — 4AC <0, 
the equation reduces to 
eu 
em 7° (353) 


under the change of variables 


f=e+y 7H=ety 


where \ and are the complex conjugate roots of the algebraic 
equation AM" + BA + C = 0. Here £ and 7 are complex variables. 
‘The general solution of Equation (353) is therefore 


u(y) = $0 +9) + We +), 
where ¢ and y are arbitrary functions. 
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PROBLEMS 
Determine whether the following equations are hyperbolic, parabolic, 
or elliptic, and determine their general solution. 
fu Ow | Ow 


L gat iy tne 

2 385422 + 1 m0. 

3 ee ee 

4. tS tae te +#- 
ot 4 1g Oe 


oH +e 0 
For what values of x and y is the following equation elliptic, 
pene or hyperbolic? 
Au Ow eu 
Bat ay tS 
7. Is the following equation hyperbolic, elliptic, or parabolic ? 
ou ou ou 
im ~455 teat ty — 3xyu = 0. 
8. Solve the equation of the vibrating string for the initial conditions 
u(x,0) = 0, 
@ 
FeO) = 4 — 2). 


9. In the derivation of the canonical form of the equations of elliptic 
type it was atated that B’ 7 0, and hence that one could divide the 
equation through by this quantity. Establish, in detail, that B’ is different 
from zero. 


10. Show in detail the statement made in Section 95 that if 








aot om 
et Pay +cR- 
is parabolic, it reduces to T= 0 under the change of variables 


€=% q=et+y, 
where A = — B/24. 
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11, From the general solution of the hyperbolic equation deduce the 
solution for the parabolic case as A, —> Ag. 
[Hint : Consider the limit of 


Haz +9) — fer +9) 





as A, + Ag] 
12. Find the solution of the parabolic equation 
au au au 
fit 4ee + 4H 0 
satisfying the boundary conditions 
HOY) = $0). Fula) = He). 
13. Show that (Section 92) 


Pw 
my tw=0 


reduces to a Bemel equation if we assume a solution of the form 
w(x, y) = fey). 

14, Write explicitly the most general linear partial differential equation 
of the second order in three independent variables. 


15. Write explicitly the most general linear partial differential equation 
of the third order in two independent variables. 


16. Find the general solution of 
Sag =o 
Boe 


where the a, are constants. Investigate the case of multiple roots. 


Google 


Original from 


oats Google UNIVERSITY OF MICHIGAN 


10. 


Te 


12 
13. 


REFERENCES 


BATEMAN, H., Partial Differential Equations of Mathematical Physics, 
Dover Publications, 1944. 

Bocuwn, S. and CHANDRASEKHARAN, K., Fourier Transforms, Prince- 
ton University Press, 1949. 

Bysniy, W. E., An Elementary Treatie on Fourier's Series and 
Spherical, Cylindrical, and Ellipsoidal Harmonics, Ginn & Company, 
1893. 

‘Campseit, G. A. and Foster, R. M., Fourier Integrals for Practical 
Applications, American Telephone & Telegraph Company, Printing 
of 1942. 

(CarsLaw, H. S., Introduction to the Theory of Fourier’s Series and 
‘Integrals, Dover Publications, 3d. ed., 1930. 

Cuuncuit, R. V., Fourier Series and Boundary Value Problems, 
McGraw-Hill Book Company, Inc., 1941. 

Firs, W. B., Advanced Calculus, The Macmillan Company, 1938. 

Gurtizmin, E. A., Communication Networks, Vol. 1, John Wiley & 
Sons, Ine., 1935. 

Hopson, E. W., The Theory of Functions of a Real Variable and the 
Theory of Fourier's Series, Vol. U1, Cambridge University Prese, 
2d ed., 1926, 

Hopoman, C. D., Mathematical tables from Handbook of Chemistry 
and Physics, Chemical Rubber Publishing Company, 6th ed., 193: 

Houmise, G., Partielle Differentialgleichungen, Walter de Gruyter & 
Company, Berlin and Leipzig, 1928. 

Ince, E. L., Ordinary Differential Equations, Dover Publications, 1944. 

Inomnsott, L. R. and Zopm., O. J., An Introduction to the Mathematical 
Theory of Heat Conduction, Ginn & Company, 1913. 

Jamas, E. and Eupe, F., Tables of Functions with Formulae and Curves, 
Dover Publications, 4th ed., 1945. 

Ketzocs, 0. D., Foundations of Potential Theory, Springer, Berlin, 
1929. 








Lams, H., Hydrodynamics, Dover Publications, 1945, 
ue 


Google 


250 REFERENCES 


17. Love, A. E. H., A Treatise on the Mathematical Theory of Elasticity, 
Cambridge University Press, 3d ed., 1920. 
18, McLacutan, N. W., Complex Variable and Operational Calculus with 
Technical Applications, Cambridge University Press, 1944. 
19, McLacutan, N. W., Loud speakers: Theory, Performance, Testing 
and Design, Oxford, Clarendon Press, 1934. 
20. McLacutan, N. W., Theory and Application of Mathieu Functions, 
Oxford, Clarendon Press, 1947. 
21. Pace, L., Introduction to Theoretical Physics, D. Van Nostrand Com- 
pany, Inc., 2d ed., 1935. 
22, Parnce, B. O., A Short Table of Integrals, Ginn and Company, 3d ed., 
1929. 
23. Satvaport,M.G., Numerical Methods in Engineering, Prentice-Hall, 
Ine., 1952. 
24, Satvapont, M. G. and Muizn, K.S., The Mathematical Solution of 
Engineering Problems, McGraw-Hill Book Company, Inc., 1948. 
25. Satvapont, M.G. and Scuwarz, R.-J., Differential Equations in 
Engineering Problems, Prentice-Hall, Inc., 1954 
26. Srmarton, J. A., Electromagnetic Theory, McGraw-Hill Book Com- 
pany, Inc., Ist ed., 1941. 
27. Timosumxo, S., Theory of Elastic Stability, McGraw-Hill Book 
Company, Inc., 1st ed., 1936. 
Timosuenxo, S., Theory of Elasticity, McGraw-Hill Book Company, 
Inc., 1934. 
\V29. Timosuenxo, S., Theory of Plates and Shells, McGraw-Hill Book 
‘Company, Inc., 1st ed., 1940. 
‘Watson, G. N., A Treatise on the Theory of Bessel Functions, Cam- 
bridge University Press, 2d ed., 1944. 
31. Wy, C. R., Advanced Engineering Mathematics, McGraw-Hill Book 
‘Company, Inc., 1951. 
32, Zyomunp, A., Trigonometrical Series, Warszawa-Lw6w, 1935. 


8 


8 


Google 


Airy stress function, 36 
Beam theory, 15 
Bending of an elastic plate, 38 
Bemoulli-Euler theory, 15 
‘Bessel functions (see Bessel’s equation) 
Beasel’s equation, 179 
application to heat flow, 209 
funetions, 202 
graph of Jx(2), 204 
graph of ¥,(2), 205 
orthogonality of Bessel functions, 206 
properties of Bessel functions, 202 
solution of, 201 
zeros of Ja(x), 209 
Beta function, 226 
Biharmonic equation, 37 
Boundary conditions, 6 
nonhomogeneous, 108 
Boundary value problems, 8 
Buckling of an elastic plate, 38 
Built-in beam, 19 
Canonical forma, 236 
Cauchy equation, 112, 115, 178 
Cauchy's problem, 241 
(Characteristic functions, 223 
Characteristic numbers, 223 
Complex Fourier series, 76 
Cosine-elliptic functions, 213 
Cosine transform, 144 
Decomposition of material, 125 
Derivation of partial differential equa- 
tions, 1 
Differential’ equations (see Bessel, 
Partial) 
Differential aystem, 8 
Direct current transmission line, 27 


251 


Google 


Dirichlet conditions, 83 
Distributed parameter system, 24 
Double Fourier series, 79 
Dynamical equations, 39 


Electrostatic potential, 29 
Elliptic coordinates, 180 
Elliptic partial differential equations, 
240, 244 
Equations (10 also Bessel, Indicial, 
Legendre, Mathieu, Nonhomo- 
geneous, Ordinary. differential, 
Partial differential); 
of compatibility, 36 
in three dimensions, 38 
of continuity, 30 
of ‘elasticity, 36, 37 
derivation of, 35 
bf equilibrium, 36 
of motion of « fluid, 39 
Error function, 151 
Even functions, 56 
Exponential form of Fourier series 
(see Fourier series, complex 
form) 


Factorial, 186 
Flerunt! rigidity, 16, 38 
Flow of heat (see Heat flow) 
Force function, 28 
Forced vibrations, 125 
Fourier integral, 142, 153 
‘cosine transform, 144, 147 
examples, 155 
inverse Fourier transform, 147 
sine transform, 144, 147 
spectrum, 147, 154 
table of transforms, 160 
transform of derivative, 161 


252 


transform pairs, 155 
translation property, 159 
Fourier transform (sce Fourier integral) 
Fourier series, 45 
change of interval, 53 
example, 54 
coefficients, 47 
complex form, 76 
example, 78 
Dirichlet conditions, 83 
double, 79 
example, 82 
even and odd functions, 56 
examples, 58 
expansion in (—L, L), 67 
even and odd functions, 69 
examples, 70 
example, 68 
prolongation of functions, 72 
‘expansion in (—x, 2), 48 
‘example, 50 
formulas for expansion, 50 
full Fourier interval, 50 
half Fourier interval, 60 
mathematical conditions on, 83 
phase angles, 74 
example, 75 
points of discontinuity, 63, 
prolongation of functions, 60, 72 
examples, 62 
trial functions, 130 
trigonometric approximation, 64 
Frobenius, method of, 190 
Functions, 
Bessel, 202 
beta, 226 
cosine-elliptic, 213 
even, 56 
gamma, 187 
Legendre, 194 
Mathieu, 213 
odd, 56 
periodic, 45 
sine-elliptic, 213 
sine integral, 172 
Fundamental interval (see Fourier 
series, full Fourier interval) 


G 





dogle 


INDEX 


Gamma function, 187 
graph of, 189 
Gravitational potential, 28 


Heat equation (see Heat flow) 
Heat flow; 
Bessel’s equation, 209 
in a semi-circular plate, 112 
in an infinite bar, 142 
in one dimension, 20 
assumptions, 20 
derivation of equation of, 20 
equation of, 22 
in three dimensions, 22 
derivation of equation of, 23 
equation of, 23 
in two dimensions, 24 
solutin by separation of variables, 
101 
solution by Fourier integral, 150 
transform of heat equation, 161 
High-frequency line equations, 27 
Hooke's law, 12 
Hydrodynamics (see Equations of con- 
tinuity) 
Hyperbolic partial differential equa- 
tions, 238, 240, 243 


Impulse function, 157 

Indicial equation, 190 

Initial conditions, 6 

Inner product, 184 

Integrals, evalustion of, 149 

Inverse Fourier transform (see Fourier 
integral) 

Kirchhoft’s laws, 25 

Kronecker delta, 46 


Laminar flow, 39 
Laplacian, 23 
in polar coordinates, 113 
in spherical coordinates. 176 
Legendre functions (see Legendre’s 
equation) 
Legendre’s equation, 178 
application to potential, 199 
generating function, 197 


INDEX 


orthogonality of Legendre polyno- 
mials, 195 
polynomials, 194 
Rodrigues’ formula, 227 
solution of, 190 
Longitudinal vibrations of a bar; 
assumptions, 12 
derivation of equation of, 12 
equation of, 14 
Lumped parameter systern, 24 
Magnetic potential, 30 
Mathematical conditions; 
‘on Fourier integral, 165 
‘on Fourier series, 83 
Mathieu functions (tee Mathieu's equa- 
tion) 
Mathieu's equation, 183 
application to vibrating membrane, 
216 
orthogonality of Mathicu functions, 
215, 219 
‘solution of, 212 
Maxwell's equations, 40 
‘Mean square error, 64 
‘Method of separation of variables, 93 


Nonhomogeneous boundary conditions 
(see Boundary conditions) 


Nonhomogencous equations, 124 
solution of, 126 


Odd functions, $6 
Ordinary differential equations, 1, 7, 27, 
95, 178, 179, 183 





of Legendre polynomials, 195 

of Mathieu functions, 215, 219 
‘Orthonormal set, 185 

complete, 186 
Oscillation in 2 loud speaker horn, 40 


Parabolic partial differential equations, 
239, 244 
Partial differential equations; 
derivations, 1 


Google 


253 


elliptic, 240, 244 
hyperbolic, 238, 240, 243 
parabolic, 239, 244 
solution by Fourier integral, 
152 
solution by separation of variables, 
93 
solution by transform, 161, 163 
transformations of, 235 
Periodic functions, 45 
Phase angles, 74 
Poisson's ratio, 38 
Polar coordinates, 113 
Potential (xe also Electrostatic, Gravi- 


150, 





Power series (see Taylor series) 
Probability integral, 151 
Problems, 41, 84, 132, 166, 225, 246 


Radioactive decay, 124 
‘solution of equation of, 126 


Separation of variables, 93 
‘general case, 128 

Serica (see Fourier, Taylor) 

Simply supported beam, 18 

Sine-elliptic functions, 213 

Sine integral function, 172 

Sine transform, 144 

Spectrum, 147 

Spherical coordinates, 176 

‘Square wave, 50 

‘Step function, 157 

Stress function, 36 

Sturm-Liouville theory, 222 

‘Submarine cable equations, 27 

Superposition of solutions, 97, 111 


Taylor series, 2, 45 

‘Telegraph equations, 24 

‘Temperature distribution in a slab (see 
Heat flow) 

Tidal waves, 32 

Torsion of a prismatical bar, 39 

‘Transform (see Fourier integral) 


254 


‘Transformations of partial differential 
equations, 235 

‘Transmission line equations, 27 
derivation of equations, 25 
solution of, 105 

‘Transverse vibrations of a beam; 
assumptions, 15 
boundary conditions, 19 
derivation of equation, 15 
equation of, 18 
sokution of equation of, 117 


Google 


boundary conditions, 6 
derivation of equation of, 3 
equation of, 6 
forced vibrations, 125 
infinite in length, 152 
‘solution by Fourier integral, 
solution by separation of variables, 

4 
solution by transform, 163 
solution in closed form, 98 
Vibration of a microphone diaphragm, 
0 


Vibrations of a bar (see Longitudinal 
vibrations, Transverse vibra- 
tions of s beam) 

Voltage in « transmission line (co 
‘Transmission line equations) 





Wave equstion, 6, 12, 99 
Cauchy problem, 241 
numerical example, 99 
‘transform of, 163 


eros of Bessel function, 209 
Zonal harmonics, 194 


Original from 


oiizetn Google UNIVERSITY OF 





Original from 


oats Google UNIVERSITY OF MICHIGAN 


Original from 


oats Google UNIVERSITY OF MICHIGAN 


uM 


Original fom 


viaietn Google = uywensity oF MICHIGAN 


tid. Deo. lo, 1964 





NST. Scy, TECH. Lip, 
3115 1.5.7, BLDG, NC 


tea 
eater 








